Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model discrepancy”

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

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

At least 1,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

MFANS 2024 - Formally Proving Characteristics of Cyber-Physical Systems

Cyber-physical systems (CPS) are engineered systems that rely on the smooth integration of computational algorithms and physical elements. This integration presents new challenges for verifying that systems will behave as expected. The goal of this presentation is to present current challenges and potential solutions for the formal verification of cyber-physical systems. For cyber systems, formal methods refer to systematically rigorous mathematical techniques employed in the specification, development, analysis, and verification of both software and hardware systems. Recent advancements in computer science have yielded sophisticated tools specifically designed to address challenges associated with formal methods in complex systems. These tools leverage various foundational concepts such as logic, formal languages, program semantics, type systems, type theory, and automata theory. A notable achievement in the application of formal methods is the seL4 microkernel, claimed to be the first general-purpose operating-system kernel to be verified. Its proof implies the absence of bugs and guarantees that the kernel meets specifications. For physical systems, dynamic and control theory has a history of using rigorous analytic techniques to prove functional correctness. Lyapunov, optimal, classical, modern, and robust control theories all provide rigorous mathematical methods both to analyze system performance and to design controller that can be guaranteed to meet certain objectives. Recent computational techniques like level set theory and reachability analysis provide assertions that a system's state will avoid unsafe regions. Even though success has been independently achieved for cyber systems and physical systems, the integration of such systems creates new challenges. In particular, there is an obvious discrepancy between finite-state machines and infinite-state systems, resulting in different approaches for modeling and analyzing these system. While it is possible to simulate hybrid systems, this provides only a demonstration of a performance and not proof. For hybrid systems, current formal methods and system analysis approaches typically require a workarounds to work on hybrid systems like CPS. This paper will outline the state of the art and limits of current practice for formally verifying CPS and will identify possible research directions that require attention.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

The effectiveness of D 2 pellet and gas injection in reducing intra-ELM tungsten erosion and heat flux in the DIII-D small angle slot divertor

Edge localized modes (ELMs) in H-mode plasmas can melt and erode plasma-facing components (PFCs) and lead to impurities in the core, reducing confinement. This study analyzes the use of D 2 pellet and gas injection for ELM mitigation on the DIII-D tokamak during the 2022 Small Angle Slot V-shaped Tungsten (W) (SAS-VW) divertor campaign, reducing W erosion and heat flux during ELMs. D α (656 nm) and WI (400.9 nm) filterscopes and Langmuir probes provide photon emission and electron density/temperature to estimate W atom erosion using the S/XB method, while surface eroding thermocouples measured ELM peak heat flux at the outer strike point (OSP). Thomson Scattering measurements of pedestal T e and n e provided input to predict W divertor erosion and heat flux during ELMs via the Free-Streaming plus Recycling Model (FSRM). While a greater D 2 mass injection rate decreased the ELM peak heat flux, the impact on W erosion was not monotonic. The average ‘large’ intra-ELM W erosion was lower for plasma shots with D 2 mass injection compared to the plasma shots without any D 2 mass injection in the SAS-VW divertor. However, plasma shots did not see significant changes in ‘large’ intra-ELM W divertor erosion as the D 2 mass injection rate increased. On average, the FSRM overestimated SAS-VW experimental intra-ELM W erosion by a factor of 6.9, but the overestimation is reduced to a factor of 3.6 with the implication of C deposition effects. Generally, experimentally measured and predicted quantities were worse for plasma shots with a lower D 2 mass injection rate and corresponding higher plasma carbon (C) impurity percentage (f C ). The discrepancies are postulated to be due to C/W material mixing, for which a simple analytic mixed-material model is presented. These results highlight the importance of incorporating and improving the robustness of a mixed-material layer model in the analysis of PFC erosion on present and future tokamak devices.

Carbon (C)↗
DOE OSTI2022

Fluid – Fluid Interfacial Area and Its Impact on Relative Permeability - A Pore Network Modeling Study

