Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonlinear experimental dynamics”

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

Droplet formation simulation using mixed finite elements

Droplet formation happens in finite time due to the surface tension force. The linear stability analysis is useful to estimate the size of a droplet but fails to approximate the shape of the droplet. This is due to a highly nonlinear flow description near the point where the first pinch-off happens. A one-dimensional axisymmetric mathematical model was first developed by Eggers and Dupont [“Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] using asymptotic analysis. This asymptotic approach to the Navier–Stokes equations leads to a universal scaling explaining the self-similar nature of the solution. Numerical models for the one-dimensional model were developed using the finite difference [Eggers and Dupont, “Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] and finite element method [Ambravaneswaran et al., “Drop formation from a capillary tube: Comparison of one-dimensional and two-dimensional analyses and occurrence of satellite drops,” Phys. Fluids 14, 2606–2621 (2002)]. The focus of this study is to provide a robust computational model for one-dimensional axisymmetric droplet formation using the Portable, Extensible Toolkit for Scientific Computation. The code is verified using the Method of Manufactured Solutions and validated using previous experimental studies done by Zhang and Basaran [“An experimental study of dynamics of drop formation,” Phys. Fluids 7, 1184–1203 (1995)]. The present model is used for simulating pendant drops of water, glycerol, and paraffin wax, with an aspiration of extending the application to simulate more complex pinch-off phenomena.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1990

Internal rotor friction instability

The analytical developments and experimental investigations performed in assessing the effect of internal friction on rotor systems dynamic performance are documented. Analytical component models for axial splines, Curvic splines, and interference fit joints commonly found in modern high speed turbomachinery were developed. Rotor systems operating above a bending critical speed were shown to exhibit unstable subsynchronous vibrations at the first natural frequency. The effect of speed, bearing stiffness, joint stiffness, external damping, torque, and coefficient of friction, was evaluated. Testing included material coefficient of friction evaluations, component joint quantity and form of damping determinations, and rotordynamic stability assessments. Under conditions similar to those in the SSME turbopumps, material interfaces experienced a coefficient of friction of approx. 0.2 for lubricated and 0.8 for unlubricated conditions. The damping observed in the component joints displayed nearly linear behavior with increasing amplitude. Thus, the measured damping, as a function of amplitude, is not represented by either linear or Coulomb friction damper models. Rotordynamic testing of an axial spline joint under 5000 in.-lb of static torque, demonstrated the presence of an extremely severe instability when the rotor was operated above its first flexible natural frequency. The presence of this instability was predicted by nonlinear rotordynamic time-transient analysis using the nonlinear component model developed under this program. Corresponding rotordynamic testing of a shaft with an interference fit joint demonstrated the presence of subsynchronous vibrations at the first natural frequency. While subsynchronous vibrations were observed, they were bounded and significantly lower in amplitude than the synchronous vibrations.

Walton, J.↗
DOE OSTI2020

Residue Specific Characterization of Electrostatics, Conformational Heterogeneity, and Dynamics of Electron Transfer Proteins with Linear and Nonlinear Infrared Spectroscopy. Final Report

The goal of this research project was to further our understanding of biological electron transfer through the development and application of new methods for studying proteins and their complexes with high spatial and temporal detail. The approach was based on incorporation of reporter bonds with distinct vibrational frequencies at specific locations in proteins and their characterization through the interaction of the bonds with pulses of infrared light. The project focused on the proteins plastocyanin (Pc) and cytochrome (cyt) f, two electron transfer partners in photosynthesis. The first thrust of the project was the incorporation of carbon deuterium bonds as non-perturbative probes at a ligand to the copper ion of Pc and their use to detect changes in the bonding due to the number of electrons at the metal and how the bonding changes upon association with cyt f. Toward developing more advanced approaches to characterize proteins, we incorporated a stronger reporter, the cyano group of cyanophenylalanine, at several locations in Pc and measured the frequency dynamics in the absence and presence of cyt f. The acquired data improved the current model for the protein complex. Notably, the data indicated that the majority of the population is involved in a loosely bound complex, providing insight into the nature of the ensembles formed by electron transfer protein complexes that is not available via other approaches. In addition to learning about the protein complex, we made strides toward developing and improving the new experimental approach to study highly dynamic behaviors of proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Angle-resolved nonresonant two-photon single ionization of argon using 9.3-eV photons produced via high-order harmonic generation

