Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupled system”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Mooring System Cost Estimates for Wave Energy Farms in Shared Mooring Arrays

As wave energy converters (WECs) become more advanced and cost-efficient, so too must their mooring systems. A key question in the development plans for WECs is the cost of the mooring system, particularly for large wave farms. WEC devices deployed in a WEC farm array, where each device can be connected by shared mooring lines in various array layouts, have potential to reduce mooring system costs significantly. This paper presents the modeling and designing of mooring systems for large WEC arrays and the calculation of the cost of each mooring system to determine the change in cost as the number of WECs in a farm increases. A baseline mooring system for a single floating oscillating water column (OWC) WEC was developed for this analysis. The mooring system utilizes four anchored mooring lines connected to a square assembly of wire rope mooring lines, supported by four floating buoys, and attached to the floating WEC by four polyester rope mooring lines. This assembly, referred to as a floating cell, can be tiled to form various rectangular WEC arrays. The objective of this analysis is to determine how the mooring system cost changes as more WECs are added to an array layout, each with their own interconnected floating cell. To do this, complete mooring systems need to be designed for each WEC array layout. To narrow down the design space of a WEC array mooring system, a couple assumptions were made. It was assumed that the floating cell parameters of the baseline design were to stay constant across all floating cells in the WEC array. It was also assumed that the anchored mooring lines would be of the baseline configuration, a predominantly chain mooring line with a short section of polyester rope near the fairlead, and a drag-embedment anchor. These anchored lines were assumed to extend from the outer edges of the WEC array, inline with the headings of the wire rope mooring lines of the floating cells, or diagonal if extending from a corner of the WEC array. Full mooring systems were designed for 2xN and NxN WEC array layouts and efficiently simulated in the mooring dynamics simulation tool, MoorDyn, to ensure all dynamic constraints were met. The system costs were calculated and then refined by shortening unnecessary chain line lengths and reducing the chain diameters of the downstream anchored mooring lines. It was found that, in general, mooring system costs per WEC decrease when WECs are installed in an array. Compared to the baseline mooring system for a single WEC, the 2x3 mooring system array had the lowest mooring system cost per WEC, reducing the cost per WEC by 59%. The 3x3 and 4x4 array mooring systems also saw significant reductions in cost per WEC but had negligible cost savings between the two designs, primarily because the larger 4x4 mooring system requires larger chain diameters, which increases cost. These results provide an interesting glimpse into modeling, designing, and calculating the cost of mooring systems for large WEC arrays.

cost↗
DOE OSTI2026

Understanding Majorana braiding in superconductors with a fixed total number of particles

Mean-field one-dimensional topological superconductors host edge Majorana zero modes (MZMs) that encode a topologically protected ground-state degeneracy and enable robust braiding operations within this subspace. Mean-field states lack definite particle number and thus cannot represent isolated systems. Nevertheless, projecting them onto fixed particle number can yield good approximations to an isolated superconductor ground state. In earlier work [Sajith et al., Phys. Rev. B 109, 184509 (2024)] we showed that the projected Kitaev wave function of a single wire preserves some important mean-field features, such as the zero-energy spectral peaks near the wire edges. However, uniqueness of the fixed-number ground state does not allow for any nontrivial operations in the ground-state subspace. To overcome this limitation, here we consider the case of multiple wires with a conserved total charge, using the same approach. In the limit of vanishing interwire coupling, this system has a macroscopic ground-state degeneracy. We show how this degeneracy is resolved by coherent single-particle tunneling, and identify special many-body states which play the same role as the MZM parity states in the mean field. In this work, we demonstrate how braiding operations can be implemented and discuss both intrinsic and extrinsic limits on their fidelity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

The winter central Arctic surface energy budget: A model evaluation using observations from the MOSAiC campaign

