SEARCH · Engineering Papers

Results for “Non linear 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.

SourceResearch connectionFrom yearThrough year

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

At least 523 records · Page 29

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Stability analysis of a one-dimensional multiphysics model of a molten salt fast reactor

Reactor designs with a fast neutron spectrum have recently moved into the research focus for possible Generation IV reactor concepts - a representative of which is the molten salt fast reactor (MSFR). In order to assess the safety and reliability of this reactor concept, it is necessary to obtain knowledge of the system behaviour. An important part of this research concerns the dynamical stability of the MSFR. Research regarding the MSFR stability so far has either been using linearized models for stability criteria or only been examining certain transients at fixed parameter values. This work delivers a comprehensive stability analysis for a non-linearized MSFR model, and a wide range of values for all relevant parameters. For this purpose, a one-dimensional MSFR model was set up, taking neutron kinetics and thermal hydraulics in the form of a system of differential equations into account. The stability of this model was investigated by means of the numerical tool MATCONT, which was used to monitor the evolution of a so-called fixed-point solution, here referring to the steady state at operating conditions of the system, while system parameters got varied. The MATCONT analyses showed no loss of stability for any of the considered parameter variations and no solutions that might co-exist in parallel with the steady-state fixed-point solution. Therefore these results indicate a stable fixed point to which all solution transients converge, making the considered MSFR model insusceptible to deviations from the equilibrium. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2025

Fused-ring isomerism modulates molecular packing and device performance in non-halogenated organic solar cells

Subtle changes in molecular backbone geometry impact intermolecular interactions and performance of organic solar cells. Here, three isomeric small-molecule acceptors (NaO1, NaO2, and NaO3) are investigated to reveal how different fused-ring configurations control molecular packing, electronic coupling, and film formation. Structural and spectroscopic analyses show that the linearly fused NaO1 forms a compact three-dimensional packing network with large and balanced electronic couplings (>24 meV) across multiple directions, while the more curved analogues exhibit excessive crystallization and phase segregation. In-situ optical measurements demonstrate that NaO1 promotes fast and continuous structural evolution during film formation, resulting in smooth morphology and homogeneous phase distribution. These structural and dynamic advantages facilitate efficient charge generation and transport, accompanied by reduced non-radiative energy loss, ultimately achieving an efficiency of 20.07% for non-halogenated ternary devices. Our findings highlight how fused-ring isomerism decisively governs structure–packing–performance relationships in organic solar cells.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Deep Koopman operators for causal discovery

Causal discovery aims to identify cause-effect mechanisms for better scientific understanding, explainable decision-making, and more accurate modeling. Standard statistical frameworks, such as Granger causality, lack the ability to quantify causal relationships in nonlinear dynamics due to the presence of complex feedback mechanisms, timescale mixing, and nonstationarity. Thus, applying these methods to study causal dynamics in real-world systems, such as the Earth, is a major challenge. Addressing this shortcoming, we leverage deep learning and a Koopman operator-theoretic formalism to present a class of causal discovery algorithms. Kausal uses deep Koopman operator methods to approximate nonlinear dynamics in a linearized vector space in which traditional causal inference methods such as Granger causality can be more easily applied. Our idealized experiments demonstrate Kausal’s superior ability in discovering and characterizing causal signals compared to existing deep learning and non-deep learning state-of-the-art approaches. Finally, the successful identification of major El Niño and La Niña events in observations showcases Kausal’s skill to handle real-world applications.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

The nonlinear behavior of generic tail fins for small wind turbines

This paper describes analysis and measurements of the yaw response of tail fins for small wind turbines. It is based on an extension of unsteady slender body theory (USBT) to cover non-slender fins and high angles of incidence, both of which make the theory nonlinear. We provide three main additions to the substantial literature on linearized USBT for tail fins. First, USBT is extended to high angles by modeling the nonlinear vortex dynamics. Second, the restriction to slender bodies is removed by modeling the chordwise load variation. Third, we consider the effect of time-varying wind speed. Further, the extended theory is compared to wind tunnel measurements of the yaw behavior of delta, elliptical, and rectangular tail fins without a rotor and nacelle. The fins were released from initial yaw angles of -40° and -80°; the latter is of sufficient magnitude to show the importance of the nonlinear yaw dynamics. Generally good agreement was found between the theory and measurements, and the theory was shown to be more accurate than a “polar” or quasi-steady model which uses only the lift and drag of a delta planform. Of the three planforms, the rectangular one showed the lowest accuracy in terms of frequency but the damping was accurately predicted. Overall, the results demonstrate the importance of nonlinearity in the response of a yawing tail fin, particularly for the higher aspect ratio fins at large yaw angles.

