Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differential Experiments”

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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Dataset_for_Molecular_Motion_Below_the_Glass_Transition_A_Solid-State_NMR_Study_of_Siloxane_Polymer_Dynamics Study

This dataset contains solid-state 1H and 13C NMR relaxometry data, differential scanning calorimetry (DSC) data, and size exclusion chromatography (SEC/GPC) data supporting the study of sub-glass-transition (sub-Tg) molecular dynamics in a composition- and sequence-controlled series of diphenyl-substituted polysiloxanes (PDMS, 14Ph, 33Ph, 50Ph, 67Ph, and 100Ph; 0–100% diphenylsiloxane content by mole).All solid-state NMR data were acquired on a 200 MHz Bruker Avance III HD spectrometer using a static 7 mm HX probe or a 4 mm HX probe under 4 kHz magic-angle spinning. Raw Bruker TopSpin experiment folders are included for: (1) variable-temperature 1H lineshape measurements used to determine linewidth (FWHM) as a function of temperature across the glass transition; (2) 1H T1 (saturation recovery with solid-echo detection), probing nanosecond-scale dynamics near the 1H Larmor frequency; (3) 1H T1rho (direct spin-lock, 62.5 kHz), probing microsecond-scale segmental dynamics; (4) 13C-detected Lee–Goldburg cross-polarization 1H T1rho (LGCPH T1rho) for 33Ph and 50Ph, resolving aromatic and aliphatic proton environments; and (5) 13C T1 relaxation for 33Ph and 50Ph. Differential scanning calorimetry data (TA Instruments DSC 25, −150 to +120 °C, up to +300 °C for 100Ph, 10 °C/min) are included for all six compositions and support the glass-transition temperatures in Table 1 and Figure 1. Size exclusion chromatography data (Agilent 1200 Series, PL-Gel 300 mixed-C column, THF mobile phase, polystyrene calibration standards) are included for the three synthesized copolymers (33Ph, 50Ph, 67Ph) and support the number-average molecular weights in Table 1. Processed data include per-composition relaxation-time summaries (Excel), curve-fitting and Bloembergen-Purcell-Pound (BPP) model analysis notebooks (Jupyter/Python), and Igor Pro (.pxp) master files used to generate the manuscript's figures.

Bloembergen-Purcell-Pound theory↗
DOE OSTI2022

Silicon alloying enhances fast heating rate combustion of aluminum particles

Alloys and composites of aluminum (Al) have shown promise in regulating and enhancing particle combustion for energy generation applications. Recent work on aluminum-silicon (Al-Si) spherical alloy particles has shown improved combustion at low heating rates through enhanced diffusion accompanying a lower melting temperature. This study extends reactivity analysis to higher heating rates, comparing oxidation of Al-Si with Al. Flame speeds of Al-Si powder mixed with molybdenum trioxide (MoO 3 ) powder (Al-Si+MoO 3 ) exhibited a faster transition to steady propagation relative to Al+MoO 3 . Bomb calorimetry experiments revealed up to 5.8% greater early temperature rise for Al-Si powder. Rapid steady propagation for Al-Si mixture and faster temperature rise for Al-Si particles were attributed to accelerated kinetics evidenced in thermal equilibrium analysis of the mixtures using a differential scanning calorimeter (DSC). Larger DSC exotherms in the early stages of oxidation (i.e., 480–720 °C) correlated with early heat release and promoted steady flame propagation for Al-Si+MoO 3 compared with Al+MoO 3 . Furthermore, in multiple heating rate DSC studies, Al-Si+MoO 3 ignited at a lower heating rate (i.e., 15 ⁰C/min) than Al+MoO 3 (i.e., 20 ⁰C/min). Both mixtures ignited at temperatures less than 625 °C which is above the melting temperature of Al-Si (574 °C) and below Al (659 °C). Thus, alloy particle fuels with enhanced diffusion-controlled kinetics promote steady flame propagation and show promise for energy generation applications. This is especially promising for technologies driven by enhanced condensed phase combustion such as some propellant additives or applications dependent on reliable burn rates such as pyrotechnic time delay formulations and primers.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Light-driven hydrogen evolution via a novel pincer/no pincer mechanism including a possible concerted proton electron transfer