Here we present an experimental study on the photoionization dynamics of nonresonant one-color two-photon single valence ionization of neutral argon atoms. Using 9.3-eV photons produced via high harmonic generation and a three-dimensional momentum imaging spectrometer, we detect the photoelectrons and ions produced from nonresonant two-photon ionization in coincidence. Photoionization from the $3p$ orbital produces a photoelectron scattering wave function with $p$- and $f$-partial-wave components, which interfere and result in a photoelectron angular distribution with peak amplitude perpendicular to the vacuum ultraviolet polarization. The comparison between the present results and the two previous sets of theoretical calculations [Phys. Rev. A 44, 324 (1991), and J. Phys. B: At. Mol. Phys. 16, 2737 (1983)] indicates that electron-electron correlation contributes appreciably to the two-photon ionization dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Spatially Resolved Near Field Spectroscopy of Vibrational Polaritons at the Small N Limit

Vibrational polaritons, which have been primarily studied in Fabry–Pérot cavities with a large number of molecules (N ~ 10 6 –10 10 ) coupled to the resonator mode, exhibit various experimentally observed effects on chemical reactions. However, the exact mechanism is elusively understood from the theoretical side, as the large number of molecules involved in an experimental strong coupling condition cannot be represented completely in simulations. This discrepancy between theory and experiment arises from computational descriptions of polariton systems typically being limited to only a few molecules, thus failing to represent the experimental conditions adequately. To address this mismatch, we used surface phonon polariton (SPhP) resonators as an alternative platform for vibrational strong coupling. SPhPs exhibit strong electromagnetic confinement on the surface and thus allow for coupling to a small number of molecules. As a result, this platform can enhance nonlinearity and slow down relaxation to the dark modes. In this study, we fabricated a pillar-shaped quartz resonator and then coated it with a thin layer of cobalt phthalocyanine (CoPc). By employing scattering-type scanning near-field optical microscopy (s-SNOM), we spatially investigated the dependency of vibrational strong coupling on the spatially varying electromagnetic field strength and demonstrated strong coupling with 38,000 molecules only–reaching to the small N limit. Through s-SNOM analysis, we found that strong coupling was observed primarily on the edge of the quartz pillar and the apex of the s-SNOM tip, where the maximum field enhancement occurs. In contrast, a weak resonance signal and lack of coupling were observed closer to the center of the pillar. This work demonstrates the importance of spatially resolved polariton systems in nanophotonic platforms and lays a foundation to explore polariton chemistry and chemical dynamics at the small N limit–one step closer to reconcile with high-level quantum calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2011

Surface Instability of Liquid Propellant under Vertical Oscillatory Forcing

Fluid motion in a fuel tank produced during thrust oscillations can circulate sub-cooled hydrogen near the liquid-vapor interface resulting in increased condensation and ullage pressure collapse. The first objective of this study is to validate the capabilities of a Computational Fluid Dynamics (CFD) tool, CFD-ACE+, in modeling the fundamental interface transition physics occurring at the propellant surface. The second objective is to use the tool to assess the effects of thrust oscillations on surface dynamics. Our technical approach is to first verify the CFD code against known theoretical solutions, and then validate against existing experiments for small scale tanks and a range of transition regimes. A 2D axisymmetric, multi-phase model of gases, liquids, and solids is used to verify that CFD-ACE+ is capable of modeling fluid-structure interaction and system resonance in a typical thrust oscillation environment. Then, the 3D mode is studied with an assumed oscillatory body force to simulate the thrust oscillating effect. The study showed that CFD modeling can capture all of the transition physics from solid body motion to standing surface wave and to droplet ejection from liquid-gas interface. Unlike the analytical solutions established during the 1960 s, CFD modeling is not limited to the small amplitude regime. It can extend solutions to the nonlinear regime to determine the amplitude of surface waves after the onset of instability. The present simulation also demonstrated consistent trends from numerical experiments through variation of physical properties from low viscous fluid to high viscous fluids, and through variation of geometry and input forcing functions. A comparison of surface wave patterns under various forcing frequencies and amplitudes showed good agreement with experimental observations. It is concluded that thrust oscillations can cause droplet formation at the interface, which results in increased surface area and enhanced heat transfer between the liquid and gas phases as the ejected droplets travel well into the warmer gas region.

