SEARCH · Engineering Papers

Results for “Constrained control”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,657 records · Page 92

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Integration of Low-Fidelity MDO and CFD-Based Redesign of Low-Boom Supersonic Transports

A mixed-fidelity low-boom multidisciplinary optimization (MDO) problem is formulated for integration of low-fidelity MDO and computational fluid dynamics (CFD) based low-boom inverse design optimization. The mixed-fidelity low-boom MDO problem aims to enforce the weight consistency for CFD-based low-boom inverse design optimization: the optimum low-boom configuration is designed for the weight at the start of cruise of the same configuration for the low-boom overland mission. Moreover, it also seeks the optimum trades among the maximum takeoff gross weight (MTOGW), cruise Mach, and range for the low-boom overland mission while meeting the requirement of flying a transatlantic overwater mission. A block coordinate optimization (BCO) method is developed to find an approximate solution of the mixed-fidelity low-boom MDO problem. The BCO method is successfully applied to generate two CFD-based low-boom configurations that closely match two reversed equivalent area targets with ground noise levels below 70 PLdB, respectively. Moreover, these CFD-based low-boom configurations can be obtained with minor wing modifications from the solutions of the low-fidelity MDO problem. The low-fidelity MDO solutions are Pareto points for constrained multiobjective optimization of MTOGW, low-boom cruise Mach, and low-boom range. The best low-boom concept has a predetermined fuselage shape tailored for passengers and main gear storage, carries 40 passengers at seat pitch of 48 in, flies a low-boom overland mission with cruise Mach of 1.8 and range of 2,950 nm, cruises overwater at Mach 1.8 with range of 3,600 nm, satisfies the specified constraints for landing/cruise/takeoff, has MTOGW of 145,164 lb, trims the low-boom cruise flight with fuel redistributions instead of control surface deflections, and has a reversed equivalent area distribution closely matching a target with ground noise level below 70 PLdB.

low-boom supersonic transports↗

NASA NTRS2021

Integration of Low-Fidelity MDO and CFD-Based Redesign of Low-Boom Supersonic Transports

low-boom supersonic transports↗

DOE OSTI2022

Biomimetic Metal-Free Hydride Donor Catalysts for CO 2 Reduction

The catalytic reduction of carbon dioxide to fuels and value-added chemicals is of significance for the development of carbon recycling technologies. One of the main challenges associated with catalytic CO 2 reduction is product selectivity: the formation of carbon monoxide, molecular hydrogen, formate, methanol, and other products occurs with similar thermodynamic driving forces, making it difficult to selectively reduce CO 2 to the target product. Significant scientific effort has been aimed at the development of catalysts that can suppress the undesired hydrogen evolution reaction and direct the reaction toward the selective formation of the desired products, which are easy to handle and store. Inspired by natural photosynthesis, where the CO 2 reduction is achieved using NADPH cofactors in the Calvin cycle, in this study we explore biomimetic metal-free hydride donors as catalysts for the selective reduction of CO 2 to formate. Here, we outline our recent findings on the thermodynamic and kinetic parameters that control the hydride transfer from metal-free hydrides to CO 2 . By experimentally measuring and theoretically calculating the thermodynamic hydricities of a range of metal-free hydride donors, we derive structural and electronic factors that affect their hydride-donating abilities. Two dominant factors that contribute to the stronger hydride donors are identified to be (i) the stabilization of the positive charge formed upon HT via aromatization or by the presence of electron-donating groups and (ii) the destabilization of hydride donors through the anomeric effect or in the presence of significant structural constrains in the hydride molecule. Hydride donors with appropriate thermodynamic hydricities were reacted with CO 2 , and the formation of the formate ion (the first reduction step in CO 2 reduction to methanol) was confirmed experimentally, providing an important proof of principle that organocatalytic CO 2 reduction is feasible. The kinetics of hydride transfer to CO 2 were found to be slow, and the sluggish kinetics were assigned in part to the large self-exchange reorganization energy associated with the organic hydrides in the DMSO solvent. Finally, we outline our approaches to the closure of the catalytic cycle via the electrochemical and photochemical regeneration of the hydride (R–H) from the conjugate hydride acceptors (R + ). We illustrate how proton-coupled electron transfer can be efficiently utilized not only to lower the electrochemical potential at which the hydride regeneration takes place but also to suppress the unwanted dimerization that neutral radical intermediates tend to undergo. Overall, this account provides a summary of important milestones achieved in organocatalytic CO 2 reduction and provides insights into the future research directions needed for the discovery of inexpensive catalysts for carbon recycling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

