Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Profiling”

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 73 records · Page 4

Data-driven analysis of neutron diffraction line profiles: application to plastically deformed Ta

Abstract Non-destructive evaluation of plastically deformed metals, particularly diffraction line profile analysis (DLPA), is valuable both to estimate dislocation densities and arrangements and to validate microstructure-aware constitutive models. To date, the interpretation of whole line diffraction profiles relies on the use of semi-analytical models such as the extended convolutional multiple whole profile (eCMWP) method. This study introduces and validates two data-driven DLPA models to extract dislocation densities from experimentally gathered whole line diffraction profiles. Using two distinct virtual diffraction models accounting for both strain and instrument induced broadening, a database of virtual diffraction whole line profiles of Ta single crystals is generated using discrete dislocation dynamics. The databases are mined to create Gaussian process regression-based surrogate models, allowing dislocation densities to be extracted from experimental profiles. The method is validated against 11 experimentally gathered whole line diffraction profiles from plastically deformed Ta polycrystals. The newly proposed model predicts dislocation densities consistent with estimates from eCMWP. Advantageously, this data driven LPA model can distinguish broadening originating from the instrument and from the dislocation content even at low dislocation densities. Finally, the data-driven model is used to explore the effect of heterogeneous dislocation densities in microstructures containing grains, which may lead to more accurate data-driven predictions of dislocation density in plastically deformed polycrystals.

36 MATERIALS SCIENCE↗

TDCOSMO - IV. Hierarchical time-delay cosmography – joint inference of the Hubble constant and galaxy density profiles

The H0LiCOW collaboration inferred via strong gravitational lensing time delays a Hubble constant value of H0 = 73.3−1.8+1.7 km s−1 Mpc−1, describing deflector mass density profiles by either a power-law or stars (constant mass-to-light ratio) plus standard dark matter halos. The mass-sheet transform (MST) that leaves the lensing observables unchanged is considered the dominant source of residual uncertainty in H0. We quantify any potential effect of the MST with a flexible family of mass models, which directly encodes it, and they are hence maximally degenerate with H0. Our calculation is based on a new hierarchical Bayesian approach in which the MST is only constrained by stellar kinematics. The approach is validated on mock lenses, which are generated from hydrodynamic simulations. We first applied the inference to the TDCOSMO sample of seven lenses, six of which are from H0LiCOW, and measured H0 = 74.5−6.1+5.6 km s−1 Mpc−1. Secondly, in order to further constrain the deflector mass density profiles, we added imaging and spectroscopy for a set of 33 strong gravitational lenses from the Sloan Lens ACS (SLACS) sample. For nine of the 33 SLAC lenses, we used resolved kinematics to constrain the stellar anisotropy. From the joint hierarchical analysis of the TDCOSMO+SLACS sample, we measured H0 = 67.4−3.2+4.1 km s−1 Mpc−1. This measurement assumes that the TDCOSMO and SLACS galaxies are drawn from the same parent population. The blind H0LiCOW, TDCOSMO-only and TDCOSMO+SLACS analyses are in mutual statistical agreement. The TDCOSMO+SLACS analysis prefers marginally shallower mass profiles than H0LiCOW or TDCOSMO-only. Without relying on the form of the mass density profile used by H0LiCOW, we achieve a ∼5% measurement of H0. While our new hierarchical analysis does not statistically invalidate the mass profile assumptions by H0LiCOW – and thus the H0 measurement relying on them – it demonstrates the importance of understanding the mass density profile of elliptical galaxies. The uncertainties on H0 derived in this paper can be reduced by physical or observational priors on the form of the mass profile, or by additional data.Key words: gravitational lensing: strong / galaxies: general / galaxies: kinematics and dynamics / distance scale / cosmological parameters / cosmology: observations⋆ The full analysis is available at https://github.com/TDCOSMO/hierarchy_analysis_2020_public.

79 ASTRONOMY AND ASTROPHYSICS↗

Probing nonlinear velocity profiles of shear-thinning, nematic platelet dispersions in Couette flow using x-ray photon correlation spectroscopy