This study evaluates the simulation of wintertime (15 October, 2019, to 15 March, 2020) statistics of the central Arctic near-surface atmosphere and surface energy budget observed during the MOSAiC campaign with short-term forecasts from 7 state-of-the-art operational and experimental forecast systems. Five of these systems are fully coupled ocean-sea ice-atmosphere models. Forecast systems need to simultaneously simulate the impact of radiative effects, turbulence, and precipitation processes on the surface energy budget and near-surface atmospheric conditions in order to produce useful forecasts of the Arctic system. This study focuses on processes unique to the Arctic, such as, the representation of liquid-bearing clouds at cold temperatures and the representation of a persistent stable boundary layer. It is found that contemporary models still struggle to maintain liquid water in clouds at cold temperatures. Given the simple balance between net longwave radiation, sensible heat flux, and conductive ground flux in the wintertime Arctic surface energy balance, a bias in one of these components manifests as a compensating bias in other terms. This study highlights the different manifestations of model bias and the potential implications on other terms. Three general types of challenges are found within the models evaluated: representing the radiative impact of clouds, representing the interaction of atmospheric heat fluxes with sub-surface fluxes (i.e., snow and ice properties), and representing the relationship between stability and turbulent heat fluxes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Strongly Coupled Plasmon and Phonon Polaritons as Seen by Photon and Electron Probes

The ability to control and modify infrared excitations in condensed matter is of both fundamental and applied interest. Here we explore a system supporting low-energy excitations, in particular, mid-infrared localized plasmon modes and phonon polaritons that are tuned to be strongly coupled. We study the coupled modes by using far-field infrared spectroscopy, state-of-the-art monochromated electron energy-loss spectroscopy, numerical simulations, and analytical modeling. Here we demonstrate that the electron probe facilitates a precise characterization of polaritons constituting the coupled system, and enables an active control over the coupling and the resulting sample response both in frequency and space. Although far-field optical spectra can be substantially different from near-field electron energy-loss spectra, we show that a direct comparison is possible via postprocessing and right positioning of the electron beam. The resulting spectra allow us to evaluate the key parameters of the coupled system, such as the coupling strength, which we demonstrate to be probe independent. Our work establishes a rigorous description of the spectral features observed in light- and localized electron-based spectroscopies, which can be extended to the analysis of analogous optical systems with applications in heat management and electromagnetic field concentration or nanofocusing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Facilitating better and faster simulations of aerosol-cloud interactions in Earth system models

Focal Area(s): 1. Predictive modeling through the use of AI techniques and AI-derived model components; the use of AI and other tools to design a prediction system comprising a hierarchy of models. 2. Insight gleaned from complex data (both observed and simulated) using AI, big data analytics, and other advanced methods, including explainable AI and physics- or knowledge-guided AI. Science Challenge: One major challenge that Earth system models (ESMs) face in providing credible prediction of the Earth system and its water cycle characteristics (e.g., mean state, variability, and extreme events) is to accurately simulate aerosol-cloud interactions (ACI). The physical, chemical, and dynamical processes affecting ACI are extremely complex and they range from nanoscale to planetary scale. In each model development cycle, scientists spend significant efforts investigating model deficiencies and uncertainties associated with aerosols (e.g., emissions, chemical processes, aerosol microphysics, and transport) and clouds (e.g., macrophysics, microphysics, turbulence, and large-scale circulation) in order to develop improved treatments. However, despite decades of active research, ACI is still a major source of uncertainty in climate projections, even though great progress has been made. Specific scientific challenges include: (i) Parameterizations are developed based on limited data; (ii) The complexity of a parameterization required for accurate predictions is not understood; (iii) Incomplete and unknown physics leads to errors in the fully coupled Earth system; and (iv) Complex physics is computationally too expensive to employ in ESMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Cascading economic losses from port disruptions under capacity constrained multimodal freight networks

This study quantifies how throughput disruptions at major seaports cascade through capacity-constrained multimodal freight networks and interregional production systems. We couple an agent-based model (ABM) multimodal freight simulation that resolves rerouting, terminal queueing, and inventory drawdown under binding modal and facility capacities with a multiregional output loss input-output (MRIIM) model that propagates realized delivery shortfalls across regions and sectors. The framework is demonstrated for the Port of Los Angeles using Freight Analysis Framework flows and Bureau of Economic Analysis input-output accounts and is evaluated over a 52-week horizon under deterministic sector targeted shocks and stochastic disruption realizations with uncertain severity and duration. Results indicate nonlinear amplification: realized national losses concentrate in manufacturing and transportation/warehousing even when exogenous port shocks are dispersed, suggesting that congestion spillback and limited short-run substitution can dominate the initial shock allocation. We further evaluate a tabular reinforcement-learning (Q-learning) intervention layer that selects among a small set of implementable system level levers (truck-to-rail and truck-to-barge shift settings) without overriding shipper routing, finding that such interventions reduce total losses for moderate disruptions but yield diminishing returns once substitute modes approach capacity. By linking operational freight behavior to system wide impacts under uncertainty, the proposed ABM-MRIIM pipeline provides a reusable workflow for port disruption stress testing, identification of structurally critical sectors/corridors, and evaluation of resilience interventions under realistic capacity limits.