EGS Stimulation Design with Uncertainty Quantification at the EGS Collab Site

Engineering a robust hydraulic connection between wells is one of the most difficult aspects of enhanced geothermal systems (EGS). Designing and constructing such hydraulic connections requires and understanding of the in situ state of stress and the heterogeneities and discontinuities that naturally exist and may control the stimulation. Even with comprehensive stress and formation characterization programs substantial uncertainty remains in these key parameters. This is especially the case in high-temperature EGS environments where drilling conditions are often difficult and a far fewer logging and testing options are available. This paper presents a new approach for explicitly quantifying the uncertainties in the state of stress using a Bayesian Markov Chain Monte Carlo method. This approach produces a probability distribution for the stress tensor, including a general 3D orientation, that reflects the uncertainties in all the observations or indicators used to constrain the stress state. This method is demonstrated on the characterization data for the EGS Collab Experiment 2 site. The output of the analysis is used to guide the design of the planned stimulations. In the case of research projects like EGS Collab, explicitly quantifying the uncertainties in the stress state allow for more rigorous hypothesis testing by allowing conclusions drawn from the experiments to be interpreted in the context of the uncertain knowledge about conditions in the test bed.

Burghardt, Jeffrey A.↗

NASA NTRS1979

Development of display design and command usage guidelines for Spacelab experiment computer applications

Individual Spacelab experiments are responsible for developing their CRT display formats and interactive command scenarios for payload crew monitoring and control of experiment operations via the Spacelab Data Display System (DDS). In order to enhance crew training and flight operations, it was important to establish some standardization of the crew/experiment interface among different experiments by providing standard methods and techniques for data presentation and experiment commanding via the DDS. In order to establish optimum usage guidelines for the Spacelab DDS, the capabilities and limitations of the hardware and Experiment Computer Operating System design had to be considered. Since the operating system software and hardware design had already been established, the Display and Command Usage Guidelines were constrained to the capabilities of the existing system design. Empirical evaluations were conducted on a DDS simulator to determine optimum operator/system interface utilization of the system capabilities. Display parameters such as information location, display density, data organization, status presentation and dynamic update effects were evaluated in terms of response times and error rates.

Dodson, D. W.↗

NASA NTRS1990

Advanced launch vehicle upper stages using liquid propulsion and metallized propellants

Metallized propellants are liquid propellants with a metal additive suspended in a gelled fuel or oxidizer. Typically, aluminum particles are the metal additives. These propellants provide increase in the density and/or the specific impulse of the propulsion system. Using metallized propellants for volume- and mass-constrained upper stages can deliver modest increases in performance for Low Earth Orbit to Geosynchronous Earth Orbit and other Earth orbital transfer missions. Metallized propellants, however, can enable very fast planetary missions with a single-stage upper stage system. Trade studies comparing metallized propellant stage performance with non-metallized upper stages and the Inertial Upper Stage are presented. These upper stages are both one- and two-stage vehicles that provide the added energy to send payloads to altitudes and onto trajectories that are unattainable with only the launch vehicle. The stage designs are controlled by the volume and the mass constraints of the Space Transportation System and Space Transportation System-Cargo launch vehicles. The influences of the density and specific impulse increases enabled by metallized propellants are examined for a variety of different stage and propellant combinations.

Palaszewski, B. A.↗

NASA NTRSDate not supplied

Evaluating and Constraining Models’ Stratocumulus and Cumulus Cloud Feedbacks in the Tropics using Satellite Observations to Reduce Uncertainties in Future Climate Projections