Abstract Relative permeability (kr) is commonly modeled as an empirical function of phase saturation. Although current empirical models can provide a good match of one or two measured relative permeabilities using saturation alone, they are unable to predict relative permeabilities well when there is hysteresis or when physical properties such as wettability change. Further, current models often result in relative permeability discontinuities that can cause convergence and accuracy problems in simulation. To overcome these problems, recent research has modeled relative permeability as a state function of both saturation (S) and phase connectivity (X). Pore network modeling (PNM) data, however, shows small differences in relative permeability for the same S-X value when approached from a different flow direction. This paper examines the impact of one additional Minkowski parameter (Mecke and Arns, 2005), the fluid-fluid interfacial area, on relative permeability to identify if that satisfactorily explains this discrepancy. We calculate the total fluid-fluid interfacial areas (IA) during two-phase (oil/water) flow in porous media using pore network modeling. The area is calculated from PNM simulations using the areas associated with corners and throats in pore elements of different shapes. The pore network is modeled after a Bentheimer sandstone, using square, triangular, and circular pore shapes. Simulations were conducted for numerous primary drainage and imbibition cycles at a constant contact angle of 0° for the wetting phase. Simultaneous measurements of capillary pressure, relative permeability, saturation, and phase connectivity are made for each displacement. Fluid-fluid interfacial area is calculated from the PNM capillary pressure, the fluid location in the pore elements, and the pore element dimensional data. The results show that differences in the relative permeability at the same (S,X) point is explained well by differences in the fluid-fluid interfacial area (IA). That is, for a larger change in IA at these intersection points, the permeability difference is greater. That difference in relative permeability approaches zero as the difference in IA approaches zero. This confirms that relative permeability can be modeled better as a unique function of S, X, and IA. The results also show that an increase in IA restricts flow decreasing the nonwetting (oil) phase permeability. This decrease is caused by an increase in the throat area fraction compared to the corner area as the total area IA increases. The wetting phase relative permeability, however, shows the inverse trend, in that its relative permeability is greater when IA becomes larger owing to a greater fraction of the total area associated with the corners. The area IA, however, impacts the nonwetting phase relative permeability more than the wetting phase relative permeability. Corner flow improves the wetting phase relative permeability because the wetting phase is continuous there. Finally, a sensitivity analysis shows that relative permeability a is more sensitive to change in S than they are for IA for the case studied implying that if only two parameters are used to model relative permeability it is better to choose S and X.

Engineering↗
DOE OSTI2022

Improved $^{95}\mathrm{Mo}$ neutron resonance parameters and astrophysical reaction rates

We report improved 95 Mo neutron resonance parameters and reaction rates are important for nuclear astrophysics, testing nuclear models, and nuclear criticality safety. However, despite many previous neutron-capture and total cross-section measurements on this nuclide, there still is much room for improvement as well as several discrepancies. For example, there are very few firm resonance spin and parity assignments; average resonance parameters are available only for each parity, the currently recommended astrophysical reaction rate results in disagreements between stellar models and meteoric isotopic anomalies, and there are substantial disagreements in the neutron-capture cross section at low energies important for nuclear criticality safety. To obtain an improved set of neutron resonance parameters and astrophysical reaction rates for 95 Mo. High-resolution neutron-capture and transmission data were measured at the Oak Ridge Electron Linear Accelerator (ORELA) using highly isotopically enriched 95 Mo samples. The neutron-capture apparatus, data reduction, and analysis were improved so that information contained in the γ-ray cascade following neutron capture were used to assign resonance J π values. Following this, simultaneous analysis of the new neutron-capture and transmission data was used to obtain resonance energies, gamma widths, and neutron widths and their uncertainties to a maximum energy of 10 keV. Accurate neutron-capture cross sections also were obtained for the unresolved resonance region to a maximum energy of 500 keV and, together with the new resonance parameters, used to calculate the astrophysical reaction rates in the temperature range from 5 to 30 keV. A vastly improved set of 95Mo neutron resonance parameters and an astrophysical reaction rate accurate to about 3% were obtained. Firm J π assignments were determined for 261 of the 314 observed resonances. This is a very large improvement over the previously published 32 firm J π assignments for 108 resonances. Also, the number of resonances having both firm J π assignments and Γ γ values was increased by almost a factor of 24—from 11 to 261. Neutron- and total-radiation-width distributions and average resonance spacings, average total radiation widths, and neutron strength functions were obtained for the six different s- and p-wave possibilities. Parameters for the lowest s-wave resonance, which is most important for criticality benchmarks, were obtained with high accuracy. Simple modification of the neutron-capture apparatus and expansion and improvement of data analysis techniques led to a large increase in firm J π assignments for 95 Mo neutron resonances. The resulting astrophysical reaction rate is 20%–30% larger than the currently recommended rate at s-process temperatures, which should lead to better agreement between stellar models and meteoric isotopic anomalies. The neutron-capture cross section at low energies is substantially larger than recommended in the latest evaluation, which is problematical for criticality benchmarks. The average resonance spacing as a function of spin and parity is significantly different from current models. The total-radiation-width distributions are significantly broader than predicted by theory and show significant departures from the expected Gaussian shapes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Meso-modeling of Carbon Fiber Composite for Crash Safety Analysis