The road to efficient molecular catalysts for artificial photosynthesis requires a great deal of basic research to address the challenges of anthropogenic global warming. An unexpected result of such research is a novel pincer/no pincer mechanism discovered for the DFT-calculated catalytic cycle of PNP-C2 complexes. Structural parameters from theoretical modeling are supported by single crystal X-ray data for the corresponding complexes. Using the examples of [M(II)Cl(PNP-C2-R)] + (M = Pd 2+ , Ni 2+ ; R = CH 3 , H; PNP-C2 = N,N-bis{(di(2-methoxy-phenyl)phosphanyl)ethyl}-N-alkylamine) this pincer/no pincer mechanism is supported by differential pulse voltammetric, cyclic voltammetric, UV-Vis spectroscopic (combined with voltammetry) and 31 P nuclear magnetic resonance spectroscopic studies. The efficiency of hydrogen evolution during artificial photosynthesis depends significantly on the ratio of photosensitiser and water reduction catalysts. At first sight, the results from irradiation experiments (gas chromatography and mass spectrometry) seem obvious, but show a considerable impact of secondary and tertiary amines within the pincer structure. The best result obtained for the palladium pincer complex was with a secondary amine as a proton relay, delivering a turnover number of 2237 after only 2.5 hours. In contrast the analogous water reduction catalyst with a methylated (tertiary) amine instead of a secondary one gives only 972 turnovers under the same experimental conditions, over the same time. For all photocatalytic measurements the amount of hydrogen is approximately 2.5-times higher when a secondary amine is used as the proton relay. The influences of these different amines on the activity of water reduction catalysts during irradiation experiments in H 2 O or D 2 O lead to evidence of a kinetic isotope effect, supporting a possible concerted proton electron transfer. DFT calculations of reaction enthalpies, activation energies and intrinsic reaction coordinates of the catalytic cycle support the high efficiency of the presented catalysts. In particular, results from secondary amines suggest new candidates as examples of the rarely known concerted proton electron transfer.

14 SOLAR ENERGY↗
DOE OSTI2025

Analyzing the Impact of Future Weather Data on Energy Consumption in Weatherization Assistant

This study supports the mission of the U.S. Department of Energy’s Weatherization Assistance Program (WAP), which aims to increase the energy efficiency of dwellings and reduce their total residential expenditures. Specifically, we examine how projected future climate conditions may affect residential building energy performance by integrating future weather data into the National Energy Audit Tool (NEAT). Since WAP evaluates the cost-effectiveness of retrofit measures over lifespans of up to 30 years, accounting for evolving climate conditions is increasingly important. To reflect future household energy demands, this study replaces historically based Typical Meteorological Year (TMY3) weather inputs with Future Typical Meteorological Year (fTMY) datasets derived from global climate model (GCM) projections. A simulation-based framework was established to enable NEAT analysis under future weather conditions. This workflow involves converting EPW-format weather files into JSON inputs compatible with NEAT and generating degree-hour metrics needed for load calculations. The fTMY dataset used in this study was developed by Oak Ridge National Laboratory through downscaling of six GCMs under different emission scenarios and covers the period from 2020 to 2100. In contrast, the TMY3 dataset is based on historical weather data from 1961 to 1990. Simulations were conducted for benchmark single-family prototype buildings across ASHRAE climate zones 1–7, which cover all regions of the U.S. except the subarctic Zone 8 in northern Alaska, evaluating both heating and cooling loads under TMY3 and fTMY conditions. Four foundation types were tested, while heating systems were standardized, as NEAT does not differentiate thermal energy load by HVAC system type in its load calculations. Results show that fTMY weather input consistently yield lower heating loads and higher cooling loads across most locations, aligning with expected climate warming trends. Notably, colder regions such as zones 6A, 6B, and 7 experience marked reductions in heating load, while warmer and transitional zones, such as 2A (Lufkin, TX) and 3C (San Francisco, CA), have substantial increases in cooling loads. Although this study does not directly assess the performance of retrofit measures under future climate conditions, it provides a critical foundation for doing so. By quantifying shifts in baseline (i.e., pre-retrofit case) energy loads between historical and future weather files, the study highlights the importance of integrating climate-responsive data into audit tools. These findings will inform future efforts to evaluate the long-term effectiveness and cost-effectiveness of weatherization measures under changing climate conditions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1971