Yang, H. Q.↗
DOE OSTI2025

Pressure-induced structural and dielectric changes in liquid water at room temperature

Understanding the pressure-dependent dielectric properties of water is crucial for a wide range of scientific and practical applications. In this study, we employ a deep neural network trained on density functional theory data to investigate the dielectric properties of liquid water at room temperature across a pressure range of 0.1–1000 MPa. We observe a nonlinear increase in the static dielectric constant ɛ 0 with increasing pressure, a trend that is qualitatively consistent with experimental observations. This increase in ɛ 0 is primarily attributed to the increase in water density under compression, which enhances collective dipole fluctuations within the hydrogen-bonding network as well as the dielectric response. Despite the increase in ɛ 0 , our results reveal a decrease in the Kirkwood correlation factor G K with increasing pressure. Furthermore, this decrease in G K is attributed to pressure-induced structural distortions in the hydrogen-bonding network, which weaken dipolar correlations by disrupting the ideal tetrahedral arrangement of water molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2004

Cross-Scale Coupling in the Inner Magnetosphere

Magnetosphere-ionosphere (MI) coupling has interested scientists for decades and, in spite of experimental and theoretical research efforts, is still one of the least well-known dynamic processes in space plasma. The reason for this is that the numerous physical processes associated with MI coupling occur over multiple spatial lengths and temporal scales. One typical example of MI coupling is small- and large-scale ring current (RC) electrodynamic coupling. In this talk, we will address the two primary issues of RC electrodynamic coupling: (1) RC self-consistent coupling with electromagnetic ion cyclotron (EMIC) waves (small-scale electrodynamic coupling) and (2) RC self-consistent MI coupling that includes calculation of the magnetospheric electric field (large-scale electrodynamic coupling). We also will emphasize the role of the heavy ions in the number of wave-particle interaction magnetospheric processes. In particular, we will discuss some of the experimental and theoretical studies that have investigated the role of the heavy ions (mainly He(+) and O(+)) in generation and propagation of electromagnetic ion cyclotron waves and their contribution to the heating of magnetospheric electrons and ions. The more recent studies have also shown that the heavy ions can greatly contribute to a generation of lower hybrid waves, ring current precipitation phenomena, and the overall energy redistribution in the inner magnetosphere. Using newly developed 2.5-dimensional particle-in-cell simulations, we study the energization and nonlinear coupling of different plasma waves in the presence of the heavy ions. We have shown that the high frequency wave modes critically depend on the heavy ion density and irrespective of the driven wave modes, both the light and heavy ions undergo significant transverse acceleration. But for the large heavy-ion densities, even the electrons are significantly accelerated in the parallel direction by the waves below the LH frequency.

Khazanov, G. V.↗
NASA NTRS2004

An experimental study of the Richtmyer-Meshkov instability in microgravity

Richtmyer-Meshkov (RM) instability occurs when a planar interface separating two fluids of different density is impulsively accelerated in the direction of its normal. It is one of the most fundamental fluid instabilities and is of importance to the fields of astrophysics and inertial confinement fusion. Because RM instability experiments are normally carried out in shock tubes, where the generation of a sharp, well-controlled interface between gases is difficult, there is a scarcity of good experimental results. The experiments presented here use a novel technique that circumvents many of the experimental difficulties that have previously limited the study of RM instability in shock tubes. In these experiments, the instability is generated incompressibly, by bouncing a rectangular tank containing two liquids off of a fixed spring. These experiments, which utilize PLIF flow visualization, yield time-motion image sequences of the nonlinear development and transition to turbulence of the instability that are of a quality unattainable in shock tube experiments. Measurements obtained from these images, therefore, provide benchmark data for the evaluation of nonlinear models for the late-time growth of the instability. Because the run time in these experiments is limited, new experiments in the NASA Glenn 2.2 second drop tower, capable of achieving longer run times, are currently under way.

Rheology/methods↗
DOE OSTI2023

Microscopic activated dynamics theory of the shear rheology and stress overshoot in ultradense glass-forming fluids and colloidal suspensions