42 ENGINEERING↗
DOE OSTI2020

Novel Fault Location Method for Power Systems Based on Attention Mechanism and Double Structure GRU Neural Network

Fault location is one of the most essential techniques to maintain the stable operation of power systems. A fast and accurate fault location allows operators to restore power grids faster and avoid economic losses. Conventional methods rely on expert knowledge to extract the necessary features (e.g. DWT, DFT). For large systems, more coupling effects of transmission lines require more complex feature engineering, and incomplete features can easily introduce large errors. To overcome this, a deep learning approach without manual feature extraction is introduced to the fault location model under big data application. Towards this end, in the proposed method, the attention mechanism, the Bi-GRU and a dual structure network are applied to analyze the current data from different perspectives. Complete information for the fault features is extracted for the fault location. Based on a time series model and benefit from the ability to internally acquire the information architecture of faulty line, the established model is adaptive to the power grids with very complex topologies. Simulation results indicate that the proposed double-structure model reduces the maximum error and is less affected by noise. In comparison with different structures and different models, the proposed method shows better performance in IEEE 39-bus system.

42 ENGINEERING↗
DOE OSTI2023

Simulations of ocean deoxygenation in the historical era: insights from forced and coupled models

Ocean deoxygenation due to anthropogenic warming represents a major threat to marine ecosystems and fisheries. Challenges remain in simulating the modern observed changes in the dissolved oxygen (O 2 ). Here, we present an analysis of upper ocean (0-700m) deoxygenation in recent decades from a suite of the Coupled Model Intercomparison Project phase 6 (CMIP6) ocean biogeochemical simulations. The physics and biogeochemical simulations include both ocean-only (the Ocean Model Intercomparison Project Phase 1 and 2, OMIP1 and OMIP2) and coupled Earth system (CMIP6 Historical) configurations. We examine simulated changes in the O 2 inventory and ocean heat content (OHC) over the past 5 decades across models. The models simulate spatially divergent evolution of O 2 trends over the past 5 decades. The trend (multi-model mean and spread) for upper ocean global O 2 inventory for each of the MIP simulations over the past 5 decades is 0.03 ± 0.39×10 14 [mol/decade] for OMIP1, -0.37 ± 0.15×10 14 [mol/decade] for OMIP2, and -1.06 ± 0.68×10 14 [mol/decade] for CMIP6 Historical, respectively. The trend in the upper ocean global O2 inventory for the latest observations based on the World Ocean Database 2018 is -0.98×10 14 [mol/decade], in line with the CMIP6 Historical multi-model mean, though this recent observations-based trend estimate is weaker than previously reported trends. A comparison across ocean-only simulations from OMIP1 and OMIP2 suggests that differences in atmospheric forcing such as surface wind explain the simulated divergence across configurations in O 2 inventory changes. Additionally, a comparison of coupled model simulations from the CMIP6 Historical configuration indicates that differences in background mean states due to differences in spin-up duration and equilibrium states result in substantial differences in the climate change response of O 2 . Finally, we discuss gaps and uncertainties in both ocean biogeochemical simulations and observations and explore possible future coordinated ocean biogeochemistry simulations to fill in gaps and unravel the mechanisms controlling the O 2 changes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Network-Cognizant Time-Coupled Aggregate Flexibility of Distribution Systems Under Uncertainties

Increasing integration of distributed energy resources (DERs) within distribution feeders provides unprecedented flexibility at the distribution-transmission interconnection. To exploit this flexibility and to use the capacity potential of aggregate DERs, feasible substation power injection trajectories need to be efficiently characterized. This paper provides an ellipsoidal inner approximation of the set of feasible power injection trajectories at the substation such that for any point in the set, there exists a feasible disaggregation strategy of DERs for any load uncertainty realization. The problem is formulated as one of finding the robust maximum volume ellipsoid inside the flexibility region under uncertainty. Though the problem is NP-hard even in the deterministic case, this paper derives novel approximations of the resulting adaptive robust optimization problem based on optimal second-stage policies. The proposed approach yields less conservative flexibility characterization than existing flexibility region approximation formulations. The efficacy of the proposed method is demonstrated on a realistic distribution feeder.

