Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model serving”

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

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

At least 1,279 records · Page 71

Probing the Distribution of Ozone on Mars

We present the application of infrared heterodyne line shapes of ozone on Mars to those produced by radiative transfer modeling of ozone profiles predicted by photochemistry-coupled general circulation models (GCM), and to contemporaneous column abundances measured by Mars Express SPICAM. Ozone is an important tracer of photochemistry in Mars' atmosphere, serving as an observable with which to test predictions of photochemical models. Infrared heterodyne measurements of ozone absorption features on Mars have been obtained at various Martian seasons from 1988 until present at the NASA Infrared Telescope Facility (IRTF) on Mauna Kea, Hawai'i [I]. The NASAiGoddard Space Flight Center spectrometers used were the Infrared Heterodyne Spectrometer (IRHS) [2, 3] and, since 2003, the Heterodyne Instrument for Planetary Wind and Composition (HIPWAC) [4]. A description the infrared heterodyne technique applied to ground-base observations of Martian ozone can be found in [I]. The most recent measurements on February 21-24 2008 UT at Ls=35deg were made by HIPWAC on or near the Mars Express orbital path with the goal of acquiring spectra that can be directly compared to nadir observations by SPICAM.

Fast, K. E.↗

Validation for Global Solar Wind Prediction Using Ulysses Comparison: Multiple Coronal and Heliospheric Models Installed at the Community Coordinated Modeling Center

The prediction of the background global solar wind is a necessary part of space weather forecasting. Several coronal and heliospheric models have been installed and/or recently upgraded at the Community Coordinated Modeling Center (CCMC), including the Wang-Sheely-Arge (WSA)-Enlil model, MHD-Around-a-Sphere (MAS)-Enlil model, Space Weather Modeling Framework (SWMF), and Heliospheric tomography using interplanetary scintillation data. Ulysses recorded the last fast latitudinal scan from southern to northern poles in 2007. By comparing the modeling results with Ulysses observations over seven Carrington rotations, we have extended our third-party validation from the previous near-Earth solar wind to middle to high latitudes, in the same late declining phase of solar cycle 23. Besides visual comparison, wehave quantitatively assessed the models capabilities in reproducing the time series, statistics, and latitudinal variations of solar wind parameters for a specific range of model parameter settings, inputs, and grid configurations available at CCMC. The WSA-Enlil model results vary with three different magnetogram inputs.The MAS-Enlil model captures the solar wind parameters well, despite its underestimation of the speed at middle to high latitudes. The new version of SWMF misses many solar wind variations probably because it uses lower grid resolution than other models. The interplanetary scintillation-tomography cannot capture the latitudinal variations of solar wind well yet. Because the model performance varies with parameter settings which are optimized for different epochs or flow states, the performance metric study provided here can serve as a template that researchers can use to validate the models for the time periods and conditions of interest to them.

Jian, L. K.↗

Stellar Activity Effects on Moist Habitable Terrestrial Atmospheres Around M Dwarfs

Transit spectroscopy of terrestrial planets around nearby M dwarfs is a primary goal of space missions in coming decades. 3-D climate modeling has shown that slow-synchronous rotating terrestrial planets may develop thick clouds at the substellar point, increasing the albedo. For M dwarfs with T(eff) > 3000 K, such planets at the inner habitable zone (IHZ) have been shown to retain moist greenhouse conditions, with enhanced stratospheric water vapor (fH2O > 10(exp -3)) and low Earth-like surface temperatures. However, M dwarfs also possess strong UV activity, which may effectively photolyze stratospheric H2O. Prior modeling efforts have not included the impact of high stellar UV activity on the H2O. Here, we employ a 1-D photochemical model with varied stellar UV, to assess whether H2O destruction driven by high stellar UV would affect its detectability in transmission spectroscopy. Temperature and water vapor profiles are taken from published 3-D climate model simulations for an IHZ Earth-sized planet around a 3300 K M dwarf with an N2-H2O atmosphere; they serve as self-consistent input profiles for the 1-D model. We explore additional chemical complexity within the 1-D model by introducing other species into the atmosphere. We find that as long as the atmosphere is well-mixed up to 1 mbar, UV activity appears to not impact detectability of H2O in the transmission spectrum. The strongest H2O features occur in the JWST MIRI instrument wavelength range and are comparable to the estimated systematic noise floor of ~50 ppm.