We formulate a particle and force level, activated dynamics-based statistical mechanical theory for the continuous startup nonlinear shear rheology of ultradense glass-forming hard sphere fluids and colloidal suspensions in the context of the elastically collective nonlinear Langevin equation approach and a generalized Maxwell model constitutive equation. Activated structural relaxation is described as a coupled local-nonlocal event involving caging and longer range collective elasticity which controls the characteristic stress relaxation time. Theoretical predictions for the deformation-induced enhancement of mobility, the onset of relaxation acceleration at remarkably low values of stress, strain, or shear rate, apparent power law thinning of the steady-state structural relaxation time and viscosity, a nonvanishing activation barrier in the shear thinning regime, an apparent Herschel–Buckley form of the shear rate dependence of the steady-state shear stress, exponential growth of different measures of a yield or flow stress with packing fraction, and reduced fragility and dynamic heterogeneity under deformation were previously shown to be in good agreement with experiments. The central new question we address here is the defining feature of the transient response—the stress overshoot. In contrast to the steady-state flow regime, understanding the transient response requires an explicit treatment of the coupled nonequilibrium evolution of structure, elastic modulus, and stress relaxation time. We formulate a new quantitative model for this aspect in a physically motivated and computationally tractable manner. Theoretical predictions for the stress overshoot are shown to be in good agreement with experimental observations in the metastable ultradense regime of hard sphere colloidal suspensions as a function of shear rate and packing fraction, and accounting for deformation-assisted activated motion appears to be crucial for both the transient and steady-state responses.

Mechanics↗
DOE OSTI2024

Validation of a Hybrid Domain Overlapping Coupling Between SAM and CFD Against the TALL-3D Transients

The System Thermal Hydraulics (STH) code SAM has been coupled to the Computational Fluid Dynamics (CFD) code Simcenter STAR-CCM+ utilizing a hybrid domain overlapping method with an explicit coupling in time. The coupling aims to extend the STH code’s applicability to scenarios where local momentum and energy transfers are important yet difficult for STH codes to capture, such as three-dimensional (3D) mixing. The coupling method’s numerical stability was verified in the past against two closed-loop configurations, and it was validated against a double T-junction experiment with 3D scalar mixing. In the present work, the coupling method is validated against the TALL-3D STH/CFD coupling benchmark facility. TALL-3D is a three-legged, liquid-metal facility with a large, pool-type enclosure (test section) that exhibits 3D flow effects to be modeled by a CFD code. The rest of the system exhibits approximately 1D behavior well-predicted by an STH code. First, the present STAR-CCM+ CFD model of the 3D test section is validated against experimental data. Then, the SAM-STARCCM+ coupled model is validated against six different TALL-3D steady states, including SAM standalone model results for comparison. Lastly, the SAM-STARCCM+ coupled model is validated against two TALL-3D transients, one exhibiting flow reversal in the test section and one exhibiting nonlinear, Limit Cycle Oscillations (LCO). For the first transient, the SAM-STARCCM+ coupled model properly predicts an increase in the test section’s inlet temperature during flow reversal, and this is not predicted by the SAM standalone model. Following flow reversal, the SAM-STARCCM+ coupled model better-predicts the initial flow recovery and following oscillations as the system approaches a final natural circulation state. For the second transient, no true final steady state is observed due to LCO. Neither the SAM-STARCCM+ coupled model nor the SAM standalone model can perfectly capture the experiment’s changing oscillation frequency during the transient. However, the SAM-STARCCM+ coupled model does reproduce the oscillatory feedback observed in the system. This is a significant achievement as the SAM-STARCCM+ coupled model only uses an explicit coupling in time, as opposed to a semi-implicit coupling. In comparison, previous STH/CFD coupling efforts of the TALL-3D facility required semi-implicit coupling to obtain similar results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

3D simulations of vertical displacement events in tokamaks: A benchmark of M3D-C 1 , NIMROD, and JOREK