aggregates↗
DOE OSTI2023

Quantum information and quantum simulation of neutrino physics

In extreme astrophysical environments such as core-collapse supernovae and binary neutron star mergers, neutrinos play a major role in driving various dynamical and microphysical phenomena, such as baryonic matter outflows, the synthesis of heavy elements, and the supernova explosion mechanism itself. The interactions of neutrinos with matter in these environments are flavor-specific, which makes it of paramount importance to understand the flavor evolution of neutrinos. Flavor evolution in these environments can be a highly nontrivial problem thanks to a multitude of collective effects in flavor space, arising due to neutrino-neutrino ($v$-$v$) interactions in regions with high neutrino densities. A neutrino ensemble undergoing flavor oscillations under the influence of significant $v$-$v$ interactions is somewhat analogous to a system of coupled spins with long-range interactions among themselves and with an external field (‘long-range’ in momentum-space in the case of neutrinos). As a result, it becomes pertinent to consider whether these interactions can give rise to significant quantum correlations among the interacting neutrinos, and whether these correlations have any consequences for the flavor evolution of the ensemble. In particular, one may seek to utilize concepts and tools from quantum information science and quantum computing to deepen our understanding of these phenomena. In this article, we attempt to summarize recent work in this field. In conclusion, we also present some new results in a three-flavor setting, considering complex initial states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

STITCHES: creating new scenarios of climate model output by stitching together pieces of existing simulations

Abstract. Climate model output emulation has long been attempted to support impact research, mainly to fill in gaps in the scenario space. Given the computational cost of running coupled earth system models (ESMs), which are usually the domain of supercomputers and require on the order of days to weeks to complete a century-long simulation, only a handful of different scenarios are usually chosen to externally force ESM simulations. An effective emulator, able to run on standard computers in times of the order of minutes rather than days could therefore be used to derive climate information under scenarios that were not run by ESMs. Lately, the necessity of accounting for internal variability has also made the availability of initial-condition ensembles, under a specific scenario, important, further increasing the computational demand. At least so far, emulators have been limited to simplified ESM-like output, either seasonal, annual, or decadal averages of basic quantities, like temperature and precipitation, often emulated independently of one another. With this work, we propose a more comprehensive solution to ESM output emulation. Our emulator, STITCHES, uses existing archives of earth system models' (ESMs) scenario experiments to construct ESM-like output under new scenarios or enrich existing initial-condition ensembles, which is what other emulators also aim to do. Importantly, however, STITCHES' output has the same characteristics of the ESM output it sets out to emulate: multivariate, spatially resolved, and high frequency, representing both the forced component and the internal variability around it. STITCHES extends the idea of time sampling – according to which climate outcomes are stratified by the global warming level at which they manifest themselves, irrespective of the scenario and time at which they occur – to the construction of a continuous history of ESM-like output over the whole 21st century, consistent with a 21st-century trajectory of global surface air temperature (GSAT) derived from the scenario that has been chosen as the target of the emulation. STITCHES does so by first splitting the target GSAT trajectory into decade-long windows, then matching each window in turn to a decade-long window within an existing model simulation from the available scenario runs according to its proximity to the target in absolute size of the temperature anomaly and its rate of change. A look-up table is therefore created of a sequence of existing experiment–time-window combinations that, when stitched together, create a GSAT trajectory “similar” to the target. Importantly, we can then stitch together much more than GSAT from these windows, i.e., any output that the ESM has saved for these existing experiment–time-window combinations, at any frequency and spatial scale available in its archive. We show that the stitching does not introduce artifacts in the great majority of cases (we look at temperature and precipitation at monthly frequency and on the native grid of the ESM and at an index of ENSO activity, the Southern Oscillation Index). This is true even if the criteria for the identification of the decades to be stitched together are chosen to work for a smoothed time series of annual GSAT, a result we expect given the larger amount of noise affecting most other variables at finer spatial scales and higher frequencies, which therefore are more “forgiving” of the stitching. We successfully test the method's performance over many ESMs and scenarios. Only a few exceptions surface, but these less-than-optimal outcomes are always associated with a scarcity of the archived simulations from which we can gather the decade-long windows that form the building blocks of the emulated time series. In the great majority of cases, STITCHES' performance is satisfactory according to metrics that reward consistency in trends, interannual and inter-ensemble variance, and autocorrelation structure of the time series stitched together. The method therefore can be used to create ESM-like output according to new scenarios, on the basis of a trajectory of GSAT produced according to that scenario, which could be easily obtained by a simple climate model. It can also be used to increase the size of existing initial-condition ensembles. There are aspects of our emulator that will immediately disqualify it for specific applications, like when climate information is needed whose characteristics result from accumulated quantities over windows of times longer than those used as pieces by STITCHES, droughts longer than a decade for example. But for many applications, we argue that a stitched product can satisfy the climate information needs of impact researchers. STITCHES cannot emulate ESM output from scenarios that result in GSAT trajectories outside of the envelope available in the archive, nor can it emulate trajectories with shapes different from existing ones (overshoots with negative derivative, for example). Therefore, the size and characteristics of the available archives of ESM output are the principal limitations for STITCHES' deployment. Thus, we argue for the possibility of designing scenario experiments within, for example, the next phase of the Coupled Model Intercomparison Project according to new principles, relieved of the need to produce a number of similar trajectories that vary only in radiative forcing strength but more strategically covering the space of temperature anomalies and rates of change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Large-Area Lasing in Nanoscale Complex Media: The Critical Role of Local Dielectric Environment