Mahmuda Afrin Badhan↗

On the road to per cent accuracy – III. Non-linear reaction of the matter power spectrum to massive neutrinos

ABSTRACT We analytically model the non-linear effects induced by massive neutrinos on the total matter power spectrum using the halo model reaction framework of Cataneo et al. In this approach, the halo model is used to determine the relative change to the matter power spectrum caused by new physics beyond the concordance cosmology. Using standard fitting functions for the halo abundance and the halo mass–concentration relation, the total matter power spectrum in the presence of massive neutrinos is predicted to per cent-level accuracy, out to $k=10 \,{ h}\,{\rm Mpc}^{-1}$. We find that refining the prescriptions for the halo properties using N-body simulations improves the recovered accuracy to better than 1 per cent. This paper serves as another demonstration for how the halo model reaction framework, in combination with a single suite of standard Λ cold dark matter (ΛCDM) simulations, can recover per cent-level accurate predictions for beyond ΛCDM matter power spectra, well into the non-linear regime.

79 ASTRONOMY AND ASTROPHYSICS↗

Improving Surface PM2.5 Forecasts in the United States Using an Ensemble of Chemical Transport Model Outputs: 2. Bias Correction with Satellite Data for Rural Areas

This work serves as the second of a two-part study to improve surface PM2.5 forecasts in the continental U.S. through the integrated use of multi satellite aerosol optical depth (AOD) products (MODIS Terra/Aqua and VIIRS DT/DB), multichemical transport model (CTM) (GEOS-Chem, WRF-Chem, and CMAQ) outputs, and ground observations. In Part I of the study, an ensemble Kalman filter (KF) technique using three CTM outputs and ground observations was developed to correct forecast bias and generate a single best forecast of PM2.5 for next day over non rural areas that have surface PM2.5 measurements in the proximity of 125 km. Here, with AOD data, we extended the bias correction into rural areas where the closest air quality monitoring station is at least 125–300 km away. First, we ensembled all of satellite AOD products to yield the single best AOD. Second, we corrected daily PM2.5 in rural areas from multiple models through the AOD spatial pattern between these areas and non rural areas, referred to as “extended ground truth” or EGT, for the present day. Lastly, we applied the KF technique to reduce the forecast bias for next day using the EGT. Our results find that the ensemble of bias-corrected daily PM2.5 from three CTMs for both today and next day show the best performance. Together, the two-part study develops a multimodel and multi-AOD bias-correction technique that has the potential to improve PM2.5 forecasts in both rural and non rural areas in near real time, and be readily implemented at state levels.

Huanxin Zhang↗

OpenFOAM Simulations of Atmospheric-Entry Capsules in the Subsonic Regime

The open-source Computational Fluid Dynamics software OpenFOAM is gaining wider acceptance in industry and academia for incompressible flow simulations. To date, there has been relatively little utilization of OpenFOAM for compressible external aerodynamic applications. The numerous turbulence models available in OpenFOAM makes it an attractive option for evaluating alternate Reynolds-Averaged Navier-Stokes (RANS) turbulent models to assess separated flow on atmospheric entry vehicles in the subsonic regime, where traditional turbulent models show reduced accuracy. This paper presents simulations of an axisymmetric capsule geometry at subsonic conditions using an OpenFOAM compressible flow solver. These results are compared with results from the NASA CFD code OVERFLOW and experimental data. These OpenFOAM simulations serve as a basis to explore OpenFOAMs extended turbulence models on compressible separated flows such as found on entry capsules.

Nikaido, Ben E.↗

NASA's Virtual Product Laboratory: An Overview