In the conventional approach, the material properties for crash safety simulations are typically obtained from standard coupon tests, where the test results only provide single layer material properties used in crash simulations. However, the lay-up effects for the failure behaviors of the real structure were not considered in numerical simulations. Hence, there was discrepancy between the crash simulations and experimental tests. Consequently, an intermediate stage is required for accurate predictions. Some component tests are required to correlate the material models in the intermediate stage. In this paper, a Mazda Tube under high-impact velocity is chosen as an example for the crash safety analysis. The tube consists of 24 layers of uni-directional (UD) carbon fiber composite materials, in which 4 layers are perpendicular to, while the other layers are parallel to the impact direction. An LS-DYNA meso-model was constructed with orthotropic material models counting for the single-layer material behaviors. Between layers, a node-based tie-break contact was used for modeling the delamination of the composite material. Since fiber directions are not single-oriented, the lay-up effects could be an important effect. From the first numerical trial, premature material failure occurred due to the use of material parameters obtained directly from the coupon tests. Some parametric studies were conducted to identify the cause of the numerical instability. The finding is that the material failure strength used in the numerical model needs to be enlarged to stabilize the numerical model. Some hypothesis was made to provide the foundation for enlarging the failure strength and the corresponding experiments will be conducted to validate the hypothesis.

Lin, Shih-Po↗
NASA NTRS1923

Preliminary Experiments to Determine Scale and Slip-Stream Effects on a 1/24th Size Model of a JN4H Biplane

This work was undertaken to obtain results on a small model of a complete airplane which might be used for comparison with corresponding tests made in full flight. Somewhat similar tests have been previously made at various other laboratories; but as certain discrepancies exist between corresponding tests in different tunnels, it has been deemed advisable to obtain a direct comparison for this particular installation. The present work covers tests on a one-twenty-fourth scale model at speeds varying from 6.7 m/sec. (15 m.p.h.) to 40.2 m/sec, (90 m.p.h.). A slip stream correction has been obtained by the use of a small belt-driven propeller mounted in front of the model, and force coefficients thus obtained are compared with the measurements of the same forces made in full flight on a geometrically similar airplane. This report gives lift, drag, and longitudinal moment values obtained in tests of a particularly accurate model over a wide range of speeds. A measure of the slip stream corrections on lift and drag forces was obtained by the use of a power-driven model propeller. Measurements were also made of forces and longitudinal moments for all angles from 0 degree to 360 degrees.

Bacon, D L↗
DOE OSTI2024

Representing cropping systems with the MEMS 2 ecosystem model