Damping characteristics of a liquid squeeze film.

Consideration of the damping characteristics of a liquid film located between two nearly parallel plane surfaces in relative normal motion which are always closely spaced compared to the dimensions describing the area of the two like surfaces. It was found that if the motion is slow enough (about 3 centipoise or less), the fluid inertia can be neglected and the resulting fluid flow conforms to Reynold's lubrication theory. Experiments were also carried out to determine how the damping factor of a liquid squeeze film varies with film viscosity, film thickness, amplitude, and frequency under free vibration of a single-degree-of-freedom system. The data are compared to computer solutions of the nonlinear differential equation of motion. With other variables held constant, the damping factor was found to decrease as the plate spacing is increased, decrease as the viscosity is decreased, and increase as the initial spacing is decreased.

Yabe, T.↗
NASA NTRS1984

Theoretical and experimental analysis of laser altimeters for barometric measurements over the ocean

The statistical characteristics and the waveforms of ocean-reflected laser pulses are studied. The received signal is found to be corrupted by shot noise and time-resolved speckle. The statistics of time-resolved speckle and its effects on the timing accuracy of the receiver are studied in the general context of laser altimetry. For estimating the differential propagation time, various receiver timing algorithms are proposed and their performances evaluated. The results indicate that, with the parameters of a realistic altimeter, a pressure measurement accuracy of a few millibars is feasible. The data obtained from the first airborne two-color laser altimeter experiment are processed and analyzed. The results are used to verify the pressure measurement concept.

Tsai, B. M.↗
NASA NTRS1990

A modified pump laser system to pump the titanium sapphire laser

As a result of the wide tunability of the titanium sapphire laser NASA has sited it to be used to perform differential absorption lidar (DIAL) measurements of H2O vapor in the upper and lower troposphere. The titanium sapphire laser can provide a spectrally narrow (0.3 to 1.0 pm), high energy (0.5 to 1.0 J) output at 727, 762, and 940 nm which are needed in the DIAL experiments. This laser performance can be obtained by addressing the line-narrowing issues in a master oscillator and the high energy requirement in a fundamental mode oscillator. By injection seeding, the single frequency property of the master oscillator can produce a line narrow high energy power oscillator. A breadboard model of the titanium sapphire laser that will ultimately be used in NASA lidar atmospheric sensing experiment is being designed. The task was to identify and solve any problem that would arise in the actual laser system. One such problem was encountered in the pump laser system. The pump laser that is designed to pump both the master oscillator and power oscillator is a Nd:YLF laser. Nd:YLF exhibits a number of properties which renders this material an attractive option to be used in the laser system. The Nd:YLF crystal is effectively athermal; it produces essentially no thermal lensing and thermally induced birefringence is generally insignificant in comparison to the material birefringence resulting from the uniaxial crystal structure. However, in application repeated fracturing of these laser rods was experience. Because Nd:YLF rods are not commercially available at the sizes needed for this application a modified pump laser system to replace the Nd:YLF laser rod was designed to include the more durable Nd:YAG laser rods. In this design, compensation for the thermal lensing effect that is introduced because of the Nd:YAG laser rods is included.

Petway, Larry B.↗
NASA NTRS1997

Feasibility Study For A Spaceborne Ozone/Aerosol Lidar System