Low-cloud feedbacks are thought to be the greatest source of uncertainty in climate projections on the centennial scale. Past research has mainly focused on the global low-cloud feedback as a whole whereas only a few studies have examined separated contributions from the two main low cloud categories: stratocumulus (Sc) and shallow cumulus (Cu). Here we aim to evaluate and constrain Sc and Cu cloud feedbacks in models using satellite observations to ultimately reduce their contribution to the spread in equilibrium climate sensitivity. We utilize the Cumulus and Stratocumulus CloudSat-CALIPSO Dataset (CASCCAD) recently developed by Cesana et al. (2019), which discriminates Sc from Cu clouds. This algorithm – based on cloud morphological characteristics at the orbital level – outperforms previous cloud-type classifications stemming from passive-sensor measurements. We investigate the observed relationship between the Sc and Cu cloud cover in tropical subsidence oceanic regions and their link to the primary low-cloud controlling factors (sea surface temperature, estimated inversion strength). To gain further insight into the Sc-Cu cloud partitioning, we explore the Sc-Cu relationship during the four meteorological seasons (DJF, MAM, JJA, SON) as well as on monthly timescale. We find a substantial negative correlation between the Sc and Cu cloud cover on annual, seasonal and monthly timescales. We finally evaluate how well climate models of the Coupled Model Intercomparison Project phase 5 (CMIP5) and phase 6 (CMIP6) agree with the observations. These results establish a framework to advance model parameterizations of the underlying physical processes driving the stratocumulus and shallow cumulus clouds (turbulence and convection).

Cloud feedbacks↗

DOE OSTI2012

Dark matter as the trigger of strong electroweak phase transition

In this paper, we propose a new possible connection between dark matter relic density and baryon asymmetry of the universe. The portal between standard model sector and dark matter not only controls the relic density and detections of dark matter, but also allows the dark matter to trigger the first order electroweak phase transition. We discuss systematically possible scalar dark matter candidates, starting from a real singlet to arbitrary high representations. We show that the simplest realization is provided by a doublet, and that strong first-order electroweak phase transition implies a lower bound on the dark matter direct detection rate. The mass of dark matter lies between 45 and 80 GeV, allowing for an appreciable invisible decay width of the Standard Model Higgs boson, which is constrained to be lighter than 130 GeV for the sake of the strong phase transition.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

Rearrangements of 2.5 Kilobases of Noncoding DNA from the Drosophila even-skipped Locus Define Predictive Rules of Genomic cis-Regulatory Logic

Rearrangements of about 2.5 kilobases of regulatory DNA located 59 of the transcription start site of the Drosophila even skipped locus generate large-scale changes in the expression of even-skipped stripes 2, 3, and 7. The most radical effects are generated by juxtaposing the minimal stripe enhancers MSE2 and MSE3 for stripes 2 and 3 with and without small ‘‘spacer’’ segments less than 360 bp in length. We placed these fusion constructs in a targeted transformation site and obtained quantitative expression data for these transformants together with their controlling transcription factors at cellular resolution. These data demonstrated that the rearrangements can alter expression levels in stripe 2 and the 2–3 interstripe by a factor of more than 10. We reasoned that this behavior would place tight constraints on possible rules of genomic cisregulatory logic. To find these constraints, we confronted our new expression data together with previously obtained data on other constructs with a computational model. The model contained representations of thermodynamic protein–DNA interactions including steric interference and cooperative binding, short-range repression, direct repression, activation, and coactivation. The model was highly constrained by the training data, which it described within the limits of experimental error. The model, so constrained, was able to correctly predict expression patterns driven by enhancers for other Drosophila genes; even-skipped enhancers not included in the training set; stripe 2, 3, and 7 enhancers from various Drosophilid and Sepsid species; and long segments of even-skipped regulatory DNA that contain multiple enhancers. The model further demonstrated that elevated expression driven by a fusion of MSE2 and MSE3 was a consequence of the recruitment of a portion of MSE3 to become a functional component of MSE2, demonstrating that cis-regulatory ‘‘elements’’ are not elementary objects.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1983

TRAGEN: Computer program to simulate an aircraft steered to follow a specified verticle profile. User's guide