Abstract Croplands have been the focus of substantial investigation due to their considerable potential for sequestering carbon. Understanding the potential for soil organic carbon (SOC) sequestration and necessary management strategies will be enabled with accurate process‐based models. Accurately representing crop growth and agricultural practices will be critical for realistic SOC modeling. The MEMS 2 model incorporates a current understanding of SOC formation and stabilization, measurable SOC pools, and deep SOC dynamics and is seen as a highly promising tool to inform management intervention for SOC sequestration. Thus far, MEMS 2 has been developed to represent grasslands. In this study, we further developed MEMS 2 to model annual grain crops and common agricultural practices, such as irrigation, fertilization, harvesting, and tillage. Using four Ameriflux sites, we demonstrated an accurate simulation of crop growth and development. Model performance was strong for simulating aboveground biomass (index of agreement [d] range of 0.89–0.98) and green leaf area index (dfrom 0.90 to 0.96) across corn, soybean, and winter wheat. Good agreement with observations was also achieved for net ecosystem CO 2 exchange (dfrom 0.90 to 0.96), evapotranspiration (dfrom 0.91 to 0.94), and soil temperature (dof 0.96), while discrepancy with the available soil water content data remain (dfrom 0.14 to 0.81 at four depths to 100 cm). While we will continue model testing and improvement, MEMS 2 (version 2.14) has now demonstrated its ability to effectively simulate the growth of common grain crops and practices.

Agriculture↗
NASA NTRS1987

A survey of nested grid techniques and their potential for use within the MASS weather prediction model

A survey of various one- and two-way interactive nested grid techniques used in hydrostatic numerical weather prediction models is presented and the advantages and disadvantages of each method are discussed. The techniques for specifying the lateral boundary conditions for each nested grid scheme are described in detail. Averaging and interpolation techniques used when applying the coarse mesh grid (CMG) and fine mesh grid (FMG) interface conditions during two-way nesting are discussed separately. The survey shows that errors are commonly generated at the boundary between the CMG and FMG due to boundary formulation or specification discrepancies. Methods used to control this noise include application of smoothers, enhanced diffusion, or damping-type time integration schemes to model variables. The results from this survey provide the information needed to decide which one-way and two-way nested grid schemes merit future testing with the Mesoscale Atmospheric Simulation System (MASS) model. An analytically specified baroclinic wave will be used to conduct systematic tests of the chosen schemes since this will allow for objective determination of the interfacial noise in the kind of meteorological setting for which MASS is designed. Sample diagnostic plots from initial tests using the analytic wave are presented to illustrate how the model-generated noise is ascertained. These plots will be used to compare the accuracy of the various nesting schemes when incorporated into the MASS model.

Koch, Steven E.↗
NASA NTRSDate not supplied

A Comprehensive Forward Model for Spaceborne Radar Instruments

We present the development and validation of a comprehensive forward model designed to enhance remote sensing capabilities of spaceborne radar instruments. To overcome limitations in existing models, we integrated a Discrete Dipole Approximation (DDA) cloud scattering database into our Radiative Transfer Model (RTM), focusing on microwave frequencies. By simulating the optical properties of non-spherical frozen hydrometeors, the DDA technique effectively reduced discrepancies between simulated and observed values, surpassing traditional Mie tables. The evaluation of DDA lookup tables involved comparisons with a collocated dataset comprising short-term forecasts and satellite microwave data, providing evidence of their superiority. Additionally, we address the challenges of assimilating active radar measurements, which offer vertically resolved insights into clouds and precipitation. We explored the assimilation of spaceborne radar measurements in Numerical Weather Prediction (NWP) models by integrating a forward radar model, along with its adjoint and tangent linear, into the data assimilation system. Evaluation using CloudSat measurements demonstrated promising agreement between simulations and observations, particularly when the input hydrometeor profiles aligned with the measured reflectivity profiles, showcasing the potential of the developed forward radar model. Moreover, we discuss other challenges in radar measurement assimilation within NWP models, including potential observation errors and biases.

Isaac Moradi↗
DOE OSTI2025

Towards a real-time computation of timelike hadronic vacuum polarization and light-by-light scattering: Schwinger Model tests

Hadronic vacuum polarization (HVP) and light-by-light scattering (HLBL) are crucial for evaluating the Standard Model predictions concerning the muon’s anomalous magnetic moment. However, direct first-principle lattice gauge theory-based calculations of these observables in the timelike region remain challenging. Discrepancies persist between lattice quantum chromodynamics (QCD) calculations in the spacelike region and dispersive approaches relying on experimental data parametrization from the timelike region. Here, we introduce a methodology employing 1+1-dimensional quantum electrodynamics (QED), i.e. the Schwinger Model, to investigate the HVP and HLBL. To that end, we use both tensor network techniques, specifically matrix product states, and classical emulators of digital quantum computers. Demonstrating feasibility in a simplified model, our approach sets the stage for future endeavors leveraging digital quantum computers.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2018