In recent years, the nonlinear 3D magnetohydrodynamic codes JOREK, M3D-C1, and NIMROD developed the capability of modeling realistic 3D vertical displacement events (VDEs) including resistive walls. In this paper, a comprehensive 3D VDE benchmark is presented between these state-of-the-art codes. The simulated case is based on an experimental NSTX plasma but with a simplified rectangular wall. There are differences between the physics models and numerical methods, and the VDE evolution leads to sensitivities on the initial conditions that cannot be avoided as can be done in edge localized modes (ELM) and sawtooth simulations (due to the non-cyclical nature of VDEs). Nonetheless, the comparison serves to quantify the level of agreement in the relevant quantities used to characterize disruptions, such as the 3D wall forces and energy decay. The results bring confidence regarding the use of the mentioned codes for disruption studies, and they distinguish aspects that are specific to the models used (e.g., reduced vs full MHD models). Here, the simulations show important 3D fea- tures for a NSTX plasma, such as the self-consistent evolution of the halo current and the origin of the wall forces. In contrast to other reduced MHD models based on an ordering in the aspect ratio, the ansatz-based JOREK reduced MHD model allows capturing many aspects of the 3D dynamics even in the spherical tokamak limit considered here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2003

Airloads and Wake Geometry Calculations for an Isolated Tiltrotor Model in a Wind Tunnel

Th tiltrotor aircraft configuration has the potential to revolutionize air transportation by providing an economical combination of vertical take-off and landing capability with efficient, high-speed cruise flight. To achieve this potential it is necessary to have validated analytical tools that will support future tiltrotor aircraft development. These analytical tools must calculate tiltrotor aeromechanical behavior, including performance, structural loads, vibration, and aeroelastic stability, with an accuracy established by correlation with measured tiltrotor data. For many years such correlation has been performed for helicopter rotors (rotors designed for edgewise flight), but correlation activities for tiltrotors have been limited, in part by the absence of appropriate measured data. The recent test of the Tilt Rotor Aeroacoustic Model (TRAM) with a single, U4-scale V-22 rotor in the German-Dutch Wind Tunnel (DNW) now provides an extensive set of aeroacoustic, performance, and structural loads data. This paper will present calculations of airloads, wake geometry, and performance, including correlation with TRAM DNW measurements. The calculations were obtained using CAMRAD II, which is a modern rotorcraft comprehensive analysis, with advanced models intended for application to tiltrotor aircraft as well as helicopters. Comprehensive analyses have received extensive correlation with performance and loads measurements on helicopter rotors. The proposed paper is part of an initial effort to perform an equally extensive correlation with tiltrotor data. The correlation will establish the level of predictive capability achievable with current technology; identify the limitations of the current aerodynamic, wake, and structural models of tiltrotors; and lead to recommendations for research to extend tiltrotor aeromechanics analysis capability. The purpose of the Tilt Rotor Aeroacoustic Model (TRAM) experimental project is to provide data necessary to validate tiltrotor performance and aeroacoustic prediction methodologies and to investigate and demonstrate advanced civil tiltrotor technologies. The TRAM project is a key part of the NASA Short Haul Civil Tiltrotor (SHCT) project. The SHCT project is an element of the Aviation Systems Capacity Initiative within NASA. In April-May 1998 the TRAM was tested in the isolated rotor configuration at the Large Low-speed Facility of the German-Dutch Wind Tunnels (DNW). A preparatory test was conducted in December 1997. These tests were the first comprehensive aeroacoustic test for a tiltrotor, including not only noise and performance data, but airload and wake measurements as well. The TRAM can also be tested in a fill-span configuration, incorporating both rotors Lnd a fuselage model. The wind tunnel installation of the TRAM isolated rotor is shown. The rotor tested in the DNW was a 1/4-scale (9.5 ft diameter) model of the right-hand V-22 proprotor. The rotor and nacelle assembly was attached to an acoustically-treated, isolated rotor test stand through a mechanical pivot (the nacelle conversion axis). The TRAM was analyzed using the rotorcraft comprehensive analysis CAMRAD II. CAMRAD II is an aeromechanical analysis of helicopters and rotorcraft that incorporates a combination of advanced technologies, including multibody dynamics, nonlinear finite elements, and rotorcraft aerodynamics. The trim task finds the equilibrium solution (constant or periodic) for a steady state operating condition, in this case a rotor operating in a wind tunnel. For wind tunnel operation, the thrust and flapping are trimmed to target values. The aerodynamic model includes a wake analysis to calculate the rotor nonuniform induced-velocities, using a free wake geometry. The paper will present the results of CAMRAD II calculations compared to the TRAM DNW measurements for hover performance, helicopter mode performance, and helicopter mode airloads. An example of the hover performance results, comparing both mearements and calculations for the JVX (large scale) and TRAM (small scale) rotors, is shown. An example of the helicopter mode performance, showing the influence of the aerodynamic model (particularly the stall delay model) on the calculated power, induced power, and profile power is also shown. An example of the helicopter mode airloads, showing the influence of various wake and aerodynamic models on the calculations, is shown. Good correlation with measured airloads is obtained using the multiple-trailer wake model. The paper will present additional results, and describe and discuss the aerodynamic behavior in detail.