In this study, we report experiments employing x-ray photon correlation spectroscopy (XPCS) to characterize the velocity profiles of complex fluids in Couette flow. The approach involves modeling the XPCS correlation functions obtained with the incident x-ray beam passing tangentially through the Couette cell gap at various distances from the inner wall. We first demonstrate the technique with measurements on a dilute colloidal dispersion in the Newtonian liquid glycerol, where the expected linear velocity profiles are recovered. We then employ the technique to map the shear-rate-dependent velocity profiles of a shear-thinning dispersion of nematically ordered Gibbsite platelets. The nonlinear velocity profiles of the Gibbsite dispersion include a narrow slip region adjacent to the outer wall and a band with a small velocity gradient in the interior of the gap that evolves into a region increasingly resembling plug flow with increasing shear rate. Variations in the velocity profile along the vorticity direction indicate an instability in the interface between this region of small velocity gradient and a region of high velocity gradient near the inner wall. The analysis of the small-angle scattering patterns provides information about the spatial and temporal variations in the nematic order of the Gibbsite dispersion and their coupling to the velocity profile. Additional potential applications of this XPCS-based technique and comparisons with established methods for characterizing velocity profiles are discussed.

42 ENGINEERING↗

Equilibrium reconstruction of DIII-D plasmas using predictive modeling of the pressure profile

New workflows have been developed for predictive modeling of magnetohydrodynamic (MHD) equilibrium in tokamak plasmas. The goal of this work is to predict the MHD equilibrium in tokamak discharges without having measurements of the kinetic profiles. The workflows include a cold start tool, which constructs all the profiles and power flows needed by transport codes; a Grad–Shafranov equilibrium solver; and various codes for the sources and sinks. For validation purposes, a database of DIII-D tokamak discharges has been constructed that is comprised of scans in the plasma current, toroidal magnetic field, and triangularity. Initial efforts focused on developing a workflow utilizing an empirically derived pressure model tuned to DIII-D discharges with monotonic safety factor profiles. This workflow shows good agreement with experimental kinetic equilibrium calculations, but is limited in that it is a single fluid (equal ion and electron temperatures) model and lacks H-mode pedestal predictions. The best agreement with the H-mode database is obtained using a theory-based workflow utilizing pressure profile predictions from a coupled TGLF turbulent transport and EPED pedestal models together with external magnetics and Motional Stark Effect (MSE) data to construct the equilibrium. Here, we obtain an average root mean square error of 5.1% in the safety factor profile when comparing the predicted and experimental kinetic equilibrium. We also find good agreement with the plasma stored energy, internal inductance, and pressure profiles. Including MSE data in the theory-based workflow results in noticeably improved agreement with the q-profiles in high triangularity discharges in comparison with the results obtained with magnetic data only. The predictive equilibrium workflow is expected to have wide applications in experimental planning, between-shot analysis, and reactor studies.

Kinsey, J. E. (ORCID:0000000193347473)↗

Reconstruction of current profiles from external magnetic and MSE measurements on DIII-D

This work introduces a new current model into the EFIT code for plasma current profile reconstruction. Using this model, edge-localized current density profiles are reconstructed solely from external magnetic measurements. Quantitative comparisons show that edge current profiles obtained with the new current model align closely with those from kinetic equilibrium reconstruction. The new model enables current density reconstruction with fewer diagnostics than kinetic equilibrium reconstruction (Lao et al 1990 Nucl. Fusion 30 1035–49), demonstrating that edge-localized current density information can be obtained from external magnetic measurements only, without the help of kinetic profiles and internal diagnostics, and without using the bootstrap current model. Furthermore, with the constraints of external magnetic and MSE measurements, current profiles in full radius obtained from the new model are in good agreement with the results of kinetic equilibrium reconstruction, which not only improves the efficiency of current reconstruction but also helps avoid possible errors caused by different kinds of diagnoses. Possible applications are also discussed, and a brief analysis of the evolution of the full-radius current density is performed. Noticeable decreases in both the magnitude and gradient of the edge current around the time of ELM bursts are observed. Moreover, analysis of the full radius current density profiles ( ρ N = 0 − 1.0, ρN denotes to square root of the normalized toroidal flux) reveals that as the edge current changes, the inner current also changes following a sinusoidal-like pattern at almost all time slices during four type-I ELMs in discharge 161493. The current evolution in the inner region ( ρ N = 0.2 − 0.6) is synchronized with the boundary region ( ρ N = 0.89 − 1.0), while the intermediate region ( ρ N = 0.6 − 0.89) shows opposite behavior. This sinusoidal-like variations in full-radius current density profile indicates strong core-pedestal coupling, suggesting ELMs significantly affect core plasmas, though current changes manifest differently across radial regions.