(GO)2-SIM: a GCM-Oriented Ground-Observation Forward-Simulator Framework for Objective Evaluation of Cloud and Precipitation Phase

General circulation model (GCM) evaluation using ground-based observations is complicated by inconsistencies in hydrometeor and phase definitions. Here we describe (GO)2-SIM, a forward simulator designed for objective hydrometeor-phase evaluation, and assess its performance over the North Slope of Alaska using a 1-year GCM simulation. For uncertainty assessment, 18 empirical relationships are used to convert model grid-average hydrometeor (liquid and ice, cloud, and precipitation) water contents to zenith polarimetric micropulse lidar and Ka-band Doppler radar measurements, producing an ensemble of 576 forward-simulation realizations. Sensor limitations are represented in forward space to objectively remove from consideration model grid cells with undetectable hydrometeor mixing ratios, some of which may correspond to numerical noise.Phase classification in forward space is complicated by the inability of sensors to measure ice and liquid signals distinctly. However, signatures exist in lidar–radar space such that thresholds on observables can be objectively estimated and related to hydrometeor phase. The proposed phase-classification technique leads to misclassification in fewer than 8% of hydrometeor-containing grid cells. Such misclassifications arise because, while the radar is capable of detecting mixed-phase conditions, it can mistake water- for ice-dominated layers. However, applying the same classification algorithm to forward-simulated and observed fields should generate hydrometeor-phase statistics with similar uncertainty. Alternatively, choosing to disregard how sensors define hydrometeor phase leads to frequency of occurrence discrepancies of up to 40%. So, while hydrometeor-phase maps determined in forward space are very different from model "reality" they capture the information sensors can provide and thereby enable objective model evaluation.

Lamer, K.↗
DOE OSTI2018

(GO) 2 -SIM: a GCM-oriented ground-observation forward-simulator framework for objective evaluation of cloud and precipitation phase

General circulation model (GCM) evaluation using ground-based observations is complicated by inconsistencies in hydrometeor and phase definitions. Here we describe (GO) 2 -SIM, a forward simulator designed for objective hydrometeor-phase evaluation, and assess its performance over the North Slope of Alaska using a 1-year GCM simulation. For uncertainty assessment, 18 empirical relationships are used to convert model grid-average hydrometeor (liquid and ice, cloud, and precipitation) water contents to zenith polarimetric micropulse lidar and Ka-band Doppler radar measurements, producing an ensemble of 576 forward-simulation realizations. Sensor limitations are represented in forward space to objectively remove from consideration model grid cells with undetectable hydrometeor mixing ratios, some of which may correspond to numerical noise. Phase classification in forward space is complicated by the inability of sensors to measure ice and liquid signals distinctly. However, signatures exist in lidar–radar space such that thresholds on observables can be objectively estimated and related to hydrometeor phase. The proposed phase-classification technique leads to misclassification in fewer than 8 % of hydrometeor-containing grid cells. Such misclassifications arise because, while the radar is capable of detecting mixed-phase conditions, it can mistake water- for ice-dominated layers. However, applying the same classification algorithm to forward-simulated and observed fields should generate hydrometeor-phase statistics with similar uncertainty. Alternatively, choosing to disregard how sensors define hydrometeor phase leads to frequency of occurrence discrepancies of up to 40 %. So, while hydrometeor-phase maps determined in forward space are very different from model reality they capture the information sensors can provide and thereby enable objective model evaluation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Time-sequence Spectroscopy of Epsilon CrA: The 518 nm Mg i Triplet Region Analyzed with Broadening Functions

