SEARCH · Search NASA

Results for “dynamic parameter estimation”

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,513 records · Page 84

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Beam-dump ceiling and its experimental implications: The case of a portable experiment

We generalize the nature of the so-called beam-dump “ceiling” beyond which the improvement on the sensitivity reach in the search for fast-decaying mediators dramatically slows down, and we point out its experimental implications that motivate tabletop-sized beam-dump experiments for the search. Light (bosonic) mediators are well-motivated new-physics particles, as they can appear in dark-sector portal scenarios and models to explain various laboratory-based anomalies. Due to their low mass and feebly interacting nature, beam-dump-type experiments, utilizing high-intensity particle beams, can play a crucial role in probing the parameter space of such visibly decaying mediators—in particular, the “prompt decay” region, where the mediators feature relatively large coupling and mass. We present a general and semianalytic proof that the ceiling effectively arises in the prompt-decay region of an experiment and show its insensitivity to data statistics, background estimates, and systematic uncertainties, considering a concrete example, the search for axion-like particles interacting with ordinary photons at three benchmark beam facilities: PIP-II at FNAL, and SPS and LHC-dump at CERN. We then identify optimal criteria to perform a cost-effective and short-term experiment to reach the ceiling, demonstrating that very short-baseline compact experiments enable access to the parameter space unreachable thus far.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2017

Boxy/Peanut/X-Shaped Bulges: Steep Inner Rotation Curve Leads to Barlens Face-on Morphology

We use stellar dynamical bulge/disk/halo simulations to study whether barlenses (lens-like structures embedded in the narrow bar component) are only the face-on counterparts of Boxy/Peanut/X-shapes (B/P/X) seen in edge-on bars, or if some additional physical parameter affects that morphology. A range of bulge-to-disk mass and size ratios are explored: our nominal parameters (B/D=0.08, r{sub eff}/h{sub r}=0.07, disk comprising two-thirds of total force at 2.2h{sub r}) correspond to typical Milky Way mass galaxies. In all models, a bar with pronounced B/P/X forms in a few Gyr, visible in the edge-on view. However, the pure barlens morphology forms only in models with sufficiently steep inner rotation curves, dV{sub cir}/dr≳5V{sub max}/h{sub r}, achieved when including a small classical bulge with B/D≳0.02 and r{sub eff}/h{sub r}≲0.1. For shallower slopes, the central structure still resembles a barlens, but shows a clear X signature even in low inclinations. A similar result holds for bulge-less simulations, where the central slope is modified by changing the halo concentration. The predicted sensitivity on the inner rotation curve is consistent with the slopes that are estimated from gravitational potentials calculated from the 3.6 μm images, for the observed barlens and X-shaped galaxies in the Spitzer Survey of Stellar Structure in Galaxies (S{sup 4}G). For inclinations <60° the galaxies with barlenses have on average twice steeper inner rotation curves than galaxies with X shapes: the limiting slope is ∼250 km s{sup −1} kpc{sup −1}. Among barred galaxies, those with barlenses have both the strongest bars and the largest relative excess of inner surface density, both in barlens regions (≲0.5h{sub r}) and near the center (≲0.1h{sub r}); this provides evidence for bar-driven secular evolution in galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Sub-percentage measure of distances to redshift of 0.1 by a new cosmic ruler

ABSTRACT Distance-redshift diagrams probe expansion history of the Universe. We show that the stellar mass-binding energy (massE) relation of galaxies proposed in our previous study offers a new distance ruler at cosmic scales. By using elliptical galaxies in the main galaxy sample of the Sloan Digital Sky Survey Data Release 7, we construct a distance-redshift diagram over the redshift range from 0.05 to 0.2 with the massE ruler. The best-fit dark energy density is 0.675 ± 0.079 for flat Λ-cold dark matter (ΛCDM) model, consistent with those by other probes. At the median redshift of 0.11, the median distance is estimated to have a fractional error of 0.34 per cent, much lower than those by supernova (SN) Ia and baryonic acoustic oscillation (BAO) and even exceeding their future capability at this redshift. The above low-$\mathit{ z}$ measurement is useful for probing dark energy that dominates at the late Universe. For a flat dark energy equation of state model (flat wCDM), the massE alone constrains w to an error that is only a factor of 2.2, 1.7, and 1.3 times larger than those by BAO, SN Ia, and cosmic microwave background (CMB), respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

