Search NASA⌕ Search

SEARCH · Search NASA

Results for “state analysis”

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 163 records · Page 9

Neutronics and Thermo-Fluids Simulation of Generic Pebble-Bed Fluoride-Salt-Cooled High-Temperature Reactor

The fluoride-salt-cooled high-temperature reactor (FHR) is one of the advanced reactors that has been attracting considerable interest from both the research community and the nuclear industry. To help facilitate the nuclear community's familiarity with the FHR, Kairos Power has developed a generic FHR (gFHR) benchmark. In the research performed here, this benchmark was used to assess innovative modeling methods that combine stochastic and deterministic computer codes to perform the design and analysis of the gFHR. Further, the Monte Carlo code Serpent 2 was used to generate few-group cross sections that were then used in the neutron diffusion and thermal-fluids code AGREE to perform full-core neutronics and thermal-fluids steady-state and transient core analysis. The Argonne National Laboratory code SAM was then used to model the gFHR system and to simulate the load-follow operation of the gFHR.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Cavitand-Mediated Photodimerization of Chalcones: The Effect of Supramolecular Influences and Temperature on Reaction Selectivity

The photocycloaddition (PCA) of chalcones represents an important reaction pathway for accessing substituted cyclobutanes, which is a molecular framework with utility in synthetic chemistry, materials science, and medicine. In the past, our group has demonstrated the utility of the large cavity of γ-CD as a container for encapsulating two photo reactants for directing the PCA of several classes of aryl alkenes with high stereo- and regioselectivity: the cavitand-mediated photodimerization (CMP) approach. The CMP of chalcones reported in this work further demonstrates the effectiveness of this approach as high yields of dimers were observed in the photoreactions, while they were non-reactive in the solid state and yielded only the isomerization product in homogeneous media. The γ-CD CMP of chalcones yielded predominantly dimerized products in very good to high yields (>70%), composed of a mixture of three dimers in different proportions with syn HH as the major product. Computational analysis of the ground state complex structures revealed a strong correlation between the stability of the complex and predominance of the stereoisomer in the mixture. Further insights were deduced from temperature-dependence studies, which showed a shift in dimer selectivity tending towards a single stereoisomer.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impacts of Irradiation Structural Behavior on Thermal Hydraulics Safety Analysis to Support MURR LEU Conversion

The University of Missouri Research Reactor (MURR) located in Columbia, Missouri is one of six U.S. High Performance Research Reactors (USHPRR), including one critical facility, that is actively collaborating with U. S. Department of Energy (DOE) National Nuclear Security Administration (NNSA) Material Management and Minimization (M3) Office of Reactor Conversion and Uranium Supply to convert from the use of highly enriched uranium (HEU; ≥20 wt% U-235) to low-enriched uranium (LEU; <20 wt% U-235) fuel. A new type of very high-density LEU fuel based on an alloy of uranium and 10 wt% molybdenum (U-10Mo) is expected to allow the conversion to LEU of MURR, as well as four other USHPRR. MURR has been working with the Reactor Conversion Pillar at Argonne to perform fuel element design and fuel cycle performance analyses, steady-state thermal hydraulics safety analysis, and accident safety analyses in preparation for the conversion of MURR and to support a preliminary safety analysis report for conversion to LEU fuels. Subsequent analyses have also been performed, including transition cycles where all-fresh LEU fuel elements are introduced upon conversion and progressing through reactor operations the core is brought to equilibrium. Thermal hydraulics safety analyses performed as part of the above have employed an assumption on channel gap reduction due to burnup-related phenomena including fuel swelling, irradiation creep, and oxide layer buildup. Recently, a series of structural analyses have been performed on the MURR LEU fuel plates and an element due to significant differences between the plate and element designs of the MURR HEU and LEU fuels. In addition, NUREG-1537 indicates that structural phenomena are to be evaluated. Two separate types of structural analyses were performed for the MURR LEU fuel element: fluid-structure interaction (FSI) and irradiation thermo-mechanical. The FSI analysis evaluated the effects of hydraulic forces on the MURR LEU fuel element to quantify the flow-induced plate deflection, and a minimal impact to the channel gap thickness was predicted under prototypic and bounding conditions. The irradiation thermo-mechanical analysis evaluated the effects of fuel swelling, irradiation creep, and thermal expansion for the MURR LEU plates and the element for prototypic thermal and irradiation conditions based on a high-fidelity approach multiphysics approach. Overall, this thermo-mechanical analysis predicts smaller gap thickness changes in previously limiting regions. Larger changes are predicted in the middle of channels, and for end channels where power density is not typically a maximum. An additional thermo-mechanical analysis was performed for the outermost HEU fuel plate, which showed a similar magnitude of deflection as the outermost LEU plate. Due to substantial differences between the channel gap reductions assumed for the previous safety analyses and those predicted by the irradiation thermo-mechanical analysis, a need to evaluate their impact on the thermal hydraulics safety analyses arose. This report presents the results from the steady-state safety analyses for normal operation as well as the accident analyses for the two most limiting accident scenarios.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Exciton Enhanced Nonlinear Optical Responses in Monolayer h-BN and MoS 2 : Insight from First-Principles Exciton-State Coupling Formalism and Calculations