The Virtual Product Laboratory (VPL) developed by NASA's Commercial Remote Sensing Program at the John C. Stennis Space Center is a tool that enables simulation, design, and verification of remote sensing systems within a software (virtual) environment. The VPL can serve industry, government, and university communities by providing a means to conduct system trade studies, optimization, visual modeling, and data product simulations entirely in a virtual environment. The VPL can serve as a complete end-to-end simulation tool capable of producing system and subsystem performance characterizations, system optimization, and simulated data products or as a means of evaluating any one component of a remote sensing system. In this paper, we present an overview of the VPL capabilities. The VPL functional areas include Requirements, Design & Analysis, Simulation, Project Management, Knowledge Base, and Help. A description of each function along with the tools and techniques used to accomplish these functions will be provided. When possible, sample VPL displays and products will be used in the presentation.

Gasser, Jerry↗

Flow matching meets biology and life science: a survey

Over the past decade, advances in generative modeling, such as generative adversarial networks, masked autoencoders, and diffusion models, have significantly transformed biological research and discovery, enabling breakthroughs in molecule design, protein generation, catalysis discovery, drug discovery, and beyond. At the same time, biological applications have served as valuable testbeds for evaluating the capabilities of generative models. Recently, flow matching has emerged as a powerful and efficient alternative to diffusion-based generative modeling, with growing interest in its application to problems in biology and life sciences. This paper presents the first comprehensive survey of recent developments in flow matching and its applications in biological domains. We begin by systematically reviewing the foundations and variants of flow matching, and then categorize its applications into three major areas: biological sequence modeling, molecule generation and design, and peptide and protein generation. For each, we provide an in-depth review of recent progress. We also summarize commonly used datasets and software tools, and conclude with a discussion of potential future directions.

59 BASIC BIOLOGICAL SCIENCES↗

Active Measurement of Plasma Parameters in Over-Dense MFE Reactor Plasma Using High Frequency Beat-Wave Resonantly Generated Thermal Bernstein Waves