The longitudinal dynamics of a medium range twin-jet or tri-jet transport aircraft are simulated. For the climbing trajectory, the thrust is constrained to maximum value, and for descent, the thrust is set at idle. For cruise, the aircraft is held in the trim condition. For climb or descent, the aircraft is steered to follow either (a) a fixed profile which is input to the program or (b) a profile computed at the beginning of that segment of the run. For climb, the aircraft is steered to maintain the given airspeed as a function of altitude. For descent, the aircraft is steered to maintain the given altitude as a function of range-to-go. In both cases, the control variable is angle-of-attack. The given output trajectory is presented and compared with the input trajectory. Step climb is treated just as climb. For cruise, the Breguet equations are used to compute the fuel burned to achieve a given range and to connect given initial and final values of altitude and Mach number.

Source record↗

DOE OSTI2014

Constraining possible variations of the fine structure constant in strong gravitational fields with the Kα iron line

In extensions of general relativity and in theories aiming at unifying gravity with the forces of the Standard Model, the value of the ''fundamental constants'' is often determined by the vacuum expectation value of new fields, which may thus change in different backgrounds. Variations of fundamental constants with respect to the values measured today in laboratories on Earth are expected to be more evident on cosmological timescales and/or in strong gravitational fields. In this paper, I show that the analysis of the Kα iron line observed in the X-ray spectrum of black holes can potentially be used to probe the fine structure constant α in gravitational potentials relative to Earth of Δφ ≈ 0.1. At present, systematic effects not fully under control prevent to get robust and stringent bounds on possible variations of the value of α with this technique, but the fact that current data can be fitted with models based on standard physics already rules out variations of the fine structure constant larger than some percent.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1986

The expanded LaGrangian system for constrained optimization problems

Smooth penalty functions can be combined with numerical continuation/bifurcation techniques to produce a class of robust and fast algorithms for constrainted optimization problems. The key to the development of these algorithms is the Expanded Lagrangian System which is derived and analyzed in this work. This parameterized system of nonlinear equations contains the penalty path as a solution, provides a smooth homotopy into the first-order necessary conditions, and yields a global optimization technique. Furthermore, the inevitable ill-conditioning present in a sequential optimization algorithm is removed for three penalty methods: the quadratic penalty function for equality constraints, and the logarithmic barrier function (an interior method) and the quadratic loss function (an interior method) for inequality constraints. Although these techniques apply to optimization in general and to linear and nonlinear programming, calculus of variations, optimal control and parameter identification in particular, the development is primarily within the context of nonlinear programming.

Poore, A. B.↗

DOE OSTI2021

STOCHASTIC OPTIMAL POWER FLOW FOR REAL-TIME MANAGEMENT OF DISTRIBUTED RENEWABLE GENERATION AND DEMAND RESPONSE (Final Report)

To meet the grand challenge of a sustainable energy future, there has been a surge of interest in renewable energy. Today, the uncertainty associated with renewable resources is handled by using operating reserves. The high penetration of renewable resources, however, introduces difficult-to-control dynamics and challenges for power system operation. Decision support tools are necessary at the bulk system operational level to recognize and efficiently utilize renewable resources and distributed demand response products in concert with traditional grid resources. It is envisaged that responsive load can potentially have very significant cost advantages over either spinning or non-spinning ramping reserve. Critical decisions are made during hour(s)-ahead and real-time power system operation regarding the commitment and dispatch of generators to ensure power delivery is both reliable and economic. These decisions are typically made by a security constrained optimal flow, which determines future generator commitments, dispatches, and ensures adequate reserves are available in the event of a contingency (unexpected outage) or if future system conditions deviate from forecasts. However, security has been always based on a pre-specified subset of contingency constraints whose enforcement does not guarantee security under all possible future possibilities while also giving little or no weight to the likelihood of each contingent event or the severity of its consequences. Existing tools, which are based exclusively on deterministic optimization models, do not yield optimal operational decisions to address these new challenges, in terms of both reliability and cost-effectiveness. This project has focused on developing a stochastic optimal power flow (SOPF) framework, which integrates renewable resource uncertainty, load uncertainty, distributed storage (DS), demand response (DR) products, in a holistic manner to address the uncertainty associated with ever-increasing renewable resources, along with the inclusion of distributed demand response products in future power systems. A proof-of-concept problem was created using the Pennsylvania-Jersey-Maryland (PJM) power system network. Synthetic wind generation was added to the system to simulate 50% wind penetration. A 1-hour test of SOPF operation indicated more than 6% operational cost savings. The project continued by adding the Midwestern Independent System Operator (MISO) as a partner, with focus shifting from SOPF to Stochastic Look-Ahead Unit Commitment (SLAC). Unlike PJM, MISO is faced with significant renewable energy resources within its footprint and is challenged with substantial uncertainty in its operations. The SLAC distinguishes itself from existing tools that operators use. At best, today’s tools solve two to three cases independently, where one or two system parameters, such as forecasted load level (e.g., a low, base, and high forecast), are varied and the resulting scenarios are analyzed independently. The stochastic-based optimization of SLAC leverages statistical information from an ensemble of potential operational scenarios and their respective likelihood. The SLAC output can be translated into valuable information to the operator such as suggested commitments, optimal scheduling and dispatch of resources, reserve requirements at both locational and zonal resolutions, ramping availability and requirements, availability of demand response including operational guidance concerning the near-term and real-time coordination between distributed energy resources, and utilization of distributed storage resources. The developed SOPF/SLAC tool, a stand-alone tool compatible with existing EMSs, will provide system operators with unprecedented visibility, flexibility and predictability to these resources and operational guidance concerning the real-time coordination between DERs and DR/DS products. The game changing and practical impact of this disruptive technology will be dramatic and will usher in a new era in the electric power industry, wherein green energy concepts are fully embraced, and electric power costs are lowered throughout the nation.