Excitons are vital in the photophysics of materials, especially in low-dimensional systems. The conceptual and quantitative understanding of excitonic effects in nonlinear optical (NLO) processes is more challenging compared to linear ones. Here, we present an ab initio approach to second-order NLO responses, incorporating excitonic effects, that employs an exciton-state coupling formalism and allows for a detailed analysis of the role of individual excitonic states. Taking monolayer h-BN and MoS 2 as two prototype 2D materials, we calculate their second harmonic generation (SHG) susceptibility and shift current conductivity tensor. We find strong excitonic enhancement in the NLO responses requires that the resonant excitons are not only optically bright themselves but also able to couple strongly to other bright excitons. Our results explain the occurrence of two strong peaks in the SHG of monolayer h-BN and why the A and B excitons of MoS 2 unexpectedly exhibit minimal excitonic enhancement in both SHG and shift current generation.

36 MATERIALS SCIENCE↗

First-Principles Cost Analysis of Advanced High-Temperature Nuclear Plants

Due to the vast number of recent nuclear reactor innovations, particularly those pertaining to generation IV reactor types like high-temperature gas cooled reactors (HTGRs) and sodium-cooled fast reactors (SFRs), and the newer deployment strategies envisioned, such as use of small modular reactors (SMRs) or even microreactors, reliable, detailed, and complete costs of these nuclear innovations are needed in wide availability. The types of models that generally achieve these objectives are those incorporating the fundamental nature of the real-world systems they aspire to predict, such as first-principles models. Furthermore, first principles models typically offer predictiveness that is not attained by most other types of individual-models. However, detailed, first-principles cost modeling of nuclear reactors and entire nuclear plants is relatively limited. To address this limitation in the availability of detailed, predictive models based on fundamentals, we recently developed a range of cost models, mostly based on first-principles methodologies, to project full lifecycle costs (LCCs) of nuclear power plants (NPPs) based on multiple parallel SM-HTG-pebble bed reactors (PBRs) and SM-SFRs.