17 WIND ENERGY↗

NASA NTRS1990

Spurious frequencies as a result of numerical boundary treatments

The stability theory for finite difference Initial Boundary-Value approximations to systems of hyperbolic partial differential equations states that the exclusion of eigenvalues and generalized eigenvalues is a sufficient condition for stability. The theory, however, does not discuss the nature of numerical approximations in the presence of such eigenvalues. In fact, as was shown previously, for the problem of vortex shedding by a 2-D cylinder in subsonic flow, stating boundary conditions in terms of the primitive (non-characteristic) variables may lead to such eigenvalues, causing perturbations that decay slowly in space and remain periodic time. Characteristic formulation of the boundary conditions avoided this problem. A more systematic study of the behavior of the (linearized) one-dimensional gas dynamic equations under various sets of oscillation-inducing legal boundary conditions is reported.

Abarbanel, Saul↗

NASA NTRS2012

Parametric Study of Decay of Homogeneous Isotropic Turbulence Using Large Eddy Simulation

Numerical simulations of decaying homogeneous isotropic turbulence are performed with both low-order and high-order spatial discretization schemes. The turbulent Mach and Reynolds numbers for the simulations are 0.2 and 250, respectively. For the low-order schemes we use either second-order central or third-order upwind biased differencing. For higher order approximations we apply weighted essentially non-oscillatory (WENO) schemes, both with linear and nonlinear weights. There are two objectives in this preliminary effort to investigate possible schemes for large eddy simulation (LES). One is to explore the capability of a widely used low-order computational fluid dynamics (CFD) code to perform LES computations. The other is to determine the effect of higher order accuracy (fifth, seventh, and ninth order) achieved with high-order upwind biased WENO-based schemes. Turbulence statistics, such as kinetic energy, dissipation, and skewness, along with the energy spectra from simulations of the decaying turbulence problem are used to assess and compare the various numerical schemes. In addition, results from the best performing schemes are compared with those from a spectral scheme. The effects of grid density, ranging from 32 cubed to 192 cubed, on the computations are also examined. The fifth-order WENO-based scheme is found to be too dissipative, especially on the coarser grids. However, with the seventh-order and ninth-order WENO-based schemes we observe a significant improvement in accuracy relative to the lower order LES schemes, as revealed by the computed peak in the energy dissipation and by the energy spectrum.

Swanson, R. C.↗

NASA NTRS2003

MHD Instability and Turbulence in the Tachocline

The focus of this project was to study the physical processes that govern tachocline dynamics and structure. Specific features explored included stratification, shear, waves, and toroidal and poloidal background fields. In order to address recent theoretical work on anisotropic mixing and dynamics in the tachocline, we were particularly interested in such anisotropic mixing for the specific tachocline processes studied. Transition to turbulence often shapes the largest-scale features that appear spontaneously in a flow during the development of turbulence. The resulting large-scale straining field can control the subsequent dynamics; therefore, anticipation of the large-scale straining field that results for individual realizations of the transition to turbulence can be important for subsequent dynamics, flow morphology, and transport characteristics. As a result, we paid particular attention to the development of turbulence in the stratified and sheared environment of the tachocline. This is complicated by the fact that the linearly stability of sheared MHD flows is non-self-adjoint, implying that normal asymptotic linear stability theory may not be relevant.

Werne, Joe↗

DOE OSTI2024

A prototype scintillator real‐time beam monitor for ultra‐high dose rate radiotherapy