42 ENGINEERING↗

NASA NTRSDate not supplied

Integrated GNC Modifications and Performance Assessments for the Mars Ascent Vehicle, Spin-Stabilized Upper Stage Configuration

The objective of the Mars Ascent Vehicle (MAV), an element of the Mars Sample Return (MSR) campaign, is to successfully transport Martian surface samples from the Martian surface to a stable Low-Martian Orbit. Complete autonomy is required throughout ascent, and orbital insertion is constrained by tight dispersion boundaries. Given mass requirements to merely transport the MAV onto the Martian surface, reducing the total mass of the MAV reduces stress on all other campaign elements. An unguided, spin-stabilized second stage MAV allows for lower overall vehicle mass, at the cost of reduced GNC capability. This paper presents GNC modifications developed and employed to improve the performance of an unguided, spin-stabilized second stage MAV configuration. These updates include a predictive calculation scheme to solve the spin-stabilized attitude necessary for successful orbital insertion, and a spin-stabilized control law for a preliminary low-rate first stage spin-up. Trades were conducted focusing on first and second stage low-rate and high-rate spin rates (and other factors) in MSFC's MANTIS toolkit through 6-degree-of-freedom simulation and Monte Carlo analysis, and are presented here.

Dane Erickson↗

DOE OSTI2017

Detecting Surface Changes from an Underground Explosion in Granite Using Unmanned Aerial System Photogrammetry

Efficient detection and high-fidelity quantification of surface changes resulting from underground activities are important national and global security efforts. In this investigation, a team performed field-based topographic characterization by gathering high-quality photographs at very low altitudes from an unmanned aerial system (UAS)-borne camera platform. The data collection occurred shortly before and after a controlled underground chemical explosion as part of the United States Department of Energy’s Source Physics Experiments (SPE-5) series. The high-resolution overlapping photographs were used to create 3D photogrammetric models of the site, which then served to map changes in the landscape down to 1-cm-scale. Separate models were created for two areas, herein referred to as the test table grid region and the nearfield grid region. The test table grid includes the region within ~40 m from surface ground zero, with photographs collected at a flight altitude of 8.5 m above ground level (AGL). The near-field grid area covered a broader area, 90–130 m from surface ground zero, and collected at a flight altitude of 22 m AGL. The photographs, processed using Agisoft Photoscan® in conjunction with 125 surveyed ground control point targets, yielded a 6-mm pixel-size digital elevation model (DEM) for the test table grid region. This provided the ≤3 cm resolution in the topographic data to map in fine detail a suite of features related to the underground explosion: uplift, subsidence, surface fractures, and morphological change detection. The near-field grid region data collection resulted in a 2-cm pixel-size DEM, enabling mapping of a broader range of features related to the explosion, including: uplift and subsidence, rock fall, and slope sloughing. This study represents one of the first works to constrain, both temporally and spatially, explosion-related surface damage using a UAS photogrammetric platform; these data will help to advance the science of underground explosion detection.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Optical Characterization of the Keck Array and BICEP3 CMB Polarimeters from 2016 to 2019