A Digital Twin for an Inverter-Based Resource Power Plant: Real-time data streaming unlocks situation awareness

Here, this study presents the development and successful implementation of a digital twin specifically designed for a grid-connected IBR power plant. By integrating a reduced-order model of the IBR system and dynamically updating the grid impedance with real-time data, the digital twin effectively captures and replicates the behavior of the physical system. Its accuracy and reliability are validated through critical test scenarios, including a three-phase fault and a line-tripping event. The results confirm that the digital twin closely emulates its physical counterpart, demonstrating its strong potential for real-time analysis, system monitoring, and predictive decision making in modern power systems.

Digital twins↗

DOE OSTI2016

Leaf chlorophyll content as a proxy for leaf photosynthetic capacity

Improving the accuracy of estimates of forest carbon exchange is a central priority for understanding ecosystem response to increased atmospheric CO 2 levels and improving carbon cycle modelling. However, the spatially continuous parameterization of photosynthetic capacity (Vcmax) at global scales and appropriate temporal intervals within terrestrial biosphere models (TBMs) remains unresolved. This research investigates the use of biochemical parameters for modelling leaf photosynthetic capacity within a deciduous forest. Particular attention is given to the impacts of seasonality on both leaf biophysical variables and physiological processes, and their interdependent relationships. Four deciduous tree species were sampled across three growing seasons (2013-2015), approximately every 10 days for leaf chlorophyll content (Chl Leaf ) and canopy structure. Leaf nitrogen (N Area ) was also measured during 2014. Leaf photosynthesis was measured during 2014-2015 using a Li-6400 gas-exchange system, with A-Ci curves to model Vcmax. Results showed that seasonality and variations between species resulted in weak relationships between Vcmax normalized to 25°C (Vcmax 25 ) and N Area ( R 2 = 0.62, P < 0.001), whereas Chl Leaf demonstrated a much stronger correlation with Vcmax 25 ( R 2 = 0.78, P < 0.001). The relationship between Chl Leaf and N Area was also weak ( R 2 = 0.47, P < 0.001), possibly due to the dynamic partitioning of nitrogen, between and within photosynthetic and nonphotosynthetic fractions. The spatial and temporal variability of Vcmax 25 was mapped using Landsat TM/ETM satellite data across the forest site, using physical models to derive Chl Leaf . TBMs largely treat photosynthetic parameters as either fixed constants or varying according to leaf nitrogen content. This research challenges assumptions that simple N Area - Vcmax 25 relationships can reliably be used to constrain photosynthetic capacity in TBMs, even within the same plant functional type. It is suggested that Chl Leaf provides a more accurate, direct proxy for Vcmax 25 and is also more easily retrievable from satellite data. These results have important implications for carbon modelling within deciduous ecosystems.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2011

On-Orbit Solar Dynamics Observatory (SDO) Star Tracker Warm Pixel Analysis