High-resolution spectroscopic observations of the W UMa-type binary ϵ CrA obtained as a time-monitoring sequence on four full and four partial nights within two weeks have been used to derive orbital elements of the system and discuss the validity of the Lucy model for description of the radial-velocity data. The observations had more extensive temporal coverage and better quality than similar time-sequence observations of the contact binary AW UMa. The two binaries share several physical properties and show very similar deviations from the Lucy model: the primary component is a rapidly rotating star almost unaffected by the presence of the secondary component, while the latter is embedded in a complex gas flow and appears to have its own rotation-velocity field, in contradiction to the model. The spectroscopic mass ratio is found to be larger than the one derived from the light-curve analysis, as in many other W UMa-type binaries, but the discrepancy for ϵ CrA is relatively minor, much smaller than for AW UMa. The presence of the complex velocity flows contradicting the assumption of solid-body rotation suggests a necessity of modification to the Lucy model, possibly along the lines outlined by Stȩpień in his concept of energy transfer between the binary components.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Solving inverse problems in stochastic models using deep neural networks and adversarial training

Inverse problems associated with stochastic models constitute a significant portion of scientific and engineering applications. In such cases the unknown quantities are distributions. The applicability of traditional methods is limited because of their demanding assumptions or prohibitive computational consumption; for example, maximum likelihood methods require closed-form density functions, and Markov Chain Monte Carlo needs a large number of simulations. We propose a new method that estimates the unknown distribution by matching the statistical properties between observed and simulated random processes. We leverage the expressive power of neural networks to approximate the unknown distribution and use a discriminative neural network for computing the statistical discrepancies between the observed and simulated random processes. Here we demonstrated numerically that the proposed methods can estimate both the model parameters and learn complicated unknown distributions.

42 ENGINEERING↗
DOE OSTI2017

Multipoint Observations of Energetic Particle Injections and Substorm Activity During a Conjunction Between Magnetospheric Multiscale (MMS) and Van Allen Probes

Here, this study examines multipoint observations during a conjunction between Magnetospheric Multiscale (MMS) and Van Allen Probes on 7 April 2016 in which a series of energetic particle injections occurred. With complementary data from Time History of Events and Macroscale Interactions during Substorms, Geotail, and Los Alamos National Laboratory spacecraft in geosynchronous orbit (16 spacecraft in total), we develop new insights on the nature of energetic particle injections associated with substorm activity. Despite this case involving only weak substorm activity (maximum AE <300 nT) during quiet geomagnetic conditions in steady, below-average solar wind, a complex series of at least six different electron injections was observed throughout the system. Intriguingly, only one corresponding ion injection was clearly observed. All ion and electron injections were observed at <600 keV only. MMS reveals detailed substructure within the largest electron injection. A relationship between injected electrons with energy <60 keV and enhanced whistler mode chorus wave activity is also established from Van Allen Probes and MMS. Drift mapping using a simplified magnetic field model provides estimates of the dispersionless injection boundary locations as a function of universal time, magnetic local time, and L shell. The analysis reveals that at least five electron injections, which were localized in magnetic local time, preceded a larger injection of both electrons and ions across nearly the entire nightside of the magnetosphere near geosynchronous orbit. The larger ion and electron injection did not penetrate to L < 6.6, but several of the smaller electron injections penetrated to L < 6.6. Due to the discrepancy between the number, penetration depth, and complexity of electron versus ion injections, this event presents challenges to the current conceptual models of energetic particle injections.

58 GEOSCIENCES↗
NASA NTRS1984

A redetermination of the succinonitrile-water phase diagram

The phase diagram of the water/succinonitrile system often used in modeling binary metallic alloys is determined experimentally using a laser scattering method. Triply distilled and zone-refined succinonitrile is placed in a jacketed water-cooled vessel with a magnetic mixer and a copper-constantan thermocouple, and the melting point is determined under vacuum and under Ar by observing the response of a photodiode placed opposite a laser source as the chamber is cooled and reheated. Then water is added, and the eutectic and monotectic temperatures and the entire phase diagram are determined under Ar. The results are presented in a table and graph and compared to those of Schreinemakers (1897): significant discrepancies are noted and attributed to impurities in the succinonitrile used by Schreinemakers. The importance of accurate phase diagrams for modeling binary alloys is stressed.

Smith, J. E., Jr.↗
NASA NTRS2003

Measurement and Modeling of Vertically Resolved Aerosol Optical Properties and Radiative Fluxes Over the ARM SGP Site