Because ozone provides a shield against harmful ultraviolet radiation, determines the temperature profile in the stratosphere, plays important roles in tropospheric chemistry and climate, and is a health risk near the surface, changes in natural ozone layers at different altitudes and their global impact are being intensively researched. Global ozone coverage is currently provided by passive optical and microwave satellite sensors that cannot deliver high spatial resolution measurements and have particular limitations in the troposphere. Vertical profiling DIfferential Absorption Lidars (DIAL) have shown excellent range-resolved capabilities, but these systems have been large, inefficient, and have required continuous technical attention for long term operations. Recently, successful, autonomous DIAL measurements have been performed from a high-altitude aircraft (LASE - Lidar Atmospheric Sensing Experiment), and a space-qualified aerosol lidar system (LITE - Laser In-space Technology Experiment) has performed well on Shuttle. Based on the above successes, NASA and the Canadian Space Agency are jointly studying the feasibility of developing ORACLE (Ozone Research with Advanced Cooperative Lidar Experiments), an autonomously operated, compact DIAL instrument to be placed in orbit using a Pegasus class launch vehicle.

Campbell, Richard E.↗
NASA NTRS2012

Application of GRACE for Monitoring Groundwater in Data Scarce Regions

In the United States, groundwater storage is somewhat well monitored (spatial and temporal data gaps notwithstanding) and abundant data are freely and easily accessible. Outside of the U.S., groundwater often is not monitored systematically and where it is the data are rarely centralized and made available. Since 2002 the Gravity Recovery and Climate Experiment (GRACE) satellite mission has delivered gravity field observations which have been used to infer variations in total terrestrial water storage, including groundwater, at regional to continental scales. Challenges to using GRACE for groundwater monitoring include its relatively coarse spatial and temporal resolutions, its inability to differentiate groundwater from other types of water on and under the land surface, and typical 2-3 month data latency. Data assimilation can be used to overcome these challenges, but uncertainty in the results remains and is difficult to quantify without independent observations. Nevertheless, the results are preferable to the alternative - no data at all- and GRACE has already revealed groundwater variability and trends in regions where only anecdotal evidence existed previously.

Rodell, Matt↗
DOE OSTI2026

APOLLO: a facility-scale differentiable virtual accelerator for Fermilab

As the design complexity of modern accelerators grows, there is more interest in using advanced simulations that have fast execution time or yield additional insights like gradients. The FAST/IOTA facility has been working on implementing and experimentally validating an end-to-end digital twin that is both fast and gradient-aware, allowing for rapid prototyping of new software and experiments with minimal beam time costs. Our framework integrates physics and ML codes for linac and ring simulation through a set of generic interfaces between surrogate and physics-based sections. To reproduce device inputs and outputs, system state is exposed as a deterministic discrete event simulator. Because Fermilab is undergoing control system transition, both EPICS and ACNET frontends are supported. Recently, we have begun transitioning to a new community lattice standard, PALS, as well as developing standardized infrastructure for data ingest and normalization to prepare for model calibration during FAST proton injector commissioning. We discuss implementation details as well as challenges, and future plans to extend modelling to main complex proton accelerators like PIPII and Booster.

Kuklev, Nikita [Fermilab]↗
DOE OSTI2025

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Second LDEF Post-Retrieval Symposium interim results of experiment A0034

Thermal control coatings and contaminant collector mirrors were exposed on the leading and trailing edge modules of Long Duration Exposure Facility (LDEF) experiment A0034 to provide a basis of comparison for investigating the role of atomic oxygen in the stimulation of volatile outgassing products. The exposure of identical thermal coatings on both the leading and trailing edges of the LDEF and the additional modified exposure of identical coatings under glass windows and metallic covers in each of the flight modules provided multiple combinations of space environmental exposure to the coatings and the contaminant collector mirrors. Investigations were made to evaluate the effects of the natural space and the induced environments on the thermal coatings and the collector mirrors to differentiate the sources of observed material degradation. Two identical flight units were fabricated for the LDEF mission, each of which included twenty-five thermal control coatings mounted in isolated compartments, each with an adjacent contaminant collector mirror mounted on the wall. The covers of the flight units included apertures for each compartment, exposing the thermal coatings directly to the space environment. Six of these compartments were sealed with ultraviolet-grade transmitting quartz windows and four other compartments were sealed with aluminum covers. One module of this passive LDEF experiment, occupying one-sixth of a full tray, was mounted in Tray C9 (leading edge), while the other identical module was mounted in Tray C3 (trailing edge).