This paper describes the process of identification and analysis of warm pixels in two autonomous star trackers on the Solar Dynamics Observatory (SDO) mission. A brief description of the mission orbit and attitude regimes is discussed and pertinent star tracker hardware specifications are given. Warm pixels are defined and the Quality Index parameter is introduced, which can be explained qualitatively as a manifestation of a possible warm pixel event. A description of the algorithm used to identify warm pixel candidates is given. Finally, analysis of dumps of on-orbit star tracker charge coupled devices (CCD) images is presented and an operational plan going forward is discussed. SDO, launched on February 11, 2010, is operated from the NASA Goddard Space Flight Center (GSFC). SDO is in a geosynchronous orbit with a 28.5 inclination. The nominal mission attitude points the spacecraft X-axis at the Sun, with the spacecraft Z-axis roughly aligned with the Solar North Pole. The spacecraft Y-axis completes the triad. In attitude, SDO moves approximately 0.04 per hour, mostly about the spacecraft Z-axis. The SDO star trackers, manufactured by Galileo Avionica, project the images of stars in their 16.4deg x 16.4deg fields-of-view onto CCD detectors consisting of 512 x 512 pixels. The trackers autonomously identify the star patterns and provide an attitude estimate. Each unit is able to track up to 9 stars. Additionally, each tracker calculates a parameter called the Quality Index, which is a measure of the quality of the attitude solution. Each pixel in the CCD measures the intensity of light and a warns pixel is defined as having a measurement consistently and significantly higher than the mean background intensity level. A warns pixel should also have lower intensity than a pixel containing a star image and will not move across the field of view as the attitude changes (as would a dim star image). It should be noted that the maximum error introduced in the star tracker attitude solution during suspected warm pixel corruptions is within the specified 36 attitude error budget requirement of [35, 70, 70] arcseconds. Thus, the star trackers provided attitude accuracy within the specification for SDO. The star tracker images are intentionally defocused so each star image is detected in more than one CCD pixel. The position of each star is calculated as an intensity-weighted average of the illuminated pixels. The exact method of finding the positions is proprietary to the tracker manufacturer. When a warm pixel happens to be in the vicinity of a star, it can corrupt the calculation of the position of that particular star, thereby corrupting the estimate of the attitude.

Felikson, Denis↗

DOE OSTI2022

Dynamic Data-Driven Multiscale Modeling for Predicting Microstructure and Mechanical Properties of 316L Stainless Steel

We have developed a long short-term memory stacked ensemble (LSTM-SE) surrogate modeling approach that can provide rapid predictions of microstructural evolution and the resultant mechanical properties of American Iron and Steel Institute (AISI) 316L series stainless steel (316LSS) fuel cladding under conditions of varying temperature and radiation dose rate. To acquire training data, we developed and implemented a kinetic Monte Carlo (KMC) model to simulate precipitation kinetics of M 23 C 6 , γ', and G phases within SS316L cladding. Experimentally reported precipitation kinetics of SS316L in literature were linked to the kinetic parameters of the simulated precipitation in our KMC model. The model was then used to simulate microstructure evolution under synthetically generated treatments of varying temperature and radiation dose rate, for periods of up to 3000 hours. Changes in volume fraction, number density, and particle size of precipitates were recorded, and particle area fractions were correlated using statistical methods to develop the surrogate model. Simultaneously, the mechanical properties of the simulated microstructures were evaluated using microstructure-based finite element method (FEM) analysis to determine the elastic modulus, yield stress, ultimate tensile strength, and elongation to failure of the aged microstructures. Using this approach, our surrogate model can predict precipitation behavior within 0.25% volume fraction and mechanical properties within 6% relative error from the values predicted by the KMC and FEM models using 50 training simulations as input. The trained recurrent neural network-based model can return estimations of precipitation kinetics and mechanical properties ~1000 times faster than the physics-based codes. This work demonstrates, as a proof of concept, that reactor material service lifetimes under variable service conditions can be predicted for a statistics-based model from a practicably obtainable dataset.

36 MATERIALS SCIENCE↗

DOE OSTI2019

initMIP-Antarctica: an ice sheet model initialization experiment of ISMIP6

Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1983

Crustal dynamics project session 4 validation and intercomparison experiments 1979-1980 report