EFIT↗

Attached Jet Velocity Profiles in Mixing Tanks

Free jets have been studied in detail over much of the last century, but the theory for offset and attached jets remains incomplete. Attached jets differ from free jets in that they lose momentum to nearby surfaces, attenuating their velocities. The velocity profiles of free circular jets are nearly Gaussian, with quantitative mathematical descriptions derived from first principles by Goertler and Tollmien (Rajaratnam, 1976). In contrast, mathematical descriptions of three-dimensional attached jets from circular nozzles remain much less mature. Agelin-Chaab and Tachie (2011) used particle imaging velocimetry of a three-dimensional attached jet to show that the scaled velocity decays with scaled distance from the nozzle with a power law exponent between -1.15 and -1.21, which is larger in magnitude than that of a free jet. However, quantitative analytical expressions for the velocity profiles of attached jets similar to those of free jets remain elusive. This paper addresses this critical gap. Here we evaluate the velocity profiles of three-dimensional offset jets emerging from circular nozzles that that become attached jets. These jets lose momentum due to interactions with nearby surfaces and are important to evaluating flows in mixing vessels and to suspending solids and trapped gases in radioactive waste tanks. Despite the importance of attached jets, prior insight has been purely experimental, limited to overly simplistic analytical models, or restricted to computationally expensive computational fluid dynamics case studies. We compare the expression of Verhoff (1963) to experimental results to find reasonable quantitative agreement. As stated by Agelin-Chaab and Tachie (2011), “detailed velocity measurements of 3D offset jets are rare.” Such remains the case. This study adds to the literature by providing information at two additional Reynolds numbers (1.43 106 and 1.87 106) and evaluating simple but accurate expressions for velocity profiles. These Reynolds numbers are an order of magnitude higher than prior reports. The semi-empirical stream wise velocity profile perpendicular to the surface proposed by Verhoff (1963) is in approximate agreement with these velocity profiles, which is surprising because these attached jets are three-dimensional instead of two-dimensional as evaluated by Verhoff. However, additional work is necessary to fully describe these profiles quantitatively.

jet, attached jet, offset jet↗

Experimental Study of TCP Throughput Profiles and Dynamics Over Dedicated Connections

A throughput profile expressed as a function of the round-trip times of a data transport infrastructure is a critical indicator of its level of optimization, particularly, over dedicated connections. We study the throughput profiles of eleven TCP versions using measurements collected over dedicated hardware-emulated connections with distances spanning round the earth. Among them, BBR has overall higher temporal variations and lower throughput profiles compared to five loss-based TCP versions, and the comparison with other TCP versions is mixed. Using the Poincare map regions of throughput time traces, we show that the richer dynamics of BBR are correlated with its lower throughput profiles. We present basic analytical results that indicate that higher temporal variations of TCP methods lead to lower throughput profiles, thereby establishing the underlying causality. Under external losses, BBR achieves sustained high throughput compared to others as the loss rate is increased. Overall, our results provide insights into the relationship between time dynamics and throughput profiles of TCP versions, in particular, show that limited time traces can be indicative of global properties of throughput profiles.

Rao, Nageswara↗

(U) A Linear Response Model Predicts Reactivity From a Density Profile