In order to meet one of its goals - to relate observations of radiative fluxes and radiances to the atmospheric composition - the Department of Energy's Atmospheric Radiation Measurement (ARM) program has pursued measurements and modeling activities that attempt to determine how aerosols impact atmospheric radiative transfer, both directly and indirectly. However, significant discrepancies between aerosol properties measured in situ or remotely remain. One of the objectives of the Aerosol Intensive Operational Period (TOP) conducted by ARM in May 2003 at the ARM Southern Great Plains (SGP) site in north central Oklahoma was to examine and hopefully reduce these differences. The IOP involved airborne measurements from two airplanes over the heavily instrumented SGP site. We give an overview of airborne results obtained aboard the Center for Interdisciplinary Remotely-Piloted Aircraft Studies (CIRPAS) Twin Otter aircraft. The Twin Otter performed 16 research flights over the SGP site. The aircraft carried instrumentation to perform in-situ measurements of aerosol absorption, scattering, extinction and particle size. This included such novel techniques as the photoacoustic and cavity ring-down methods for in-situ absorption (675 nm) and extinction (675 and 1550 nm) and a new multiwavelength, filter-based absorption photometer (467, 530, 660 nm). A newly developed instrument measured cloud condensation nucleus concentration (CCN) concentrations at two supersaturation levels. Aerosol optical depth and extinction (354-2139 nm) were measured with the NASA Ames Airborne Tracking 14-channel sunphotometer. Furthermore, up-and downwelling solar (broadband and spectral) and infrared radiation were measured using seven individual radiometers. Three up-looking radiometers werer mounted on a newly developed stabilized platform, keeping the instruments level up to aircraft pitch and roll angles of approximately 10(exp 0). This resulted in unprecedented continuous vertical profiles of radiative fluxes, which we will compare to modeled fluxes using the aforementioned data as input.

Schmid, B.↗
DOE OSTI2024

Relating flow resistance to equivalent roughness

Describing flow resistance using the physical properties of an underlying surface is a recalcitrant problem in overland flow models. If discharge measurements are available, an equivalent roughness (e.g., Manning’s n) can be calibrated to represent the effects of surface properties within the domain with a single numerical value. Alternatively, the flow resistance can be estimated from discharge and velocity measured at a point, typically a runoff plot outlet. However, such experimental estimates are often inconsistent with the equivalent roughness determined from calibration to discharge, even if both derive from the same dataset. For example, if Manning’s equation is used to parameterize flow resistance, the Manning’s n obtained by calibrating a model to discharge differs from the value of n calculated from measured flow and velocity at the hillslope outlet. Here, this discrepancy is resolved by deriving a correction factor relating experimentally-determined flow resistance to the equivalent roughness. The derived correction factor is tested for four commonly-used resistance formulations using 129 rainfall simulator experiments. The correction factor is necessary to reproduce measured velocities, and yields minor improvements in discharge prediction. Plain Language Summary: Accurate runoff prediction is needed for land and water management in dryland regions, where sporadic and limited rainfall necessitate efficient water use and drought mitigation strategies. The skill of runoff models is known to be hindered by out ability to estimate flow resistance, which is the quantity that describes how energy is lost from flowing water to the underlying surface. Typically, models represent flow resistance with an equivalent roughness, e.g., Manning’s n, that is adjusted until the model can reproduce available discharge observations at watershed scale. However, the flow resistance measured in plot-scale experiments (1–10 m) often exceeds equivalent roughness coefficients by a factor of 10. This means that the direct use of plot-scale experimental data to parameterize runoff models could cause errors in discharge and runoff velocity predictions. Here, we resolve these differences by deriving an analytic correction factor that relates flow resistance to the equivalent roughness required for models to reproduce experimental velocity and discharge data. This correction factor is tested using rainfall simulator data from 129 experiments performed in the US Southwest covering a wide range of precipitation intensities, soil textures and vegetation types. Use of the correction factor substantially improves model prediction of flow velocity, which is needed for reproducing the timing of flood events and the estimation of erosion.

54 ENVIRONMENTAL SCIENCES↗