Controlling how electromagnetic waves interact with complex media is critical for applications in imaging and focusing. Such lightwave interactions with complex media can lead to dramatic optical effects like lasing. While much work in random lasing focus on understanding how gain and scattering co-operatively generate lasing, little work has focused on how to manipulate the lasing threshold without modifying the structural disorder. Here, a simple, mostly unexplored, strategy is demonstrated that employs atomic layer deposition (ALD) to tune the local near-field environment while preserving the underpinning disorder—controlling lasing in a nanoscale complex medium on a large scale (>cm 2 ). The nanoscale complex medium is a quasi-2D system of coupled zinc oxide nanospheres with overall thickness deep in the sub-wavelength regime (≈λ/4). Near-ultraviolet femtosecond spectroscopy probes the broadband response of the gain nanomaterial, details how ALD process fundamentally modifies the fast-picosecond and slow-nanosecond carrier dynamics, and informs on the relevant timescales critical for lasing. Full-field electromagnetic simulations provide critical insights about how near-field dielectric environment modifies the nanostructure's scattering cross-section, which ultimately results in enhanced lasing. Importantly, these results highlight a simple path to control how electromagnetic waves interact in a complex medium, a key step toward large-scale implementation of complex lasers.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dispersion-Engineered Deep Sub-Wavelength Plasmonic Metasurfaces for Broadband Seira Applications

Surface enhanced infrared absorption (SEIRA) spectroscopy is a powerful tool in which plasmonically enhanced electromagnetic fields provide high-sensitivity molecular detection. Most SEIRA platforms operate at a single resonant frequency, which must be tuned to match that of the target molecule, and commonly rely on time-consuming lithographic techniques. This study presents a high-throughput and cost-effective plasmonic metasurface for broadband, tunable, and strong infrared signal enhancement. The platform is built upon the principle of dispersion-engineered plasmonic Fabry–Pérot (FP) nanocavity arrays. It offers 1) tight squeezing of infrared (IR) photons into deep sub-wavelength nano-volumes and 2) spectrally tunable near-field enhancements of up to ≈10 6 , two to three orders of magnitude higher than most optical metasurface systems. By coupling multilayer nano-thin film deposition and nanoskiving fabrication techniques, the dispersive FP metasurfaces can be rapidly and reproducibly constructed in a scalable and lithography-free manner. Using IR spectroscopy, the selective and sensitive label-free detection of a molecular monolayer is achieved at a range of frequencies. An enhancement factor of nearly 10 5 is measured at the carbonyl (C = O) vibrational marker band of the molecule. The confluence of high field enhancement, broadband plasmonic response, and facile fabrication makes this metasurface a promising platform for SEIRA spectroscopy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Time‐Dependent Cation Selectivity of Titanium Carbide MXene in Aqueous Solution