Background: FLASH Radiotherapy (RT) is an emergent cancer RT modality where an entire therapeutic dose is delivered at more than 1000 times higher dose rate than conventional RT. For clinical trials to be conducted safely, a precise and fast beam monitor that can generate out-of-tolerance beam interrupts is required. This paper describes the overall concept and provides results from a prototype ultra-fast, scintillator-based beam monitor for both proton and electron beam FLASH applications. Purpose: A FLASH Beam Scintillator Monitor (FBSM) is being developed that employs a novel proprietary scintillator material. The FBSM has capabilities that conventional RT detector technologies are unable to simultaneously provide: (1) large area coverage; (2) a low mass profile; (3) a linear response over a broad dynamic range; (4) radiation hardness; (5) real-time analysis to provide an IEC-compliant fast beam-interrupt signal based on true two-dimensional beam imaging, radiation dosimetry and excellent spatial resolution. Methods: The FBSM uses a proprietary low mass, less than 0.5 mm water equivalent, non-hygroscopic, radiation tolerant scintillator material (designated HM: hybrid material) that is viewed by high frame rate CMOS cameras. Folded optics using mirrors enable a thin monitor profile of ∼10 cm. A field programmable gate array (FPGA) data acquisition system generates real-time analysis on a time scale appropriate to the FLASH RT beam modality: 100–1000 Hz for pulsed electrons and 10–20 kHz for quasi-continuous scanning proton pencil beams. An ion beam monitor served as the initial development platform for this work and was tested in low energy heavy-ion beams ( 86 Kr +26 and protons). A prototype FBSM was fabricated and then tested in various radiation beams that included FLASH level dose per pulse electron beams, and a hospital RT clinic with electron beams. Results: Results presented in this report include image quality, response linearity, radiation hardness, spatial resolution, and real-time data processing. Furthermore, the HM scintillator was found to be highly radiation damage resistant. It exhibited a small 0.025%/kGy signal decrease from a 216 kGy cumulative dose resulting from continuous exposure for 15 min at a FLASH compatible dose rate of 237 Gy/s. Measurements of the signal amplitude versus beam fluence demonstrate linear response of the FBSM at FLASH compatible dose rates of >40 Gy/s. Comparison with commercial Gafchromic film indicates that the FBSM produces a high resolution 2D beam image and can reproduce a nearly identical beam profile, including primary beam tails. The spatial resolution was measured at 35–40 µm. Tests of the firmware beta version show successful operation at 20 000 Hz frame rate or 50 µs/frame, where the real-time analysis of the beam parameters is achieved in less than 1 µs. Conclusions: The FBSM is designed to provide real-time beam profile monitoring over a large active area without significantly degrading the beam quality. A prototype device has been staged in particle beams at currents of single particles up to FLASH level dose rates, using both continuous ion beams and pulsed electron beams. Using a novel scintillator, beam profiling has been demonstrated for currents extending from single particles to 10 nA currents. Radiation damage is minimal and even under FLASH conditions would require ≥50 kGy of accumulated exposure in a single spot to result in a 1% decrease in signal output. Beam imaging is comparable to radiochromic films, and provides immediate images without hours of processing. Real-time data processing, taking less than 50 µs (combined data transfer and analysis times), has been implemented in firmware for 20 kHz frame rates for continuous proton beams.

2D beam imaging↗

DOE OSTI2020

Dynamically morphing microchannels in liquid crystal elastomer coatings containing disclinations

Liquid crystal elastomers (LCEs) hold a major promise as a versatile material platform for smart soft coatings, since their orientational order can be predesigned to program a desired dynamic profile. In this work, we introduce temperature-responsive dynamic coatings based on LCEs with arrays of singular defects-disclinations that run parallel to the surface. The disclinations form in response to antagonistic patterns of the molecular orientation at the top and bottom surfaces, imposed by the plasmonic mask photoalignment. Upon heating, an initially flat LCE coating develops linear microchannels located above each disclination. The stimulus that causes a non-flat profile of LCE coatings upon heating is the activation force induced by the gradients of molecular orientation around disclinations. To describe the formation of microchannels and their thermal response, we adopt a Frank-Oseen model of disclinations in patterned director field and propose a linear elasticity theory to connect the complex spatially-varying molecular orientation to the displacements of the LCE. The thermo-responsive surface profiles predicted by the theory and by the finite element modeling are in good agreement with the experimental data; in particular, higher gradients of molecular orientation produce a stronger modulation of the coating profile. The elastic theory and the finite element simulations allow us to estimate the material parameter that characterizes the elastomer coating’s response to the thermal activation. The disclinations containing LCEs show potential as soft dynamic coatings with a predesigned responsive surface profile.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Life-Cycle Assessment Integration into Scalable Open-Source Numerical Models (LiAISON) for Prospective Impact Analysis of Novel Technologies