As part of the Crustal Dynamics Project, an experiment was performed to verify the ability of Satellite Laser Ranging (SLR), Very Long Baseline interferometry (VLBI) and Doppler Satellite Positioning System (Doppler) techniques to estimate the baseline distances between several locations. The Goddard Space Flight Center (GSFC) lasers were in operation at all five sites available to them. The ten baselines involved were analyzed using monthly orbits and various methods of selecting data. The standard deviation of the monthly SLR baseline lengths was at the 7 cm level. The GSFC VLBI (Mark III) data was obtained during three separate experiments. November 1979 at Haystack and Owens Valley, and April and July 1980 at Haystack, Owens Valley, and Fort Davis. Repeatability of the VLBI in determining baseline lengths was calculated to be at the 2 cm level. Jet Propulsion Laboratory (JPL) VLBI (Mark II) data was acquired on the Owens Valley to Goldstone baseline on ten occasions between August 1979 and November 1980. The repeatability of these baseline length determinations was calculated to be at the 5 cm level. National Geodetic Survey (NGS) Doppler data was acquired at all five sites in January 1980. Repeatability of the Doppler determined baseline lengths results were calculated at approximately 30 cm. An intercomparison between baseline distances and associated parameters was made utilizing SLR, VLBI, and Doppler results on all available baselines. The VLBI and SLR length determinations were compared on four baselines with a resultant mean difference of -1 cm and a maximum difference of 12 cm. The SLR and Doppler length determinations were compared on ten baselines with a resultant mean difference of about 30 cm and a maximum difference of about 60 cm. The VLBI and Doppler lengths from seven baselines showed a resultant mean difference of about 30 cm and maximum difference of about 1 meter. The intercomparison of baseline orientation parameters were consistent with past analysis.

Liebrecht, P.↗

NASA NTRS2003

Merging Ocean Color Data From Multiple Missions

Oceanic phytoplankton may play an important role in the cycling of carbon on the Earth, through the uptake of carbon dioxide in the process of photosynthesis. Although they are ubiquitous in the global oceans, their abundances and dynamics are difficult to estimate, primarily due to the vast spatial extent of the oceans and the short time scales over which their abundances can change. Consequently, the effects of oceanic phytoplankton on biogeochemical cycling, climate change, and fisheries are not well known. In response to the potential importance of phytoplankton in the global carbon cycle and the lack of comprehensive data, NASA and the international community have established high priority satellite missions designed to acquire and produce high quality ocean color data (Table 6.1). Ten of the missions are routine global observational missions: the Ocean Color and Temperature Sensor (OCTS), the Polarization and Directionality of the Earth's Reflectances sensor (POLDER), Sea-viewing Wide Field-of-view Sensor (SeaWiFS), Moderate Resolution Imaging Spectrometer-AM (MODIS-AM), Medium Resolution Imaging Spectrometer (MERIS), Global Imager (GLI), MODIS-PM, Super-GLI (S-GLI), and the Visible/Infrared Imager and Radiometer Suite (VIIRS) on the NPOESS Preparatory Project (NPP) and the National Polar-orbiting Operational Environmental Satellite System (NPOESS). In addition, there are several other missions capable of providing ocean color data on smaller scales. Most of these missions contain the spectral band complement considered necessary to derive oceanic chlorophyll concentrations and other related parameters. Many contain additional bands that can provide important ancillary information about the optical and biological state of the oceans.

Gregg, Watson W.↗

NASA NTRS2019

initMIP-Antarctica: an Ice Sheet Model Initialization Experiment of ISMIP6

ice sheet model↗

NASA NTRS2020

A Coupled Analysis of Atmospheric Mass Loss and Tidal Evolution in XUV Irradiated Exoplanets: The TRAPPIST-1 Case Study