Abstract Electrochemical ion separation is a promising technology to recover valuable ionic species from water. Pseudocapacitive materials, especially 2D materials, are up‐and‐coming electrodes for electrochemical ion separation. For implementation, it is essential to understand the interplay of the intrinsic preference of a specific ion (by charge/size), kinetic ion preference (by mobility), and crystal structure changes. Ti 3 C 2 T z MXene is chosen here to investigate its selective behavior toward alkali and alkaline earth cations. Utilizing an online inductively coupled plasma system, it is found that Ti 3 C 2 T z shows a time‐dependent selectivity feature. In the early stage of charging (up to about 50 min), K + is preferred, while ultimately Ca 2+ and Mg 2+ uptake dominate; this unique phenomenon is related to dehydration energy barriers and the ion exchange effect between divalent and monovalent cations. Given the wide variety of MXenes, this work opens the door to a new avenue where selective ion‐separation with MXene can be further engineered and optimized.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

A non‐intrusive domain‐decomposition model reduction method for linear steady‐state partial differential equations with random coefficients

Abstract Domain decomposition methods have been proved to be an effective strategy to reduce the dimension of parametric partial differential equations (PDEs). However, existing domain decomposition methods for parametric PDEs are usually intrusive, which means domain decomposition based solvers need to be implemented from scratch for each target parametric PDE. To address this issue, we develop a new non‐intrusive domain‐decomposition model reduction method for linear steady‐state PDEs with random‐field coefficients. As a variant of our previous work by Mu and Zhang, the new method only needs access to the final linear system, that is, the global stiffness matrix and the right hand side, of a deterministic PDE solver, in order to build a domain‐decomposition‐based reduced model without intrusive implementation from scratch. The key idea is to remove the interface condition between sub‐domains and rely on the correlation between columns of the linear system to couple the sub‐domains. The non‐intrusive feature enables the applicability of the proposed method to a broader class of uncertainty quantification problems, where many legacy codes/solvers can be fully reused by our method. Two numerical examples including diffusion equations with random diffusivity and convection‐dominated transport with random velocity, are provided to demonstrate the effectiveness and efficiency of our method.

Zhang, Guannan↗
DOE OSTI2025

Central Compact Finite‐Difference Scheme With High Spectral Resolution for KdV Equation

This work presents a combination of cell‐node and cell‐centered compact finite difference scheme for the approximation of third derivatives involved in Korteweg–de Vries (KdV) equations. This approach employs a half‐shifted derivative construction at cell centers, avoiding the need for compact interpolation, thereby removing transfer errors; hence, it improves spectral resolution and maintains high‐order accuracy. Fourier analysis is performed to show the spectral properties of the proposed formulation, which provides higher spectral resolutions as compared to node‐based compact schemes. A filtering strategy is incorporated to suppress high‐frequency oscillations without compromising the accuracy of the numerical scheme, and the total variation diminishing Runge Kutta (TVDRK3) method is applied for time integration. Numerical experiments on linear, nonlinear, and coupled KdV systems are conducted, and a comparative analysis with cell‐node compact schemes confirms that the proposed scheme consistently reduces errors by up to an order of magnitude and achieves high spectral resolution properties.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

The holar wind

String theory in AdS 3 with purely NS-NS fluxes and vanishing RR moduli has a continuum of winding string excitations in radial plane wave states. BTZ black holes can emit such strings, which then flow out toward the AdS 3 boundary as a stream of massive quanta, and form a black hole analogue of the solar wind. The winding string sector thus provides a decay channel for the black hole to evaporate without having either to couple the system to an external reservoir or to match the AdS 3 throat onto an asymptotically flat region. We compute the emission amplitude of this “holar wind” in the semi-classical approximation, and consider the associated version of the black hole information paradox.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Stability Analysis of Coupled Advection-Diffusion Models with Bulk Interface Condition

Numerical stability is of critical importance in general circulation models because it affects the design of algorithms, time to solution, and computational costs associated with the simulations, which are very expensive in practice. In this paper we extend the stability analysis for ocean-atmosphere coupling proposed in [Zhang et al., J. Sci. Comput. 84, 44(2020)] to a more realistic model that includes horizontal advection. We analyze various time-stepping strategies. We find that advection has a stabilizing effect in scenarios common to climate models when bulk interface condition and explicit flux coupling are used. We also show that our method can be used to study the stability impact of advection for other interface conditions such as Dirichlet-Neumann conditions.

97 MATHEMATICS AND COMPUTING↗