The BICEP/Keck experiment (BK) is a series of small-aperture refracting telescopes observing degree-scale cosmic microwave background (CMB) polarization from the South Pole in search of a primordial B-mode signature. This B-mode signal arises from primordial gravitational waves interacting with the CMB and has amplitude parametrized by the tensor-to-scalar ratio r. Since 2016, BICEP3 and the Keck Array have been observing with 4800 total antenna-coupled transition-edge sensor detectors, with frequency bands spanning 95, 150, 220, and 270 GHz. Here we present the optical performance of these receivers from 2016 to 2019, including far-field beams measured in situ with an improved chopped thermal source and instrument spectral response measured with a field-deployable Fourier transform spectrometer. As a pair differencing experiment, an important systematic that must be controlled is the differential beam response between the co-located, orthogonally polarized detectors. We generate per-detector far-field beam maps and the corresponding differential beam mismatch that is used to estimate the temperature-to-polarization leakage in our CMB maps and to give feedback on detector and optics fabrication. The differential beam parameters presented here were estimated using improved low-level beam map analysis techniques, including efficient removal of non-Gaussian noise as well as improved spatial masking. These techniques help minimize systematic uncertainty in the beam analysis, with the goal of constraining the bias on r induced by temperature-to-polarization leakage to be subdominant to the statistical uncertainty. This is essential as we progress to higher detector counts in the next generation of CMB experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Controlling Thin-Film Stress and Wrinkling during Perovskite Film Formation

Significant effort has focused on controlling the deposition of perovskite films to enable uniform films, enabling efficiencies to climb dramatically. However, little attention has been paid to the evolution of thin-film stresses during deposition and the consequent effect on film morphology. While a textured surface topology has potential benefits for light scattering, a smooth surface is desirable to enable the pinhole-free deposition of contact layers. In this paper, we show that the highly textured morphology made by popular antisolvent conversion methods arises because of in-plane compressive stress experienced during the intermediate phase of film formation where the substrate constrains the film from expanding—leading to energy release in the form of wrinkling, resulting in trenches that can be hundreds of nanometers deep with periods of several micrometers. Finally, we demonstrate that the extent of wrinkling is correlated with the rate of film conversion and that ultrasmooth films are obtained by slowing the rate of film formation.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Quantifying Soil Phosphorus Dynamics: A Data Assimilation Approach

The dynamics of soil phosphorus (P) control its bioavailability. Yet it remains a challenge to quantify soil P dynamics. In this work we developed a soil P dynamics (SPD) model. We then assimilated eight data sets of 426-day changes in Hedley P fractions into the SPD model, to quantify the dynamics of six major P pools in eight soil samples that are representative of a wide type of soils. The performance of our SPD model was better for labile P, secondary mineral P, and occluded P than for nonoccluded organic P (Po) and primary mineral P. All parameters describing soil P dynamics were approximately constrained by the data sets. The average turnover rates were labile P 0.040 g g -1 day-1, nonoccluded Po 0.051 g g -1 day -1 , secondary mineral P 0.023 g g -1 day -1 , primary mineral P 0.00088 g g -1 day -1 , occluded Po 0.0066 g g -1 day -1 , and occluded inorganic P 0.0065 g g -1 day -1 , in the greenhouse environment studied. Labile P was transferred on average more to nonoccluded Po (transfer coefficient of 0.42) and secondary mineral P (0.38) than to plants (0.20). Soil pH and organic C concentration were the key soil properties regulating the competition for P between plants and soil secondary minerals. The turnover rate of labile P was positively correlated with that of nonoccluded Po and secondary mineral P. The pool size of labile P was most sensitive to its turnover rate. Overall, we suggest data assimilation can contribute significantly to an improved understanding of soil P dynamics.

54 ENVIRONMENTAL SCIENCES↗