Linton, Roger C.↗
DOE OSTI2025

APOLLO: a facility-scale differentiable virtual accelerator at Fermilab FAST/IOTA

As the design complexity of modern accelerators grows, there is more interest in using advanced simulations that have fast execution time or yield additional insights like gradients. The FAST/IOTA facility has been working on implementing and experimentally validating an end-to-end digital twin that is both fast and gradient-aware, allowing for rapid prototyping of new software and experiments with minimal beam time costs. Our framework integrates physics and ML codes for linac and ring simulation through a set of generic interfaces between surrogate and physics-based sections. To reproduce device inputs and outputs, system state is exposed as a deterministic event loop in a specialized discrete event simulator architecture. Because Fermilab is undergoing control system transition, several APIs were implemented as final user interfaces - a fully asynchronous EPICS soft IOC, a gRPC-based Data Pool Manager (DPM), and legacy ACNET protocols. We discuss implementation details as well as challenges handling live data assimilation and future plans to extend modelling to main complex proton accelerators like PIPII and Booster.

Kuklev, Nikita [Fermilab]↗
NASA NTRS1973

Lunar surface cosmic ray experiment S-152, Apollo 16

This investigation was directed at determining the energy spectra and abundances of low energy heavy cosmic rays (0.03 E or = 150 MeV/nucleon). The cosmic rays were detected using plastic and glass particle track detectors. Particles emitted during the 17 April 1972 solar flare dominated the spectra for energies below about 70 MeV/nucleon. Two conclusions emerge from the low energy data: (1) The differential energy spectra for solar particles vary rapidly for energies as low as 0.05 MeV/nucleon for iron-group nuclei. (2) The abundance ratio of heavy elements changes with energy at low energies; heavy elements are enhanced relative to higher elements increasingly as the energy decreases. Galactic particle fluxes recorded within the spacecraft are in agreement with those predicted taking into account solar modulation and spacecraft shielding. The composition of the nuclei at energies above 70 MeV/nucleon imply that these particles originate outside the solar system and hence are galactic cosmic rays.

Fleischer, R. L.↗
DOE OSTI2021

Reducing Sample Loss in Measurement of Heat of Vaporization of Ethanol/Gasoline Blends by Differential Scanning Calorimetry/Thermogravimetric Analysis

Higher gasoline Heat of Vaporization (HOV) can enable higher compression-ratio direct injection spark ignition engines by providing evaporative cooling that effectively increases fuel knock resistance. Methods to directly measure this fuel property in complex gasoline samples are not well developed. The objective of the present study was to further improve a Differential Scanning Calorimetry/Thermogravimetric Analyzer (DSC/TGA) method to measure the total and partial HOV of gasoline. Ten market gasoline samples were chosen to have a wide range of properties to assess the method's capability across the entire volatility range with an emphasis on understanding how well the method captures the initial 10 percent (%) of sample evaporation and how much sample is left unevaporated at the end of the experiment. Modifications to the sample preparation/introduction method as well as to the instrument itself were made to reduce initial sample losses which included measurements of the HOV at 10 degrees C and 5 degrees C (in a cold chamber) as well as under (uncontrolled) ambient conditions. Experimental results from the DSC/TGA were compared to calculated total HOV results based on Detailed Hydrocarbon Analysis (DHA). Results from the two methods agreed very well with the difference being 5% or less in almost every case. In addition, the repeatability of the method was investigated by analyzing samples in triplicate at the three temperatures investigated. One valuable conclusion from the study was that the lower temperatures of 10 degrees C and 5 degrees C enabled more reproducible measurements for both total and partial HOV. This improved precision may be caused by the fact that temperature control in the cold chamber was more reliable than the ambient laboratory temperature control. Additionally, cold chamber experiments, due to ergonomic limitations, did not allow for the use of a lid on the sample pan. The reproducibility of the evaporation rate was found to be highly dependent on the pan/lid fit which can vary significantly such that elimination of the lid improved measurement precision while operation at sub-ambient temperature slowed the evaporation rate. Results detailing sample preparation and instrument modifications as well as a detailed comparison of total and partial HOV results are presented.