Decarbonizing the industrial sector is a significant challenge in achieving a net-zero greenhouse gas (GHG) emissions economy by 2050 and the Paris Agreement, i.e., a global climate change mitigation target of achieving a maximum average temperature change potential of 1.5 Degrees Celsius or less by 2100 with respect to pre-industrial levels. In the United States (US), the industrial sector accounts for 23% of total GHG emissions and is home to a number of hard-to-electrify activities. The chemicals subsector has the single largest subsector emissions profile after direct emissions from fossil fuel combustion and leakage from fossil fuel distribution systems. Within the chemicals subsector, many processes depend on hydrogen or ammonia precursors. Decarbonizing these two commodities would contribute significantly to decarbonizing the industrial sector as hydrogen could also be used for low carbon steel production (e.g., hydrogen-based direct reduction of iron) and other industrial applications. Emerging technologies require the application of prospective life cycle assessment (LCA), which can account for technology (foreground) scaling and process improvements via learning-by-doing, among others. In many cases, the future system context (background) in which the technologies are assumed to operate in is equally relevant. Background scenarios generated by integrated assessment models (IAM) can coherently incorporate potential future dynamics of the energy-climate-human-land system. Further, IAM scenarios are harmonized across socioeconomic and climate change mitigation pathways, which facilitates the comparability of prospective LCAs using different IAMs. We introduce an open source prospective LCA framework, the Life-cycle Assessment Integration into Scalable Open-source Numerical models (LiAISON), to analyze the non-linear relationships between technology foreground and the future energy system background across a series of midpoint and resource use metrics. The integration of LCA and IAM data is achieved using prospective environmental Impact assessment (PREMISE). We showcase it by assessing two Power-to-Hydrogen (PtH2) processes, namely Solid Oxide Electrolysis (SOE) and Polymer Electrolyte Membrane Electrolysis (PEME). We compare the technologies to a baseline of hydrogen production via natural gas-based Steam Methane Reforming (SMR) in a US context of multiple energy system and climate change mitigation futures. Besides providing an analysis that specifies the LCA results ranges with temporal and geospatial explicitness across the two technologies, metrics, and impact assessment methods, this research also aims to establish a base framework that can be expanded to use other IAM generated scenarios and US open-source life cycle inventory (LCI) databases. We find that the temporal environmental performance of either technology or their difference to SMR is directly influenced by the underlying background dynamics. Additionally we compare our results by linking two other prospective models with LiAISON - GCAM (Global Change Assessment Model) and ReEDS (Regional Energy Deployment System) to analyze the effect of changing background scenarios using varying predictions in life cycle analysis.

decarbonizing↗

DOE OSTI2024

Life-Cycle Assessment Integration into Scalable Open-Source Numerical Models (LiAISON) for Analyzing Emerging Low-Carbon Technologies

Decarbonizing the industrial sector is a significant challenge in achieving a net-zero greenhouse gas (GHG) emissions economy by 2050 and the Paris Agreement, i.e., a global climate change mitigation target of achieving a maximum average temperature change potential of 1.5 degrees Celsius or less by 2100 with respect to pre-industrial levels. In the United States (US), the industrial sector accounts for 23% of total GHG emissions and is home to a number of hard-to-electrify activities. The chemicals subsector has the single largest subsector emissions profile after direct emissions from fossil fuel combustion and leakage from fossil fuel distribution systems. Within the chemicals subsector, many processes depend on hydrogen or ammonia precursors. Decarbonizing these two commodities would contribute significantly to decarbonizing the industrial sector as hydrogen could also be used for low carbon steel production (e.g., hydrogen-based direct reduction of iron) and other industrial applications. Emerging technologies require the application of prospective life cycle assessment (LCA), which can account for technology (foreground) scaling and process improvements via learning-by-doing, among others. In many cases, the future system context (background) in which the technologies are assumed to operate in is equally relevant. Background scenarios generated by integrated assessment models (IAM) can coherently incorporate potential future dynamics of the energy-climate-human-land system. Further, IAM scenarios are harmonized across socioeconomic and climate change mitigation pathways, which facilitates the comparability of prospective LCAs using different IAMs. We introduce an open source prospective LCA framework, the Life-cycle Assessment Integration into Scalable Open-source Numerical models (LiAISON), to analyze the non-linear relationships between technology foreground and the future energy system background across a series of midpoint and resource use metrics. The integration of LCA and IAM data is achieved using prospective environmental Impact assessment (PREMISE). We showcase it by assessing two Power-to-Hydrogen (PtH2) processes, namely Solid Oxide Electrolysis (SOE) and Polymer Electrolyte Membrane Electrolysis (PEME). We compare the technologies to a baseline of hydrogen production via natural gas-based Steam Methane Reforming (SMR) in a US context of multiple energy system and climate change mitigation futures. Besides providing an analysis that specifies the LCA results ranges with temporal and geospatial explicitness across the two technologies, metrics, and impact assessment methods, this research also aims to establish a base framework that can be expanded to use other IAM generated scenarios and US open-source life cycle inventory (LCI) databases. We find that the temporal environmental performance of either technology or their difference to SMR is directly influenced by the underlying background dynamics. Under baseline projections (i.e., no decarbonization goals), neither process reaches parity with the incumbent technology across several environmental metrics. Under the decarbonization scenarios, the underlying sectoral shifts result in declining impacts over time, compared to 2020 levels, except for metal depletion levels, which increase. The background shifts postulate a heavily decarbonized economy and energy system, which help technologies reach parity with SMR between 2040-2050 (RCP2.6) and 2030-2040 (RCP1.9) for global warming. Despite declines across several other metrics over time, neither PtH2 technology break even with SMR by 2100 besides for global warming.