UC Davis was funded to explore the combination of beat-wave excitation and laser scattering to make completely non-invasive determination of local plasma conditions. This involves the use of two high frequency mm-wave sources of differing frequencies (the pump waves) to excite plasma waves at the difference of the frequencies, which is known as the beat frequency. The difference frequency or the frequency of the beat wave is chosen so that we produce plasma waves whose fluctuating electron density is a sensitive function of plasma parameters at the location where the source waves overlap. This fluctuation can be detected from outside the plasma, by measuring the deflection of a probe beam directed at the overlap region between the two pump waves (i.e., collective scattering). This work is an extension of mm-wave frequencies of a 10 GHz beat-wave experiment conducted on the DDT toroidal device (see Fig. 1) in a magnetized plasma where the excited Langmuir wave is replaced by the magnetized upper hybrid wave. Here, two high power (> 1 kW) ~94 GHz pump sources will be employed to match the upper hybrid resonance. The primary envisioned experiment is to drive the electron Bernstein wave (EBW) in the interior of a magnetized plasma. The EBW is a hot-plasma electrostatic wave, which cannot propagate in vacuum. Roots of the EBW dispersion relation occur at frequencies slightly higher than harmonics of the electron cyclotron frequency. Use of the EBW has particular advantages for measuring plasma magnetic field in plasmas with relatively high density and low magnetic field, as is the case in spherical tokamaks. This is the “over dense” condition of future efficient high beta magnetic confinement fusion reactors such as the spherical tokamak (STs). The overarching goal of the proposed project is to demonstrate beat wave EBW excitation in low density and temperature parameters that are directly scalable to over-dense spherical tokamaks such as NSTX or MAST, where measurement of magnetic field by electron cyclotron emission (ECE) is experimentally complex, and requires models of plasma profiles in order to interpret results. In such experiments, availability of point measurements of magnetic field strength provided by a beat-wave diagnostic would serve as a powerful independent check on the ECE results, serving to validate the assumed profiles and wave-propagation models used. The first step was to assemble/modify the target chamber in which the beat-wave plasma experiments will take place. Here, the target chamber of the CTIX compact-toroid (spheromak) injector was chosen. At the time that the proposal was submitted, the aim was to conduct the experiments in the CTIX experimental facilities located at the Livermore campus of UC Davis using microwave/mm-wave equipment brought from the main campus., However, by the time DoE funds for this project were awarded, the CTIX experimental facilities located at the Livermore campus of UC Davis were in the process of being shut down due to the closure of the Livermore campus by Lawrence Livermore National Laboratory (LLNL). Consequently, the equipment needed for this project was crated up and shipped from Livermore to Davis as part of the laboratory shutdown after space on campus was secured and College of Engineering approval secured. Once it arrived in Davis on December 20, 2018, work on the project began in earnest with the initial efforts directed toward developing a new plasma target source. Because of the shutdown of the CTIX facility, the scope of the project thus had to be reduced due to the abbreviated amount of time allocated to the project once the shutdown/move had been completed. Chamber modifications have been completed, including the fabrication and installation of a 38 GHz interferometer, as described in Sec. 3. The thermionic plasma source has been tested and characterized. The magnetic field required for the EBW experiments will be generated by permanent magnets mounted on plates positioned above and below the chamber. Magnetic simulations have been performed, and show that a 3-ring magnet arrangement on each plate, with a plate separation of 50 cm, is required to achieve the ~100 G uniform magnetic field in the interaction region that we feel is necessary to ensure good EBW coupling. We proved unable to successfully revive either of Prof. Luhmann’s two high power (≥20 kW) 94 GHz millimeter-wave sources after their extended time in storage and decided to postpone further efforts and concentrate on slightly lower power sources. Consequently, we were successful in reviving a ~2 kW extended klystron amplifier (EIK), which together with a ~2 kW extended interaction oscillator (EIO) will provide the two pump sources required for the project. Safe operation of the EIO requires extensive modifications to one of Prof. Luhmann’s pulse modulators - these modulator modifications are ongoing, which has unfortunately delayed the planned beat-wave experiments. Details of this work are provided in Sec. 3. Even though project funding has now ended, we still plan to proceed to complete the EIO pulse modulator modifications, install the two ~94 GHz pump wave sources on the target chamber, and conduct unmagnetized tests of the beat-wave concept which involve coupling to unmagnetized electron plasma waves rather than magnetized EBWs. Should these experiments prove successful, we plan to submit a follow-up proposal for additional funding which would allow us to undertake the magnetized EBW experiments.

42 ENGINEERING↗

Gas Bubble Evolution in Polycrystalline UMo Fuels Under Elastic-Plastic Deformation: A Phase-Field Model With Crystal-Plasticity

In monolithic UMo fuels, the interaction between the Al cladding and large gas bubble volumetric swelling causes both elastic-plastic and creep deformation. In this work, a phase-field model of gas bubble evolution in polycrystalline UMo under elastic-plastic deformation was developed for studying the dynamic interaction between evolving gas bubble/voids and deformation. A crystal plasticity model, which assumes that the plastic strain rate is proportional to resolved shear stresses of dislocation slip systems on their slip planes, was used to describe plastic deformation in polycrystalline UMo. Xe diffusion and gas bubble evolution are driven by the minimization of chemical and deformation energies in the phase-field model, while evolving gas bubble structure was used to update the mechanical properties in the crystal plasticity model. With the developed model, we simulated the effect of gas bubble structures (different volume fractions and internal gas pressures) on stress-strain curves and the effect of local stresses on gas bubble evolution. The results show that 1) the effective Young’s modulus and yield stress decrease with the increase of gas bubble volume fraction; 2) the hardening coefficient increases with the increase of gas bubble volume fraction, especially for gas bubbles with higher internal pressure; and 3) the pressure dependence of Xe thermodynamic and kinetic properties in addition to the local stress state determine gas bubble growth or shrinkage. The simulated results can serve as a guide to improve material property models for macroscale fuel performance modeling.

36 MATERIALS SCIENCE↗