Prosser, Jacob H. [Strategic Analysis, Inc., Arlin↗

Bridging the Gap on Data and Analysis for Distribution System Planning: Information That Utilities Can Provide Regulators, State Energy Offices and Other Stakeholders

Electric utilities conduct planning annually to ensure their distribution system meets technical standards, policies, and regulations; addresses forecasted grid conditions; satisfies customer needs; and advances utility priorities. The plan identifies grid deficiencies, analyzes potential solutions, and prioritizes capital investments and other expenditures. About 20 U.S. states and jurisdictions require regulated utilities to file some type of distribution system plan with the public utility commission for review. Requirements for sharing distribution system data and analyses vary widely, from few specific requirements to a detailed list of information that must be provided. While utilities conduct extensive analysis to develop distribution system plans, in most jurisdictions regulators and stakeholders do not know what data are available and how the utility uses the data in planning and investing. This report aims to bridge the gap by increasing understanding of the types of data and analyses utilities employ to develop distribution system plans and how the information affects their decision-making. The report describes information that states and stakeholders can ask for related to 11 data categories: -Forecasting loads and distributed energy resources (DERs) -Scenario analysis -Worst-performing circuits -Asset management strategy -Hosting capacity analysis -Value of DERs -Grid needs assessment -Cost-effectiveness framework for investments -Distribution system investment strategy and implementation -Geotargeted programs -Non-wires alternatives procurements.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Transforming US agriculture for carbon removal with enhanced weathering

Abstract Enhanced weathering (EW) with agriculture uses crushed silicate rocks to drive carbon dioxide removal (CDR) 1,2 . If widely adopted on farmlands, it could help achieve net-zero emissions by 2050 2–4 . Here we show, with a detailed US state-specific carbon cycle analysis constrained by resource provision, that EW deployed on agricultural land could sequester 0.16–0.30 GtCO 2 yr −1 by 2050, rising to 0.25–0.49 GtCO 2 yr −1 by 2070. Geochemical assessment of rivers and oceans suggests effective transport of dissolved products from EW from soils, offering CDR on intergenerational timescales. Our analysis further indicates that EW may temporarily help lower ground-level ozone and concentrations of secondary aerosols in agricultural regions. Geospatially mapped CDR costs show heterogeneity across the USA, reflecting a combination of cropland distance from basalt source regions, timing of EW deployment and evolving CDR rates. CDR costs are highest in the first two decades before declining to about US$100–150 tCO 2 −1 by 2050, including for states that contribute most to total national CDR. Although EW cannot be a substitute for emission reductions, our assessment strengthens the case for EW as an overlooked practical innovation for helping the USA meet net-zero 2050 goals 5,6 . Public awareness of EW and equity impacts of EW deployment across the USA require further exploration 7,8 and we note that mobilizing an EW industry at the necessary scale could take decades.

Science & Technology - Other Topics↗

NSTXU Diagnostic Disruption Dynamic Loading Represented by Response Spectra

This article presents the results of transient dynamic simulations of loads due to disruption eddy currents on the NSTXU vacuum vessel. Dynamic loading at diagnostic mounting locations is expressed as response spectra derived from the time history results of the dynamic structural simulations of a variety of disruption scenarios. The disruption simulations draw on a history of the project assessments of worst case disruptions for specific components. Major efforts to assess disruption loading have included the vacuum vessel which is the major structural support for the machine, as well as the passive plates (PPs), high harmonic fast wave (HHFW) antenna, and centerstack casing. Each one of these efforts included transient electromagnetic simulations producing time-dependent eddy current Lorentz loads (and in some cases halo loads) which then were applied to time-dependent structural dynamic analyses intended to obtain the proper dynamic amplification factors. In some instances, the EM model and structural model were identical allowing direct transfer of EM forces to the structural model. In other cases, the EM and structural model were not identical and the vector potential (VP) transfer method was used. The results files from these analyses were available (or re-run) to post process in ANSYS Classic time history postprocessor. In conclusion, the ANSYS command is used to create response spectra from time history data at desired points on the vessel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Systems Analysis and Optimization of Circular PET Packaging Supply Chains in the United States: Environmental and Socioeconomic Impacts

Many actions are underway at global, national, and local levels to address the plastic waste problem and transition toward a circular economy of plastics. Studies evaluating environmental and socioeconomic impacts of such a transition are lacking. Here, the purpose of this study is to conduct a national systems analysis of polyethylene terephthalate (PET) packaging supply chains in the United States. Material flow data was combined with environmental and socioeconomic indicators to evaluate and compare the sustainability of the linear PET packaging supply chain, current (2019) supply chain, and possible future circular supply chain options in the United States. Environmentally optimal circular US PET packaging material flows showed 31% and 38% savings of GHG emissions and energy demand, respectively, with a circularity of 77% when compared with a linear supply chain. Additionally, the environmentally optimal system showed higher employment (29%) and wages (31%) than a linear system, but with a 5% decrease in revenue generation. A socioeconomically optimal circular PET supply chain showed increased employment (by 52%), wages (by 67%), and revenues (by 1%), with a circularity of 59% when compared with the linear system. However, it showed 14% higher GHG emissions than a linear system, indicating a trade-off between environmentally and socioeconomically optimal circular PET packaging systems. Overall, linear-to-circular material flow transition may not necessarily lead to increased revenues and decreased environmental impacts of the entire system, but it does benefit society due to increased employment and wages. Future systems analysis work should focus on improving data quality for environmental and socioeconomic dimensions.

PET↗

Preliminary Seismic Yield Estimates of the July 1, 2025 Explosions near Esparto California

A series of large damaging explosions involving fireworks storage occurred near Esparto, Yolo County California on July 1, 2025. Three explosions were located and reported by the University of California Berkeley Seismology Laboratory (UCB/BSL) and the United States Geological Survey (USGS). Analysis of local distance (< 30 km) seismic recordings of the first blast indicates that there were actually three explosions with later blasts delayed by about 3 and 30 seconds. We measured the first arriving P-wave amplitudes on four of these events with good signal-to-noise ratios.

58 GEOSCIENCES↗

Hydrogen Production Cost from Proton-Conducting Solid Oxide Electrolysis

Rigorous stakeholder-vetted techno-economic analysis (TEA) was conducted to estimate the cost of hydrogen (H 2 ) production using Proton-Conducting Solid Oxide (PSO) electrolysis. The analysis evaluates Current (2025) and Future (2035) technology cases at centralized plant scales of 50 and 500 metric tonnes per day (MTD), assuming electricity, water, and air as the only system inputs. Untaxed, unsubsidized levelized cost of hydrogen (LCOH) is projected to range from 2020 $\$$1.81 to $\$$2.47/kg H 2 at an electricity price of $\$$0.03/kWh and 97% capacity factor under Nth-of-a-kind (NOAK) deployment assumptions. Total installed capital cost was developed using bottom-up Design for Manufacture and Assembly (DFMA) stack cost modeling and detailed balance-of-plant estimates, including mechanical and electrical subsystems, installation, site preparation, engineering, and contingency. Stack performance assumptions include thermoneutral operation, degradation over time, and periodic replacement. LCOH was calculated using the Hydrogen Analysis (H2A) discounted cash flow model in constant 2020 dollars. Results indicate PSO electrolysis has potential for competitive hydrogen production costs under low-cost electricity and mature manufacturing conditions.

08 HYDROGEN↗

Plastic additives in the ocean: Use of a comprehensive dataset for meta-analysis and method development

In excess of 13,000 chemicals are added to plastics (‘additives’) to improve performance, durability, and production of plastic products. They are categorized into numerous chemical classes including flame retardants, light stabilizers, antioxidants, and plasticizers. While research on plastic additives in the marine environment has increased over the past decade, there is a lack of methodological standardization. To direct future measurement of plastic additives, we compiled a first-of-its-kind dataset of literature assessing plastic additives in marine environments, delineated by sample type (plastic debris, seawater, sediment, biota). Using this dataset, we performed a meta-analysis to summarize the state of the science. Currently, our dataset includes 217 publications published between 1978 and May 2023. The majority of publications analyzed plastic additives in biota collected from Europe and Asia. Analyses concentrated on plasticizers, brominated flame retardants, and bisphenols. Common sample preparation techniques included Solvent - Agitation extraction for plastic, sediment, and biota samples, and Solid Phase Extraction for seawater samples with dichloromethane and solvent mixtures including dichloromethane as the organic extraction solvent. Finally, most analyses were performed utilizing gas chromatography/mass spectrometry. There are a variety of data gaps illuminated by this meta-analysis, most notably the small number of compounds that have been targeted for detection compared to the large number of additives used in plastic production. The provided dataset facilitates future investigation of trends in plastic additive concentration data in the marine environment (allowing for comparison to toxicity thresholds) and acts as a starting point for optimizing and harmonizing plastic additive analytical methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exploring Nonmotorized Travel in New York State Using 2017 National Household Travel Survey Data

This report presents systematic analysis of nonmotorized travel behavior within New York State (NYS) utilizing 2017 National Household Travel Survey data. As walking and biking assume increasingly prominent roles in advancing active transportation objectives, accessible mobility, and sustainable transportation systems, comprehensive understanding of these modal patterns becomes essential for evidence-based policy development and strategic planning initiatives. Recognizing the substantial geographic, demographic, and socioeconomic heterogeneity characterizing NYS—spanning from the concentrated urban fabric of New York City (NYC) to dispersed suburban and rural contexts—this analysis employs the structured "4Ws" analytical framework to examine participant demographics, spatial and temporal distributions, and motivational factors underlying active transportation choices relative to motorized alternatives. This methodological approach captures modal behavior patterns, user characteristics, trip purposes, and temporal variations across the state's diverse contexts. Additionally, the research examines behavioral differences among distinct user classifications, including walk-only and bike-only travelers, while conducting comparative analysis of nonmotorized travel patterns between ALICE (Asset Limited, Income Constrained, Employed) and non-ALICE household categories.

99 GENERAL AND MISCELLANEOUS↗

R ( 3780 ) Resonance Interpreted as the 1 D 1 3 -Wave Dominant State of Charmonium from Precise Measurements of the Cross Section of e + e − → Hadrons

We report the precise measurements of the cross section of e + e − → hadrons at center-of-mass energies from 3.645 to 3.871 GeV. We thereby perform the most precise study of the cross sections and find a complex system composed of three resonances of R ( 3760 ) , R ( 3780 ) , and R ( 3810 ) . For the first time, we measure the R ( 3810 ) electronic width to be ( 19.4 ± 7.4 ± 12.1 ) eV . For the R ( 3760 ) resonance, we measure the mass to be ( 3751.9 ± 3.8 ± 2.8 ) MeV / c 2 , the total width to be ( 32.8 ± 5.8 ± 8.7 ) MeV , and the electronic width to be ( 184 ± 75 ± 86 ) eV . For the R ( 3780 ) resonance, we measure its mass to be ( 3778.7 ± 0.5 ± 0.3 ) MeV / c 2 , total width to be ( 20.3 ± 0.8 ± 1.7 ) MeV , and electronic width to be ( 265 ± 67 ± 83 ) eV . Forty-seven years ago, the ψ ( 3770 ) resonance was discovered, and was subsequently interpreted as the 1 D 1 3 -wave dominant state of charmonium. However, our analysis of the total-hadron cross sections indicates that the ψ ( 3770 ) is not a single state, but a complex system composed of the R ( 3760 ) , R ( 3780 ) , and R ( 3810 ) resonances. Among these, we interpret the R ( 3780 ) is a resonance dominated by the 1 D 1 3 charmonium state. Published by the American Physical Society 2024

Ablikim, M.↗

STAT7 v2.0 Verification and Validation Report

This report documents the software testing which has been performed for STAT7 Version 2.0 (hereafter STAT7 v2.0), a code for the statistical propagation of uncertainties in steady-state thermal hydraulics (TH) analysis of plate-fueled research and test reactors. The key capabilities addressed in this testing were selected based on the use of the software for analysis work supporting the conversion to low-enriched uranium fuel of U.S. High-Performance Research Reactors (USHPRR) such as the Massachusetts Institute of Technology Research Reactor (MITR-II). This testing scope includes comparison of the code’s statistical processing capabilities using statistical tests, comparison of coolant property evaluations to NIST data tables, and systematic comparison of key TH solver capabilities to the results of hand calculations and to predictions from PLTEMP/ANL Version 4.4. The

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Structural and Spectroscopic Characterization of Plutonium and Other Tetravalent Metals Complexed to a Keggin Ion

Here, we report the isolation of the first plutonium(IV) complex with a Keggin ion chelator: Cs 20 [Pu(PW 11 O 39 ) 2 ] 2 ·13H 2 O. Single crystal XRD and solid-state UV–vis absorbance analysis demonstrate the stabilization of Pu 4+ by the Keggin ligand. The unit cell contains two [Pu(PW 11 O 39 ) 2 ] 10– complexes (Pu(PW 11 ) 2 ) bridged by Cs + . Raman and 31 P NMR spectra of Pu(PW 11 ) 2 are consistent with the analogous Zr 4+ , Hf 4+ , Ce 4+ , and Th 4+ complexes. The Pu–O bond distances at the two Pu sites are 2.35(3) and 2.34(3) Å, matching the value extrapolated from the bonding trend built with the other 8-coordinated tetravalent cations. However, the long-range arrangement of the Pu(PW 11 ) 2 complexes within the lattice is unique in the series of M IV (PW 11 ) 2 compounds: pairs of Pu(PW 11 ) 2 are organized perpendicular to each other. Based on solution-state UV–visible absorbance, small-angle X-ray scattering (SAXS), and 31 P NMR, the tetravalent cations quantitatively form the 1:2 species in solution ([Pu(PW 11 O 39 ) 2 ] 10– (aq) ) and no 1:1 species ([Pu(PW 11 O 39 )(H 2 O) x ] 3– (aq) ). Finally, a linear correlation exists between the metal–oxygen distances in the M IV (PW 11 ) 2 compounds and the corresponding metal dioxides, allowing for extrapolation for Pa 4+ , Am 4+ , and Bk 4+ . The results indicate that our microscale POM approach represents a viable pathway to probe properties of rare actinide ions in discrete molecules, beyond the traditional oxide extended solids.

and nuclear chemistry↗

Site-specific Design Case Study for Wet Waste Hydrothermal Liquefaction and Biocrude Upgrading to Hydrocarbon Fuels

Hydrothermal liquefaction (HTL) is a thermal process that converts wet biomass to renewable hydrocarbon fuel blendstocks (i.e., renewable naphtha, renewable diesel, and sustainable aviation fuel (SAF)). It can utilize a wide range of pure and blended wet feedstocks, including sewage sludge from water resource recovery facilities (WRRF), food and agriculture wastes, algae, fats, oils and greases (FOG) and blends of dry and wet wastes/feedstocks. Historically, techno-economic analysis (TEA) and annual state of technology (SOT) assessments with standard economic assumptions used by the Bioenergy Technologies Office (BETO) were conducted for the wet waste HTL pathway leveraging experimental data collected from Pacific Northwest National Laboratory’s (PNNL) continuous flow reactor systems. The objective of the SOT assessment has been to guide and track progress of BETO’s HTL research and development (R&D) toward reduced cost and greenhouse gas (GHG) emissions for the pathway. However, gaps exist between BETO’s traditional SOT updates and the needs of key external stakeholders that – if addressed – will accelerate technology adoption. This Business Case Study aims to bridge this gap by providing an updated design, TEA, and LCA based on PNNL’s FY23 R&D with added analyses and information that provide enhanced relevance for stakeholders of the HTL technology. This includes specific siting, regional wet waste resource inventory and transportation cost analyses, fuel market information, sustainable fuel policy impacts, economic metrics of net present value (NPV) and internal rate of return (IRR), greenhouse gas (GHG) emissions analysis, and statistical analysis of cost and technical uncertainties of the HTL plant design. The study focuses on the “Detroit combined statistical area (CSA)” region for siting of a wet waste HTL plant adjacent to the Great Lakes Water Authority (GLWA) facility with guidance from industry participants. Regional resource and siting analyses were conducted to identify feedstock availability, scale, and cost, as well as a beneficial site location. TEA with detailed rigorous capital cost estimation for the specific site application was conducted to evaluate the key economic metrics of most value to industrial partners. These include total capital investment, operating costs, minimum fuel selling price (MFSP) of the biocrude and fuel blendstock, and NPV and internal rate of return IRR with sustainable fuel credits. Life cycle analysis was conducted to evaluate the supply chain greenhouse gas (GHG) emissions for the wet waste HTL process as compared with petroleum derived diesel. This study is also informed by years of R&D and process de-risking learnings and was conducted with a basic engineering HTL plant design and costing that akin to a “first-of-a-kind” plant economics. This differs from our conventional “nth plant ” SOT assessments. Specifically, the HTL process model has been updated with more operationally reliable methods for feed heating and phase separations. Further, we have implemented additional spare equipment for redundancy, a more rigorous installed equipment cost estimation approach, and additional costs associated with feed formatting and delivery, building, piping and site development. An Excel-based cost sheet based on the basic engineering design is also released alongside the report that allows users to conduct customized TEA with their own feed composition and financial assumptions.

09 BIOMASS FUELS↗

Preserving nonlinear constraints in variational flow filtering data assimilation

Data assimilation aims to estimate the states of a dynamical system by optimally combining sparse and noisy observations of the physical system with uncertain forecasts produced by a computational model. The states of many dynamical systems of interest obey nonlinear physical constraints, and the corresponding dynamics is confined to a certain sub-manifold of the state space. Standard data assimilation techniques applied to such systems yield posterior states lying outside the manifold, violating the physical constraints. This work focuses on particle flow filters which use stochastic differential equations to evolve state samples from a prior distribution to samples from an observation-informed posterior distribution. The variational Fokker-Planck (VFP)—a generic particle flow filtering framework—is extended to incorporate non-linear, equality state constraints in the analysis. To this end, two algorithmic approaches that modify the VFP stochastic differential equation are discussed: (i) VFPSTAB, to inexactly preserve constraints with the addition of a stabilizing drift term, and (ii) VFPDAE, to exactly preserve constraints by treating the VFP dynamics as a stochastic differential-algebraic equation (SDAE). Additionally, an implicit-explicit time integrator is developed to evolve the VFPDAE dynamics. The strength of the proposed approach for constraint preservation in data assimilation is demonstrated on three test problems: the double pendulum, Korteweg-de-Vries, and the incompressible Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