decarbonizing↗

DOE OSTI2024

Life-Cycle Assessment Integration into Scalable Open-Source Numerical Models (LiAISON) for Analyzing Emerging Low-Carbon Technologies

Decarbonizing the industrial sector is a significant challenge in achieving a net-zero greenhouse gas (GHG) emissions economy by 2050 and the Paris Agreement, i.e., a global climate change mitigation target of achieving a maximum average temperature change potential of 1.5 degrees C or less by 2100 with respect to pre-industrial levels. In the United States (US), the industrial sector accounts for 23% of total GHG emissions and is home to a number of hard-to-electrify activities. The chemicals subsector has the single largest subsector emissions profile after direct emissions from fossil fuel combustion and leakage from fossil fuel distribution systems. Within the chemicals subsector, many processes depend on hydrogen or ammonia precursors. Decarbonizing these two commodities would contribute significantly to decarbonizing the industrial sector as hydrogen could also be used for low carbon steel production (e.g., hydrogen-based direct reduction of iron) and other industrial applications. Emerging technologies require the application of prospective life cycle assessment (LCA), which can account for technology (foreground) scaling and process improvements via learning-by-doing, among others. In many cases, the future system context (background) in which the technologies are assumed to operate in is equally relevant. Background scenarios generated by integrated assessment models (IAM) can coherently incorporate potential future dynamics of the energy-climate-human-land system. Further, IAM scenarios are harmonized across socioeconomic and climate change mitigation pathways, which facilitates the comparability of prospective LCAs using different IAMs. We introduce an open source prospective LCA framework, the Life-cycle Assessment Integration into Scalable Open-source Numerical models (LiAISON), to analyze the non-linear relationships between technology foreground and the future energy system background across a series of midpoint and resource use metrics The integration of LCA and IAM data is achieved using prospective environmental Impact assessment (PREMISE). We showcase it by assessing two Power-to-Hydrogen (PtH2) processes, namely Solid Oxide Electrolysis (SOE) and Polymer Electrolyte Membrane Electrolysis (PEME). We compare the technologies to a baseline of hydrogen production via natural gas-based Steam Methane Reforming (SMR) in a US context of multiple energy system and climate change mitigation futures. Besides providing an analysis that specifies the LCA results ranges with temporal and geospatial explicitness across the two technologies, metrics, and impact assessment methods, this research also aims to establish a base framework that can be expanded to use other IAM generated scenarios and US open-source life cycle inventory (LCI) databases. We find that the temporal environmental performance of either technology or their difference to SMR is directly influenced by the underlying background dynamics. Under baseline projections (i.e., no decarbonization goals), neither process reaches parity with the incumbent technology across several environmental metrics. Under the decarbonization scenarios, the underlying sectoral shifts result in declining impacts over time, compared to 2020 levels, except for metal depletion levels, which increase. The background shifts postulate a heavily decarbonized economy and energy system, which help technologies reach parity with SMR between 2040-2050 (RCP2.6) and 2030-2040 (RCP1.9) for global warming. Despite declines across several other metrics over time, neither PtH2 technology break even with SMR by 2100 besides for global warming.