Johnson, Wayne↗
DOE OSTI2025

Extreme-ultraviolet spatiotemporal vortices via high harmonic generation

Spatiotemporal optical vortices (STOVs) are space–time structured light pulses with a unique topology that couples spatial and temporal domains and carry transverse orbital angular momentum (OAM). Up to now, their generation has been limited to the visible and infrared regions of the spectrum. During the last decade, it was shown that through the process of high-order harmonic generation, it is possible to upconvert spatial optical vortices that carry longitudinal OAM from the near-infrared into the extreme-ultraviolet (EUV), thereby producing vortices with distinct femtosecond and attosecond structure. In this work, we demonstrate theoretically and experimentally the generation of EUV spatiotemporal and spatiospectral vortices using near-infrared STOV driving laser pulses. Here, we use analytical expressions for focused STOVs to perform macroscopic calculations of high-order harmonic generation that are directly compared to the experimental results. As STOV beams are not eigenmodes of propagation, we characterize the highly charged EUV STOVs in both the near and far fields to show that they represent conjugated spatiotemporal and spatiospectral vortex pairs. Our work provides high-frequency light beams topologically coupled at the nanometre/attosecond scales domains with transverse OAM that could be suitable to explore electronic dynamics in magnetic materials, chiral media and nanostructures.

High-harmonic generation↗
DOE OSTI2022

Isotherm Modeling and Techno-Economic Analysis of Contactor Technologies for New Tetraamine-Appended MOF for NGCC Applications

Recently, a family of tetraamine-functionalized metal-organic frameworks (MOF) has been reported as promising sorbent materials for capturing CO2 from flue gas conditions relevant to natural gas combined cycle (NGCC) applications. The main advantages of these materials are their two-step cooperative CO2 adsorption, which gives rise to unusual two step-shaped CO2 adsorption profiles and their high thermal stability. This work presents the modelling of the two-transition isotherm of the tetraamine-appended MOF, N,N'-bis(3-aminopropyl)-1,4-diaminobutane (3-4-3)-appended Mg2(dobpdc), and the techno-economic analysis (TEA) of carbon capture processes utilizing this sorbent. Due to the unusual isotherm shapes of the experimental CO2 adsorption data for tetraamine-appended Mg2(dobpdc) and the strong nonlinearity of CO2 loading with respect to temperature and pressure, we tested two different models which use logistic functions for representing the different isotherm behaviors in the different pressure ranges. The first model uses the quadratic isotherm model in the low-pressure region, the Langmuir isotherm model in the middle pressure range, and the dual site Langmuir isotherm model in the high-pressure range. To model the transition between regions we used the arctangent functions independent of temperature and the thermal effect was accounted by using the Clausius-Clapeyron relation]. The second model is an extension of the weighted dual-site Langmuir isotherm model presented by Hughes et al.. In this extended model, the dual-site Langmuir isotherm is employed in the three transition regions, using temperature-dependent logistic functions to activate or deactivate the isotherm model in the low, middle, and high-pressure ranges. Both models fit the experimental data quite well with root mean squared errors (RMSE’s) of 0.41 and 0.17 for model 1 and model 2, respectively. Since model 2 resulted in a lower RMSE, it was leveraged for the development of the gas/solid contactor models used by the TEA. Specifically, two different contactor models, an axial-flow fixed bed and moving bed contactor, were developed as part of this work. These models are dynamic, pressure-driven, and consist of mass, energy, and momentum conservation equations. A kinetic model was also developed by performing parameter estimation using experimental fixed bed breakthrough data. These models are then used to simulate CO2 capture processes from the flue gas generated from a ~600 gross MW NGCC power plant. A cost model was developed which considers the capital cost of the reactors and the significant operating costs such as steam and electricity. Using NETL’s Framework for the Optimization and Quantification of Uncertainty of Uncertainty and Surrogates tool (FOQUS), which has the capability of linking models built using numerous modelling platforms with derivative-free optimization solvers, a techno-economic optimization of the carbon capture processes was performed which minimizes the cost of capture.