We tested the ability to predict the system reactivity, described by alpha, given a density profile using a simple linear system response. We generated a suite of 1-dimensional density profiles that consisted of nominal density, a discontinuity, and a decay. These profiles were prescribed a functional form and the mass was conserved in all cases. From these density profiles, we calculated the alpha value of the 3-dimensional system.We calculated a linear response function given a training set of the 1-dimensional density profiles, and the system reactivity described by alpha. We tested the robustness of the response function using the remaining test data. Our results showed very good agreement between the predicted and calculated test values, where the distribution of alpha differences was centered about zero and had a standard deviation of 0.005 gens/shake. The predicted and calculated alpha values did not significantly differ (t=-0.0009 p<0.99). We used Singular Value Decomposition (SVD) to reduce the matrix rank by retaining95% of the cumulative singular value contributions. This reduced the matrix rank by 91.7%. We generated the linear response matrix and calculated the difference between the predicted and calculated alpha values. Using the reduced order matrix, we showed good agreement between the predicted and calculated alpha values where the distribution of differences was centered near zero, the standard deviation was 0.006 gens/shake, and the statistical t-test showed good agreement (t=0.02, p<0.98). These results show a linear relationship between a series of 1-dimensional density profiles,where the mass was conserved, and the system reactivity. The next steps of this work will be to investigate the linear response using 2-dimensional density profiles.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Demonstrating Load-Shaping Capabilities of Cost Minimizing Heat Pump Water Heater Controls with Varying Price Profiles

Increased penetration of photovoltaics and electrification of traditionally gas appliances are exacerbating existing challenges in cost-effectively balancing electricity grid supply and demand. Decarbonization without incurring expensive transmission and distribution system capacity increases requires shifting building loads from peak demand times to peak renewable production times. Grid operators are evaluating new ways of encouraging load shifting, including using time-varying price structures to provide a financial incentive. If devices incorporate price-responsive controls, a price profile could be designed to yield a wide variety of load curves as needed to optimize grid functionality. Heat pump water heaters (HPWHs) are an ideal device for price-responsive controls because the storage tank enables them to optimize the timing of electricity consumption without impacting hot water delivery service. This paper presents work demonstrating how price-responsive controls for HPWHs can provide different load profiles, as needed to stabilize the grid, in response to different price profiles. HPWH manufacturers now include web API and CTA-2045 communication capabilities which enable sending load shaping control signals. Pilot studies and preliminary programs have utilized these capabilities with uniform control strategies to reduce 4-9 PM electricity consumption. However, no studies have developed flexible controls capable of both a) responding to constantly varying price profiles and b) customizing logic to match the needs of each HPWH. Berkeley Lab's CalFlexHub project is pioneering price-driven load flexibility by developing and deploying cost-minimizing controls utilizing setpoint setting signals for fleets of HPWHs in response to varying price profiles. Control development is based on simulations using the Flexible Heat Pump Water Heater Performance Predictor which captures the control decisions of a residential, integrated HPWH manufacturer’s on-board controller. The proposed cost-reducing controls respond to constantly changing price profiles, providing the ability to change the price profile to generate load curves as needed to maintain grid stability. Preliminary simulations studying the load shaping capabilities of price-responsive controls on a fleet of 60 HPWHs have demonstrated an average of a) 135.4% increases in load during low-price periods, b) >36.8% reductions in electricity peak-price period, and c) 6.2% electricity cost savings.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A typology of organizational readiness for change based on a latent profile analysis

Companies have to undergo many change processes to succeed in the transforming economy. However, many change processes fail because employees are insufficiently accompanied through the process in a targeted manner. This study of N = 427 employees from a steel industry company undergoing a transformation process examines whether the organizational readiness for change (ORC) of highly affected employees can be classified into profiles, how these profiles can be predicted by various antecedents, and whether outcome variables such as job satisfaction can be explained by profile membership. Based on five facets of ORC (i.e., individual valence and positive affect), a total of six ORC profiles were identified: Proactives , Acceptors , Opens , Neutrals , Reluctants and Deniers . Employees’ optimism and the degree of perceived interpersonal and informational fairness can predict profile membership. It was shown that profiles significantly differ in relevant outcome variables satisfaction and intention to leave. These results contribute to the basic understanding of ORC and provide an initial approach for improving ORC profiles which could increase the success rate of change processes in companies.

Köhler, Alina↗

Characterizing patterns and variability of building electric load profiles in time and frequency domains