ADVANCED PROPULSION SYSTEMS,BIOMASS FUELS↗
DOE OSTI2020

Far-forward neutrinos at the Large Hadron Collider

We present a new calculation of the energy distribution of high-energy neutrinos from the decay of charm and bottom hadrons produced at the Large Hadron Collider (LHC). In the kinematical region of very forward rapidities, heavy-flavor production and decay is a source of tau neutrinos that leads to thousands of charged-current tau neutrino events in a 1 m long, 1 m radius lead neutrino detector at a distance of 480 m from the interaction region. In our computation, next-to-leading order QCD radiative corrections are accounted for in the production cross-sections. Non-perturbative intrinsic-kT effects are approximated by a simple phenomenological model introducing a Gaussian kT-smearing of the parton distribution functions, which might also mimic perturbative effects due to multiple initial-state soft-gluon emissions. The transition from partonic to hadronic states is described by phenomenological fragmentation functions. To study the effect of various input parameters, theoretical predictions for D± s production are compared with LHCb data on double-differential cross-sections in transverse momentum and rapidity. The uncertainties related to the choice of the input parameter values, ultimately affecting the predictions of the tau neutrino event distributions, are discussed. We consider a 3+1 neutrino mixing scenario to illustrate the potential for a neutrino experiment to constrain the 3+1 parameter space using tau neutrinos and antineutrinos. We find large theoretical uncertainties in the predictions of the neutrino fluxes in the far-forward region. Untangling the effects of tau neutrino oscillations into sterile neutrinos and distinguishing a 3+1 scenario from the standard scenario with three active neutrino flavours, will be challenging due to the large theoretical uncertainties from QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Number of constituent quark scaling of elliptic flows in high multiplicity p-Pb collisions at $\sqrt{s_{NN}}$=5.02TeV

In this work, we briefly summarize our recent study on the number of constituent quark (NCQ) scaling of hadron elliptic flows in high multiplicity p-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. With the inclusion of hadron production via the quark coalescence model at intermediate p T , the viscous hydrodynamics at low p T , and jet fragmentation at high p T , our Hydro – Coal – Frag model provides a nice description of the p T -spectra and differential elliptic flow v 2 (p T ) of pions, kaons and protons over the pT range from 0 to 6 GeV. Our results demonstrate that including the quark coalescence is essential for reproducing the observed approximate NCQ scaling of hadron v 2 at intermediate p T in experiments, indicating strongly the existence of partonic degrees of freedom and the formation of quark-gluon plasma in high multiplicity p-Pb collisions at the LHC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Learning Constitutive Relations From Soil Moisture Data via Physically Constrained Neural Networks

Abstract The constitutive relations of the Richardson‐Richards equation encode the macroscopic properties of soil water retention and conductivity. These soil hydraulic functions are commonly represented by models with a handful of parameters. The limited degrees of freedom of such soil hydraulic models constrain our ability to extract soil hydraulic properties from soil moisture data via inverse modeling. We present a new free‐form approach to learning the constitutive relations using physically constrained neural networks. We implemented the inverse modeling framework in a differentiable modeling framework, JAX, to ensure scalability and extensibility. For efficient gradient computations, we implemented implicit differentiation through a nonlinear solver for the Richardson‐Richards equation. We tested the framework against synthetic noisy data and demonstrated its robustness against varying magnitudes of noise and degrees of freedom of the neural networks. We applied the framework to soil moisture data from an upward infiltration experiment and demonstrated that the neural network‐based approach was better fitted to the experimental data than a parametric model and that the framework can learn the constitutive relations.

54 ENVIRONMENTAL SCIENCES↗