Terra-Populus v0.1: A Python Library for LandScan High-Definition Population Analysis and Modeling

The terra-populus library is designed for use by the LandScan HD technical team, offering a streamlined set of tools for generating and updating LandScan HD datasets from foundational building-level data, referred to as 'molecules,' provided by the building-level attribution team. This document serves as the primary technical documentation for terra-populus. Version 0.1 of the library includes the core modeling components necessary for LandScan HD production. It enables the generation of the LandScan HD Baseline dataset as well as corresponding confidence measures for the occupancy rates used. Parameters have been included for incorporating damaged building indicators and changes in population, to faciliate the creation of rapid updates for LandScan HD. Future iterations of terra-populus will introduce tools for creating a confidence index, and quantifying and propagating uncertainty, facilitating the creation of probabilistic LandScan HD outputs. This report provides an overview of the tools available in the library and the corresponding code implementations. One of the key advancements implemented in terra-populus is a redefinition of the atomic modeling unit for LandScan HD. Traditionally, the LandScan HD vector analytical framework has generated population estimates at the building sub-component (molecule) level. However, terra-populus adopts a building-level modeling approach. This shift is an operational decision aimed at aligning LandScan HD outputs with confidence measures, which are computed and validated at the building level (confidence measures are not included in this version of terra-populus, aside from those associated with the occupancy rates). Additional advancements to the LandScan HD modeling, as implemented by terra-populus, include a minimum population value parameter and an auto assignment of building floor counts. The population minimum value was implemented to prevent buildings and subsequent LandScan HD pixels that contained small values that may not rasterize in production. An 'auto' value has been included as a method for dealing with buildings lacking floor count information, where it is the average floor count of all other buildings with a residential building use type tag. The logic behind this is to remain consistent with the current logic employed for dealing with building use type null instances, where a null use type is defaulted to residential since it is the most common building type. The auto logic is intended to apply the most common building floor count of the most common type of buildings. The tools provided in terra-populus represent a significant step forward in improving the efficiency, reproducibility, and transparency of the LandScan HD modeling process. As the library evolves, it will continue to serve as a foundational resource for high-resolution population modeling.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

The ezSpectra suite: An easy‐to‐use toolkit for spectroscopy modeling

Abstract A molecule's spectrum encodes information about its structure and electronic properties. It is a unique fingerprint that can serve as a molecular ID. Quantum chemistry calculations provide key ingredients for interpreting spectra, but modeling the spectra rarely ends there; it requires additional steps that entail combined treatments of electronic and nuclear degrees of freedom and account for specifics of the experimental setup (light energy, polarization, averaging over molecular orientations, temperature, etc.). This Software Focus article describes the ezSpectra suite, which currently comprises two stand‐alone open‐source codes: ezFCF and ezDyson . ezFCF calculates Franck–Condon factors, which yield vibrational progressions for polyatomic molecules, within the double‐harmonic approximation. ezDyson calculates absolute cross‐sections for photodetachment/photoionization processes and photoelectron angular distributions using Dyson orbitals computed by a quantum chemistry program. This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods Theoretical and Physical Chemistry > Spectroscopy Software > Simulation Methods

Gozem, Samer↗

NMME Monthly / Seasonal Forecasts for NASA SERVIR Applications Science