Exoplanets residing close to their stars can experience evolution of both their physical structures and their orbits due to the influence of their host stars. In this work, we present a coupled analysis of dynamical tidal dissipation and atmospheric mass loss for exoplanets in X-ray and ultraviolet (XUV) irradiated environments. As our primary application, we use this model to study the TRAPPIST-1 system and place constraints on the interior structure and orbital evolution of the planets. We start by reporting on an ultraviolet continuum flux measurement (centered around ∼1900 Å) for the star TRAPPIST-1, based on 300 ks of Neil Gehrels Swift Observatory data, and which enables an estimate of the XUV-driven thermal escape arising from XUV photodissociation for each planet. We find that the X-ray flaring luminosity, measured from our X-ray detections, of TRAPPIST-1 is 5.6 × 10(exp −4) L(*), while the full flux including non-flaring periods is 6.1 × 10(exp −5) L*, when L(*) is TRAPPIST-1ʼs bolometric luminosity. We then construct a model that includes both atmospheric mass loss and tidal evolution and requires the planets to attain their present-day orbital elements during this coupled evolution. We use this model to constrain the ratio Q' = 3Q/2k2 for each planet. Finally, we use additional numerical models implemented with the Virtual Planet Simulator VPLanet to study ocean retention for these planets using our derived system parameters.

Juliette Becker↗

DOE OSTI2022

Matrix Approach to Land Carbon Cycle Modeling

Land ecosystems contribute to climate change mitigation by taking up approximately 30% of anthropogenically emitted carbon. However, estimates of the amount and distribution of carbon uptake across the world's ecosystems or biomes display great uncertainty. The latter hinders a full understanding of the mechanisms and drivers of land carbon uptake, and predictions of the future fate of the land carbon sink. The latter is needed as evidence to inform climate mitigation strategies such as afforestation schemes. To advance land carbon cycle modeling, we have developed a matrix approach. Land carbon cycle models use carbon balance equations to represent carbon exchanges among pools. Our approach organizes this set of equations into a single matrix equation without altering any processes of the original model. The matrix equation enables the development of a theoretical framework for understanding the general, transient behavior of the land carbon cycle. While carbon input and residence time are used to quantify carbon storage capacity at steady state, a third quantity, carbon storage potential, integrates fluxes with time to define dynamic disequilibrium of the carbon cycle under global change. The matrix approach can help address critical contemporary issues in modeling, including pinpointing sources of model uncertainty and accelerating spin-up of land carbon cycle models by tens of times. The accelerated spin-up liberates models from the computational burden that hinders comprehensive parameter sensitivity analysis and assimilation of observational data to improve model accuracy. Such computational efficiency offered by the matrix approach enables substantial improvement of model predictions using ever-increasing data availability. Overall, the matrix approach offers a step change forward for understanding and modeling the land carbon cycle.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1999

Observations of Land Surface Variability Using Passive Microwave Sensing

Understanding the global variability of land surface wetness (soil moisture), skin temperature, and related surface fluxes of heat and moisture is key to assessing the importance of the land surface in influencing climate. The feasibility of producing model estimates of these quantities is being studied as part of the International Satellite Land Surface Climatology Project (ISLSCP) Global Soil Wetness Project (GSWP). In the GSWP approach, meteorological observations and analyses are used to drive global circulation models. Satellite measurements can provide independent estimates of key land surface parameters that are needed for initializing and validating the climate models and for monitoring long-term change. Satellite observations of the land surface can also be assimilated into soil models to estimate moisture in the root zone. In our research, passive microwave satellite data recorded during 1978-1987 from the Nimbus-7 Scanning Multichannel Microwave Radiometer (SMMR) are being used to examine spatial and temporal trends in surface soil moisture, vegetation, and temperature. These data include observations at C and X bands (6.6 and 10.7 GHz), which are not available on the current Special Sensor Microwave/Imager (SSM/I) and are precursors to data that will become available from the Advanced Microwave Scanning Radiometer (AMSR) on Advanced Earth Observing Satellite (ADEOS-II) and Earth Observing System (EOS) PM1 in the year 2000. A chart shows a time-series of SMMR-derived surface temperature, T-e and surface soil moisture M, retrieved on a 0.5 deg x 0.5 deg grid and further averaged over a 4 deg x 10 deg study region in the African Sahel. Also shown are National Center for Environmental Prediction (NCEP) model outputs of surface temperature, T-sfc, and soil wetness, Soil-w. The variables have been scaled to have similar dynamic ranges on the plots. The NCEP data from the NCEP Reanalysis Project are monthly averages on a 2.5 deg x 2.5 deg grid averaged over the 4 deg x 10 deg study area. Comparisons of SMMR retrievals with forecast model output show the potential of the satellite data for validating model output and monitoring long-term trends. Continuing work will extend these results to other regions to validate the retrievals more quantitatively. In preparation for the launch of AMSR, field experiments are planned in collaboration with the Global Energy and Water Cycle Experiment (GEWEX) Coordinated Enhanced Observing Period (CEOP) experiments to evaluate the satellite-derived soil moisture measurements and to demonstrate their usefulness for land surface hydrology and climate. Additional information is contained in the original.