The rapid development of advanced metering infrastructure provides a new data source—building electrical load profiles with high temporal resolution. Electric load profile characterization can generate useful information to enhance building energy modeling and provide metrics to represent patterns and variability of load profiles. Such characterizations can be used to identify changes to building electricity demand due to operations or faulty equipment and controls. In this study, we proposed a two-path approach to analyze high temporal resolution building electrical load profiles: (1) time-domain analysis and (2) frequency-domain analysis. Furthermore, the commonly adopted time-domain analysis can extract and quantify the distribution of key parameters characterizing load shape such as peak-base load ratio and morning rise time, while a frequency-domain analysis can identify major periodic fluctuations and quantify load variability. We implemented and evaluated both paths using whole-year 15-minute interval smart meter data of 188 commercial office building in Northern California. The results from these two paths are consistent with each other and complementary to represent full dynamics of load profiles. The time- and frequency-domain analyses can be used to enhance building energy modeling by: (1) providing more realistic assumptions about building operation schedules, and (2) validating the simulated electric load profiles using the developed variability metrics against the real building load data.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Two-Step Time-Series Data Clustering Method for Building-Level Load Profile

Residential and commercial buildings have huge potential to contribute value to improve grid resilience by participating grid services. To reveal the significant value, it is critical to estimate the grid service capability from these buildings. Unlike the large-scale distributed energy resources such as wind and solar farms, those buildings need to participate grid services in aggregation, not by individual. Therefore, it is important to appropriately group buildings for aggregation. The load profiles in the same group will have similar characteristics at the same time step, so grid operators can send the grid service signal to the customer group with a higher chance to respond at that time step. In this paper, we develop a load profile clustering method to classify the building-level load profiles for grid service capability estimation. In our two-step clustering approach, we first calculate the total load consumption for each building, clustering the load profiles based on energy consumption level. Then, we further cluster the load profiles in each energy cluster based on the load shape. The parameter selection for each clustering step is discussed. The proposed method is applied on actual building-level load profiles, and the results have proved the effectiveness of this method.

advanced metering infrastructure (AMI)↗

Profiling the BLAST bioinformatics application for load balancing on high-performance computing clusters

Abstract Background The Basic Local Alignment Search Tool (BLAST) is a suite of commonly used algorithms for identifying matches between biological sequences. The user supplies a database file and query file of sequences for BLAST to find identical sequences between the two. The typical millions of database and query sequences make BLAST computationally challenging but also well suited for parallelization on high-performance computing clusters. The efficacy of parallelization depends on the data partitioning, where the optimal data partitioning relies on an accurate performance model. In previous studies, a BLAST job was sped up by 27 times by partitioning the database and query among thousands of processor nodes. However, the optimality of the partitioning method was not studied. Unlike BLAST performance models proposed in the literature that usually have problem size and hardware configuration as the only variables, the execution time of a BLAST job is a function of database size, query size, and hardware capability. In this work, the nucleotide BLAST application BLASTN was profiled using three methods: shell-level profiling with the Unix “time” command, code-level profiling with the built-in “profiler” module, and system-level profiling with the Unix “gprof” program. The runtimes were measured for six node types, using six different database files and 15 query files, on a heterogeneous HPC cluster with 500+ nodes. The empirical measurement data were fitted with quadratic functions to develop performance models that were used to guide the data parallelization for BLASTN jobs. Results Profiling results showed that BLASTN contains more than 34,500 different functions, but a single function, RunMTBySplitDB, takes 99.12% of the total runtime. Among its 53 child functions, five core functions were identified to make up 92.12% of the overall BLASTN runtime. Based on the performance models, static load balancing algorithms can be applied to the BLASTN input data to minimize the runtime of the longest job on an HPC cluster. Four test cases being run on homogeneous and heterogeneous clusters were tested. Experiment results showed that the runtime can be reduced by 81% on a homogeneous cluster and by 20% on a heterogeneous cluster by re-distributing the workload. Discussion Optimal data partitioning can improve BLASTN’s overall runtime 5.4-fold in comparison with dividing the database and query into the same number of fragments. The proposed methodology can be used in the other applications in the BLAST+ suite or any other application as long as source code is available.

59 BASIC BIOLOGICAL SCIENCES↗

Synergistic Retrieval of Temperature and Humidity Profiles from Space-Based and Ground-Based Infrared Sounders Using an Optimal Estimation Method