This work details use of the North American Multi-Model Ensemble (NMME) experimental forecasts as drivers for Decision Support Systems (DSSs) in the NASA / USAID initiative, SERVIR (a Spanish acronym meaning "to serve"). SERVIR integrates satellite observations, ground-based data and forecast models to monitor and forecast environmental changes and to improve response to natural disasters. Through the use of DSSs whose "front ends" are physically based models, the SERVIR activity provides a natural testbed to determine the extent to which NMME monthly to seasonal projections enable scientists, educators, project managers and policy implementers in developing countries to better use probabilistic outlooks of seasonal hydrologic anomalies in assessing agricultural / food security impacts, water availability, and risk to societal infrastructure. The multi-model NMME framework provides a "best practices" approach to probabilistic forecasting. The NMME forecasts are generated at resolution more coarse than that required to support DSS models; downscaling in both space and time is necessary. The methodology adopted here applied model output statistics where we use NMME ensemble monthly projections of sea-surface temperature (SST) and precipitation from 30 years of hindcasts with observations of precipitation and temperature for target regions. Since raw model forecasts are well-known to have structural biases, a cross-validated multivariate regression methodology (CCA) is used to link the model projected states as predictors to the predictands of the target region. The target regions include a number of basins in East and South Africa as well as the Ganges / Baramaputra / Meghna basin complex. The MOS approach used address spatial downscaling. Temporal disaggregation of monthly seasonal forecasts is achieved through use of a tercile bootstrapping approach. We interpret the results of these studies, the levels of skill by several metrics, and key uncertainties.

Robertson, Franklin R.↗

Impact of Defect Distributions on the Reliability of Structural Materials

The purpose of this study is to develop a probabilistic model to estimate the fatigue life of gas turbine disks made of powder superalloy metals. Such metals are known to have high strength and ductility under extremely high temperatures and known to have high fatigue lives. However, despite the high quality controls instituted during metal processing, sufficiently high numbers of ceramic inclusions materialize within the metal matrix creating weak links that lead to reduced fatigue lives. To study the effect of such ceramic inclusions on the fatigue life of gas turbine disks, an experimental program has been initiated at the NASA Glenn Research. In this phase of the study, a statistical analysis of the test results is undertaken which served as the basis of the probabilistic simulation model. The study is a part of the Crack Resistant Disk Materials task within the Ultra Safe Propulsion Project.

Ghosn, Michel J.↗

An experimental and kinetic modeling study of cyclopentane and dimethyl ether blends

Cyclopentane is a suitable naphthene, or cycloalkane, in a palette for multi-component gasoline surrogate fuels due to its presence in market fuels and its relevance to alkyl substituted cyclopentanes also present. However, the previous oxidation studies of cyclopentane have primarily focused on neat mixtures. Blending cyclopentane with dimethyl ether in this work therefore serves to inform our understanding of, and improve predictive models for, multi-component mixtures. Here, the auto-ignition of cyclopentane/dimethyl ether blends was studied in a high-pressure shock tube and in a rapid compression machine. A wide range of temperatures (650 – 1350 K) and elevated pressures of 20 and 40 bar were studied at equivalence ratios of 0.5, 1.0 and 2.0 in air for two blending ratios (30/70 and 70/30 mole% cyclopentane/di-methyl ether mixtures). A detailed kinetic model for cyclopentane was revised to capture the measured ignition delay times and apparent heat release rates in this study. Literature ignition delay time, jet-stirred reactor, and laminar burning velocity measurements of neat cyclopentane were used as additional validation. Improvements to the kinetic model were based on recent literature studies related to sub-models including cyclopentene and cyclopentadiene which allowed the removal of previous local rate-constant optimizations. Low temperature reactivity of cyclopentane was found to be controlled by the branching ratio between concerted elimination of HO͘ 2 and the strained formation of Q&#775OOH radicals in agreement with previous studies. In this study, the low branching ratio of Q&#775OOH formation increases the influence of a competing consumption pathway for cyclopentyl-peroxy (CPTO͘ 2 J) radicals. The sensitivity of the simulated ignition delay times to the formation of cyclopentyl hydroperoxide (CPTO 2 H), from CPTO͘ 2 J and HO͘ 2 , is discussed. The current model is used to analyze the influence of dimethyl ether on the reactivity of cyclopentane in the context of previous literature studies of dimethyl ether binary blends with ethanol and toluene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transport coefficient sensitivities in a semi-analytic model for magnetized liner inertial fusion