Njoku, Eni G.↗

DOE OSTI2024

Dynamic data-driven multiscale modeling for predicting the degradation of a 316L stainless steel nuclear cladding material

Here, we have developed a long short-term memory stacked ensemble (LSTM-SE) surrogate modeling approach that can provide rapid predictions of microstructural evolution and the resultant mechanical properties of American Iron and Steel Institute (AISI) 316L series stainless steel (316LSS) fuel cladding under conditions of varying temperature and radiation dose rate. To acquire training data, we developed and implemented a kinetic Monte Carlo (KMC) model to simulate precipitation kinetics of M 23 C 6 , γ', and G phases within SS316L cladding. Experimentally reported precipitation kinetics of SS316L in literature were linked to the kinetic parameters of the simulated precipitation in our KMC model. The model was then used to simulate microstructure evolution under synthetically generated treatments of varying temperature and radiation dose rate, for periods of up to 3000 hours. Changes in volume fraction, number density, and particle size of precipitates were recorded, and particle area fractions were correlated using statistical methods to develop the surrogate model. Simultaneously, the mechanical properties of the simulated microstructures were evaluated using microstructure-based finite element method (FEM) analysis to determine the elastic modulus, yield stress, ultimate tensile strength, and elongation to failure of the aged microstructures. Using this approach, our surrogate model can predict precipitation behavior within 0.25% volume fraction and mechanical properties within 6% relative error from the values predicted by the KMC and FEM models using 50 training simulations as input. The trained recurrent neural network-based model can return estimations of precipitation kinetics and mechanical properties ~1000 times faster than the physics-based codes. This work demonstrates, as a proof of concept, that reactor material service lifetimes under variable service conditions can be predicted for a statistics-based model from a practicably obtainable dataset.

36 MATERIALS SCIENCE↗

NASA NTRS1994

TOPEX/POSEIDON tides estimated using a global inverse model

Altimetric data from the TOPEX/POSEIDON mission will be used for studies of global ocean circulation and marine geophysics. However, it is first necessary to remove the ocean tides, which are aliased in the raw data. The tides are constrained by the two distinct types of information: the hydrodynamic equations which the tidal fields of elevations and velocities must satisfy, and direct observational data from tide gauges and satellite altimetry. Here we develop and apply a generalized inverse method, which allows us to combine rationally all of this information into global tidal fields best fitting both the data and the dynamics, in a least squares sense. The resulting inverse solution is a sum of the direct solution to the astronomically forced Laplace tidal equations and a linear combination of the representers for the data functionals. The representer functions (one for each datum) are determined by the dynamical equations, and by our prior estimates of the statistics or errors in these equations. Our major task is a direct numerical calculation of these representers. This task is computationally intensive, but well suited to massively parallel processing. By calculating the representers we reduce the full (infinite dimensional) problem to a relatively low-dimensional problem at the outset, allowing full control over the conditioning and hence the stability of the inverse solution. With the representers calculated we can easily update our model as additional TOPEX/POSEIDON data become available. As an initial illustration we invert harmonic constants from a set of 80 open-ocean tide gauges. We then present a practical scheme for direct inversion of TOPEX/POSEIDON crossover data. We apply this method to 38 cycles of geophysical data records (GDR) data, computing preliminary global estimates of the four principal tidal constituents, M(sub 2), S(sub 2), K(sub 1) and O(sub 1). The inverse solution yields tidal fields which are simultaneously smoother, and in better agreement with altimetric and ground truth data, than previously proposed tidal models. Relative to the 'default' tidal corrections provided with the TOPEX/POSEIDON GDR, the inverse solution reduces crossover difference variances significantly (approximately 20-30%), even though only a small number of free parameters (approximately equal to 1000) are actually fit to the crossover data.