The atmospheric temperature and humidity profiles of the troposphere are generally measured by radiosondes and satellites, which are essential for analyzing and predicting weather. Nevertheless, the insufficient observation frequencies and low detection accuracy of the boundary layer restricts the description of atmospheric state changes by the temperature and humidity profiles. Therefore, this work focus on retrieving the temperature and humidity profiles using observations of the FengYun-4 (FY-4) Geostationary Interferometric Infrared Sounder (GIIRS) combined with ground-based infrared spectral observations from the Atmospheric Emitted Radiance Interferometer (AERI), which are more accurate than space-based individual retrieval results and have a wider effective retrieval range than ground-based individual retrieval results. Based on the synergistic observations, which are made by matching the space-based and ground-based data with those of different spatial and temporal resolutions, a synergistic retrieval process is proposed to obtain the temperature and humidity profiles at a high frequency under clear-sky conditions based on the optimal estimation method. In this research, using the line-by-line radiative transfer model (LBLRTM) as the forward model for observing simulations, a retrieval experiment was carried out in Qingdao, China, where an AERI is situated. Taking radiosonde data as a reference for comparing the retrieval results of the temperature and humidity profiles of the troposphere, the root-mean-square error (RMSE) of the synergistic retrieval algorithm below 400 hPa is within 2 K for temperature and within 12% for relative humidity. Compared with the GIIRS individual retrieval, the RMSE of temperature and relative humidity for the synergistic method is reduced by 0.13–1.5 K and 2.7–4.4% at 500 hPa, and 0.13–2.1 K and 2.5–7.2% at 900 hPa. Moreover, the forecast index (FI) calculated from the retrieval results shows reasonable consistency with the FIs calculated from the ERA5 reanalysis and from radiosonde data. The synergistic retrieval results have higher temporal resolution than space-based retrieval results and can reflect the changes in the atmospheric state more accurately. Overall, the results demonstrated the promising potential of the synergistic retrieval of temperature and humidity profiles at high accuracy and high temporal resolution under clear-sky conditions from FY-4/GIIRS and AERI.

54 ENVIRONMENTAL SCIENCES↗

Cross Inference of Throughput Profiles Using Micro Kernel Network Method

Dedicated network connections are being increasingly deployed in cloud, centralized and edge computing and data infrastructures, whose throughput profiles are critical indicators of the underlying data transfer performance. Due to the cost and disruptions to physical infrastructures, network emulators, such as Mininet, are often used to generate measurements needed to estimate throughput profiles, typically expressed as a function of the connection round trip time. The profiles estimated using measurements from such emulated networks are usually inaccurate for high bandwidth and high latency connections, since they do not accurately reflect the critical network transport dynamics mainly due to computing and memory constraints of the host. We present a machine learning (ML) method to estimate the throughput profiles using emulation measurements to closely match the testbed and production network profiles. In particular, we propose a micro Kernel Network (mKN) that provides baseline throughput measurements on the host running Mininet emulations, which are used to learn a regression map that converts them to the corresponding testbed measurement estimates. Once initially learned, this map is applied to measurements from subsequent network emulations on the same host. We present experimental measurements to illustrate this approach, and derive generalization equations for the proposed mKN-ML method. Using a four-site scenario emulation, we show the effectiveness of this method in providing accurate concave throughput profiles from inaccurate convex or non-smooth ones indicated by Mininet emulation.

Rao, Nageswara↗

Revisiting the definition of field capacity as a functional parameter in a layered agronomic soil profile beneath irrigated maize

The soil water content at the condition of field capacity (θ FC ) is a key parameter in irrigation scheduling and has been suggested to be determined by running a synthetic drainage experiment until the flux rate (q) at the bottom of the soil profile achieves a predefined negligible value (q FC ). We question the impact of q FC on the assessment of field capacity. Moreover, calculating θ FC as the integral mean of the water content profile when q is equal to q FC is strictly valid only for uniform soil profiles. By contrast, this practice is ambiguous and biased for stratified soil profiles due to the soil water content discontinuity at the layer interfaces. In this study, the concept of field capacity was revisited and adapted to practical agronomic heuristics. By resorting to the assessment of root-zone water storage capacity (W), we envision field capacity as a functional hydraulic parameter derived from synthetic irrigation scheduling scenarios to minimize drought stress, drainage, and nitrate leachate below the root zone. A functional analysis was carried out on a 135-cm-thick layered soil profile beneath maize in eastern Nebraska. On-farm irrigation scheduling applications and agricultural practices were recorded for 20 years (2001–2020) at a daily time step. Hydrus-1D was calibrated and validated with direct measurements of the soil water retention curve and soil water content data, respectively, in each soil layer. A set of functional field capacity values was derived from 24 irrigation scheduling scenarios, and the optimal water storage capacity at field capacity (W FC ) was approximately 50 cm (corresponding to about 80% saturation in the soil profile). An average irrigation amount of 217.5 mm distributed over 21 events was obtained by using optimal irrigation scheduling, which was initiated when the matric pressure head took on a value of –700 cm and the irrigation rate was set at 1.0 cm d –1 . This irrigation practice ensured water storage at approximately the same level (ideally at W FC ) by sustaining only evapotranspiration fluxes in the uppermost portion of the root zone and by limiting excessive drainage. This protocol can be transferred to other agricultural fields.