Performance of magnetized liner inertial fusion (MagLIF) experiments is highly dependent on transport processes including magnetized heat flows and magnetic flux losses. Magnetohydrodynamic simulations used to model these experiments require a choice of model for the transport coefficients, which are the constants of proportionality relating driving terms, such as temperature gradients and currents, to the associated heat and magnetic field transport. The coefficients have been the subject of repeated recalculation using various methods throughout the years. Using a semi-analytic MagLIF model, we compare models for the transport coefficients. The choice of model modifies magnetic-flux losses caused by the Nernst thermoelectric effect and thermal conduction losses. We present simulated results from parameter scans conducted in order to compare the effects of the different models on parameters of interest in MagLIF. In some regions of parameter space, discrepancies of up to 38% are found in integrated quantities like the fusion yield. These results may serve as a guide for experimental validation of the various models, particularly as laser preheat energies and initial axial field strengths are increased on MagLIF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation of the 2003 Foss Barge - Point Wells Oil Spill: A Comparison between BLOSOM and GNOME Oil Spill Models

The Department of Energy’s (DOE’s) National Energy Technology Laboratory’s (NETL’s) Blowout and Spill Occurrence Model (BLOSOM), and the National Oceanic and Atmospheric Administration’s (NOAA’s) General NOAA Operational Modeling Environment (GNOME) are compared. Increasingly complex simulations are used to assess similarities and differences between the two models’ components. The simulations presented here are forced by ocean currents from a Finite Volume Community Ocean Model (FVCOM) implementation that has excellent skill in representing tidal motion, and with observed wind data that compensates for a coarse vertical ocean model resolution. The comprehensive comparison between GNOME and BLOSOM presented here, should aid modelers in interpreting their results. Beyond many similarities, aspects where both models are distinct are highlighted. Some suggestions for improvement are included, e.g., the inclusion of temporal interpolation of the forcing fields (BLOSOM) or the inclusion of a deflection angle option when parameterizing wind-driven processes (GNOME). Overall, GNOME and BLOSOM perform similarly, and are found to be complementary oil spill models. This paper also sheds light on what drove the historical Point Wells spill, and serves the additional purpose of being a learning resource for those interested in oil spill modeling. The increasingly complex approach used for the comparison is also used, in parallel, to illustrate the approach an oil spill modeler would typically follow when trying to hindcast or forecast an oil spill, including detailed technical information on basic aspects, like choosing a computational time step. We discuss our successful hindcast of the 2003 Point Wells oil spill that, to our knowledge, had remained unexplained. The oil spill models’ solutions are compared to the historical Point Wells’ oil trajectory, in time and space, as determined from overflight information. Our hindcast broadly replicates the correct locations at the correct times, using accurate tide and wind forcing. While the choice of wind coefficient we use is unconventional, a simplified analytic model supported by observations, suggests that it is justified under this study’s circumstances. We highlight some of the key oceanographic findings as they may relate to other oil spills, and to the regional oceanography of the Salish Sea, including recommendations for future studies.

42 ENGINEERING↗

A system identification approach for non-intrusive reduced order modeling of radiation-induced photocurrents

In this study, development of compact photocurrent models is currently dominated by analytical techniques that rely on physical assumptions to render the governing equations solvable in a closed form. Violation of these assumptions can reduce the accuracy of the models and/or limit their scope. In this paper we show that system identification of nonlinear state-space systems can serve as an alternative numerical basis for non-intrusive reduced order modeling of photocurrent effects. To that end we develop a compact gray box photocurrent model (GBPM) by using a state-space representation with a low-dimensional latent state equation that mimics a mathematical model for the response of an idealized class of devices to ionizing radiation. In so doing we obtain a model that learns the dynamics of a quantity of interest directly from its measurements without requiring snapshots of the internal device state or its discretized model, and can be inferred from very small data sets. To demonstrate the approach we train the GBPM using a small experimental data set for a Z5236 Zener diode and a small synthetic data set obtained by simulating a synthetic pn-junction device. We then compare the GBPMs with black box models trained on the same data and show that performance of the latter is limited by the size of the data set, while the former are able to achieve excellent performance in both the reproductive and the predictive regimes.

97 MATHEMATICS AND COMPUTING↗