Caballero, Daison↗
DOE OSTI2025

Floquet engineering of polaritonic amplification in dispersive photonic time crystals

In this study, we investigate the dynamics of dispersive photonic time crystals (PTCs) and their potential applications for controlling light-matter interaction. Using the Lorentz-Drude as dispersive material model, we analyze the effects of periodic modulation of dispersion parameters through both theoretical and numerical analysis, revealing the emergence of hybrid bandgaps from the interaction of polaritonic branches with distinctive characteristics. Our study demonstrates that dispersive PTCs offer novel excitation channels and amplification possibilities, that require lower modulation frequencies compared to nondispersive systems thus alleviating experimental challenges for the realization of PTCs in the optical regime. These findings pave the way for advancements in polaritonic lasing and resonant Raman scattering.

Classical optics↗
DOE OSTI2023

Time-dependent SOLPS-ITER simulations of the tokamak plasma boundary for model predictive control using SINDy *

Abstract Time-dependent SOLPS-ITER simulations have been used to identify reduced models with the sparse identification of nonlinear dynamics (SINDy) method and develop model-predictive control of the boundary plasma state using main ion gas puff actuation. A series of gas actuation sequences are input into SOLPS-ITER to produce a dynamic response in upstream and divertor plasma quantities. The SINDy method is applied to identify reduced linear and nonlinear models for the electron density at the outboard midplane n e , s e p O M P and the electron temperature at the outer divertor T e , s e p d i v . Note that T e , s e p d i v is not necessarily the peak value of T e along the divertor. The identified reduced models are interpretable by construction (i.e. not black box), and have the form of coupled ordinary differential equations. Despite significant noise in T e , s e p d i v , the reduced models can be used to predict the response over a range of actuation levels to a maximum deviation of 0.5% in n e , s e p O M P and 5%–10% in T e , s e p d i v for the cases considered. Model retraining using time history data triggered by a preset error threshold is also demonstrated. A model predictive control strategy for nonlinear models is developed and used to perform feedback control of a SOLPS-ITER simulation to produce a setpoint trajectory in n e , s e p O M P using the integrated plasma simulator framework. The developed techniques are general and can be applied to time-dependent data from other boundary simulations or experimental data. Ongoing work is extending the approach to model identification and control for divertor detachment, which will present transient nonlinear behavior from impurity seeding, including realistic latency and synthetic diagnostic signals derived from the full SOLPS-ITER output.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1996

Nonlinear Bubble Interactions in Acoustic Pressure Fields

The systems consisting of a two-phase mixture, as clouds of bubbles or drops, have shown many common features in their responses to different external force fields. One of particular interest is the effect of an unsteady pressure field applied to these systems, case in which the coupling of the vibrations induced in two neighboring components (two drops or two bubbles) may result in an interaction force between them. This behavior was explained by Bjerknes by postulating that every body that is moving in an accelerating fluid is subjected to a 'kinetic buoyancy' equal with the product of the acceleration of the fluid multiplied by the mass of the fluid displaced by the body. The external sound wave applied to a system of drops/bubbles triggers secondary sound waves from each component of the system. These secondary pressure fields integrated over the surface of the neighboring drop/bubble may result in a force additional to the effect of the primary sound wave on each component of the system. In certain conditions, the magnitude of these secondary forces may result in significant changes in the dynamics of each component, thus in the behavior of the entire system. In a system containing bubbles, the sound wave radiated by one bubble at the location of a neighboring one is dominated by the volume oscillation mode and its effects can be important for a large range of frequencies. The interaction forces in a system consisting of drops are much smaller than those consisting of bubbles. Therefore, as a first step towards the understanding of the drop-drop interaction subject to external pressure fluctuations, it is more convenient to study the bubble interactions. This paper presents experimental results and theoretical predictions concerning the interaction and the motion of two levitated air bubbles in water in the presence of an acoustic field at high frequencies (22-23 KHz).

Barbat, Tiberiu↗