Egbert, Gary D.↗

DOE OSTI2022

Spatial resolution of different discretizations over long-time for the Dirac equation with small potentials

In this report we compare the long-time error bounds and spatial resolution of finite difference methods with different spatial discretizations for the Dirac equation with small electromagnetic potentials characterized by $\varepsilon \in (0, 1]$ a dimensionless parameter. We begin with the simple and widely used finite difference time domain (FDTD) methods, and establish rigorous error bounds of them, which are valid up to the time at $O(1/\varepsilon)$. In the error estimates, we pay particular attention to how the errors depend explicitly on the mesh size $h$ and time step $\tau$ as well as the small parameter $\varepsilon$. Based on the results, in order to obtain "correct" numerical solutions up to the time at $O(1/\varepsilon)$, the $\varepsilon$-scalability (or meshing strategy requirement) of the FDTD methods should be taken as $h = O(\varepsilon^{1/2})$ and $\tau = O(\varepsilon^{1/2})$. To improve the spatial resolution capacity, we apply the Fourier spectral method to discretize the Dirac equation in space. Error bounds of the resulting finite difference Fourier pseudospectral (FDFP) methods show that they exhibit uniform spatial errors in the long-time regime, which are optimal in space as suggested by the Shannon's sampling theorem. Extensive numerical results are reported to confirm the error bounds and demonstrate that they are sharp.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Simulation of decelerating streamers in inhomogeneous atmosphere with implications for runaway electron generation

The dynamics of positive and negative streamers is numerically simulated in atmospheric pressure air in the range of parameters corresponding to the streamer deceleration and termination in the middle of the discharge gap. A detailed comparison with experiments in air at constant and variable density demonstrates good agreement between the 2D simulation results and the observations. It is shown that positive and negative streamers behave in radically different ways when decelerating and stopping. When the head potential drops, the negative streamer transits to the mode in which the propagation is due to the forward electron drift. In this case, the radius of the ionization wave front increases, whereas the electric field at the streamer head decreases further and the streamer stops. Its head diameter continues to increase due to the slow drift of free electrons in the residual under-breakdown field. On the contrary, the only advancement mechanism for a positive streamer with a decreasing head potential is a decrease in the effective radius of the ionization wave, leading to a local increase in the electric field. This mechanism makes it possible to compensate for the reduction in the efficiency of gas photoionization at small head diameters. A qualitative 1D model is suggested to describe streamer deceleration and stopping for different discharge polarities. Estimates show that, during positive streamer stopping, the local electric field at the streamer head can exceed the threshold corresponding to the transition of electrons to the runaway mode when the head potential (relative to the surrounding space) decreases to~1.2 kV in atmospheric pressure air. In this case, pulsed generation of a beam of runaway electrons directed into the channel of a stopping positive streamer can occur. The energy of the formed pulsed electron beam depends on the intensity of photoionization in front of the streamer head. This energy can vary from 700 V (when increasing the photoionization rate by a factor of 10 with respect to the value in atmospheric pressure air) to 2.6 kV (when decreasing the photoionization rate by a factor of 1000). It is possible that this behavior of decelerating positive streamers can explain the observed bursts of x-ray radiation during the streamer propagation in long air gaps.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