decarbonizing↗

DOE OSTI2024

Towards Prospective LCA Using Life-Cycle Assessment Integration into Scalable Open-Source Numerical Models (LiAISON) Framework for Analyzing Emerging Low-Carbon Technologies

Decarbonizing the industrial sector is a significant challenge in achieving a net-zero greenhouse gas (GHG) emissions economy by 2050 and the Paris Agreement, i.e., a global climate change mitigation target of achieving a maximum average temperature change potential of 1.5 Degrees Celsius or less by 2100 with respect to pre-industrial levels. In the United States (US), the industrial sector accounts for 23% of total GHG emissions and is home to a number of hard-to-electrify activities. The chemicals subsector has the single largest subsector emissions profile after direct emissions from fossil fuel combustion and leakage from fossil fuel distribution systems. Within the chemicals subsector, many processes depend on hydrogen or ammonia precursors. Decarbonizing these two commodities would contribute significantly to decarbonizing the industrial sector as hydrogen could also be used for low carbon steel production (e.g., hydrogen-based direct reduction of iron) and other industrial applications. Emerging technologies require the application of prospective life cycle assessment (LCA), which can account for technology (foreground) scaling and process improvements via learning-by-doing, among others. In many cases, the future system context (background) in which the technologies are assumed to operate in is equally relevant. Background scenarios generated by integrated assessment models (IAM) can coherently incorporate potential future dynamics of the energy-climate-human-land system. Further, IAM scenarios are harmonized across socioeconomic and climate change mitigation pathways, which facilitates the comparability of prospective LCAs using different IAMs. We introduce an open source prospective LCA framework, the Life-cycle Assessment Integration into Scalable Open-source Numerical models (LiAISON), to analyze the non-linear relationships between technology foreground and the future energy system background across a series of midpoint and resource use metrics The integration of LCA and IAM data is achieved using prospective environmental Impact assessment (PREMISE). We showcase it by assessing two Power-to-Hydrogen (PtH2) processes, namely Solid Oxide Electrolysis (SOE) and Polymer Electrolyte Membrane Electrolysis (PEME). We compare the technologies to a baseline of hydrogen production via natural gas-based Steam Methane Reforming (SMR) in a US context of multiple energy system and climate change mitigation futures. Besides providing an analysis that specifies the LCA results ranges with temporal and geospatial explicitness across the two technologies, metrics, and impact assessment methods, this research also aims to establish a base framework that can be expanded to use other IAM generated scenarios and US open-source life cycle inventory (LCI) databases. We find that the temporal environmental performance of either technology or their difference to SMR is directly influenced by the underlying background dynamics. Additionally we compare our results by linking two other prospective models with LiAISON - GCAM(Global Change Assessment Model) and ReEDS (Regional Energy Deployment System) to analyze the effect of changing background scenarios using varying predictions in life cycle analysis.

emissions↗

DOE OSTI2026

Advances in 3D transient plasma dynamics and control through MHD and hybrid fluid-kinetic simulations with JOREK

Transient phenomena and their control are of high relevance in magnetic confinement fusion plasmas to guarantee a stable and safe plasma operation. Interpretative simulations can maximize the insights gained from experiments on present machines and predictive simulations can help in the preparation of design, mitigation techniques and operational scenarios for future devices. In this article, we provide an overview of recent advances and novel scientific results obtained with the 3D non-linear hybrid fluid-kinetic code JOREK, covering physics of plasma transients from the core to the scrape-off layer (SOL) both for tokamak and stellarator devices. Substantial progress was made in the physics understanding, model validation with experiments and experiment interpretation, thus, giving confidence for predictions to devices like DTT, ITER and DEMO. The topics addressed comprise a wide range: the edge physics of new operation scenarios and edge localized mode suppression; major disruptions with a focus on runaway electrons and vertical displacement events as well as disruption mitigation by shattered pellet injection; the physics mechanisms and operational limits of the flux pumping regime for sawtooth control; MHD limits of stellarators and work towards incorporating advanced edge/SOL/exhaust dynamics; continuing improvements of the code for more efficient hybrid simulations on conventional and accelerated high performance computing architectures.