54 ENVIRONMENTAL SCIENCES↗

The effects of low- Z shielding on uranium isotope discrimination using the time-emission profiles of long-lived delayed neutrons

Characteristic radiation signals can often be obscured or distorted by relatively small amounts of shielding, as is the case for highly-enriched uranium (HEU). Methods that can detect and identify special nuclear material (SNM) even in the presence of shielding are therefore of high value in many nuclear security and nonproliferation applications. It has been previously shown that the buildup and decay time profiles of long-lived delayed neutron emission from induced fission are unique and can be used to successfully discriminate between uranium isotopes and infer enrichment without the need for a calibration standard. Because the long-lived delayed neutron precursors exhibit decay times on the order of tens of seconds, their time profiles are not sensitive to the delays associated with scattering or diffusion in shielding materials, which occur on much shorter time scales. The use of capture-based neutron detectors can also mitigate the effects of shielding, because unlike recoil-based detectors, they do not have a low-energy detection threshold, so that neutrons that lose their energy in the shielding are still detectable. In this work, the effects of low-Z shielding on the discrimination of uranium isotopes using delayed neutron time profile measurements are investigated for the first time. Furthermore, the delayed neutron buildup and decay time profiles of polyethylene-shielded HEU and depleted uranium objects are measured using two different capture-gated detector designs. The recorded profiles show good agreement with the shapes predicted from six-group delayed neutron models even with the shielding in place, and the enrichment of the HEU and DU objects was estimated accurately to within 5%. While this method shows potential for discriminating and characterizing extreme cases of enrichment when a moderate amount of shielding is present, simulation results suggest that significant limitations exist for intermediate enrichment levels due to distortions to the delayed neutron time profile caused by neutron reflection in the shielding.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The ellipticity parameterization for an NFW profile: An overlooked angular structure in strong lens modeling

Galaxy-scale gravitational lenses are often modeled with two-component mass profiles where one component represents the stellar mass and the second is a Navarro Frenk White (NFW) profile representing the dark matter. Outside of the spherical case, the NFW profile is costly to implement, and so it is approximated via two different methods; ellipticity can be introduced via the lensing potential (NFWp) or via the mass by approximating the NFW profile as a sum of analytical profiles (NFWm). While the NFWp method has been the default for lensing applications, it gives a different prescription of the azimuthal structure, which we show introduces ubiquitous gradients in ellipticity and boxiness in the mass distribution rather than having a constant elliptical shape. Because an unmodeled azimuthal structure has been shown to be able to bias lens model results, we explored the degree to which this azimuthal structure that was introduced can affect the model accuracy. We constructed input profiles using composite models using both the NFWp and NFWm methods and fit these mocks with a power-law elliptical mass distribution (PEMD) model with external shear. As a measure of the accuracy of the recovered lensing potential, we calculated the value of the Hubble parameter H 0 one would determine from the lensing fit. We found that the fits to the NFWp input return H 0 values that are systematically biased by about 3% lower than the NFWm counterparts. We explored whether such an effect is attributable to the mass sheet transformation (MST) by using an MST-independent quantity, ξ 2 . We show that, as expected, the NFWm mocks are degenerate with PEMD through an MST. For the NFWp, an additional bias was found beyond the MST due to the azimuthal structure exterior to the Einstein radius. We recommend modelers use an NFWm prescription in the future, such that the azimuthal structure can be introduced explicitly rather than implicitly.

79 ASTRONOMY AND ASTROPHYSICS↗