disruptions↗

DOE OSTI2025

Ten-moment fluid modeling of the Weibel instability

We investigate the one-dimensional non-relativistic Weibel instability through the capture of anisotropic pressure tensor dynamics using an implicit 10-moment fluid model that employs the electromagnetic Darwin approximation. The results obtained from the 10-moment model are compared with an implicit particle-in-cell simulation. The linear growth rates obtained from the numerical simulations are in good agreement with the theoretical fluid and kinetic dispersion relations. The fluid dispersion relations are derived using Maxwell’s equations and the Darwin approximation. We also show that the magnetohydrodynamic approximation can be used to model the Weibel instability if one accounts for an anisotropic pressure tensor and unsteady terms in the generalised Ohm’s law. In addition, we develop a preliminary theory for the saturation magnetic field strength of the Weibel instability, showing good agreement with the numerical results.

Kuldinow, D. A. (ORCID:0000000319730196)↗

NASA NTRS2002

Modal and Impact Dynamics Analysis of an Aluminum Cylinder

This paper presents analyses for the modal characteristics and impact response of an all-aluminum cylinder. The analyses were performed in preparation for impact tests of the cylinder at The Impact Dynamics Research Facility (IDRF) at the NASA Langley Research Center. Mode shapes and frequencies were computed using NASTRAN and compared with existing experimental data to assess the overall accuracy of the mass and stiffness of the finite element model. A series of non-linear impact analyses were then performed using MSC Dytran in which the weight distribution on the floor and the impact velocity of the cylinder were varied. The effects of impact velocity and mass on the rebound and gross deformation of the cylinder were studied in this investigation.

Lessard, Wendy B.↗

NASA NTRS2016

Yield and Failure Behavior Investigated for Cross-linked Phenolic Resins Using Molecular Dynamics

Molecular dynamics simulations were conducted to fundamentally evaluate the yield and failure behavior of cross-linked phenolic resins at temperatures below the glass transition. Yield stress was investigated at various temperatures, strain rates, and degrees of cross-linking. The onset of non-linear behavior in the cross-linked phenolic structures was caused by localized irreversible molecular rearrangements through the rotation of methylene linkers followed by the formation or annihilation of neighboring hydrogen bonds. The yield stress results, with respect to temperature and strain rate, could be fit by existing models used to describe yield behavior of amorphous glasses. The degree of cross-linking only indirectly influences the maximum yield stress through its influence on the glass transition temperature (Tg), however there is a strong relationship between the degree of cross-linking and the failure mechanism. Low cross-linked samples were able to separate through void formation, whereas the highly cross-linked structures exhibited bond scission.

Stress-strain curves↗

DOE OSTI2022

The Mira–Titan Universe – IV. High-precision power spectrum emulation

Modern cosmological surveys are delivering data sets characterized by unprecedented quality and statistical completeness; this trend is expected to continue in the future as new ground- and space-based surveys come online. In order to maximally extract cosmological information from these observations, matching theoretical predictions are needed. At low redshifts, the surveys probe the non-linear regime of structure formation where cosmological simulations are the primary means of obtaining the required information. The computational cost of sufficiently resolved large-volume simulations makes it prohibitive to run very large ensembles. Nevertheless, precision emulators built on a tractable number of high-quality simulations can be used to build very fast prediction schemes to enable a variety of cosmological inference studies. We have recently introduced the Mira–Titan Universe simulation suite designed to construct emulators for a range of cosmological probes. This gravity-only set of simulations covers the standard six cosmological parameters {ω m , ω b , σ 8 , $h, n_s, w_0$} and, in addition, includes massive neutrinos and a dynamical dark energy equation of state {ω ν , $w_a$}. In this paper, we present the final emulator for the matter power spectrum based on 111 cosmological simulations, each covering a (2.1 Gpc) 3 volume and evolving 3200 3 particles. In this work, an additional set of 1776 lower resolution simulations and TimeRG perturbation theory results for the power spectrum are used to cover scales straddling the linear to mildly non-linear regimes (maximum wavenumber $\textit{k}$ = 5 Mpc –1 ). The emulator provides predictions at the 2–3 percent level of accuracy over a wide range of cosmological parameters and is publicly released as part of this paper.

79 ASTRONOMY AND ASTROPHYSICS↗