Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel simulation”

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

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

At least 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Integrating Marine Hydrokinetic and Offshore Wind Energy: A Review of Technologies, Deployment, and Challenges

Together, offshore wind (OSW) and marine hydrokinetic (MHK) technologies have vast potential to expand the world’s access to abundant energy resource. With more than 60 GW of offshore wind energy capacity and 527 MW of ocean energy deployed globally by 2023, there is a significant amount of available resources; however, technical and non-technical challenges prevent the combined large-scale deployment of these technologies. There is still a lack of research that provides a parallel review of both MHK and OSW technologies in order to better understand their synergistic working principles. This paper aimed to address that research gap by presenting a comprehensive side-by-side review of the worldwide technological landscape, global deployment trends, integration strategies, and modeling approaches for MHK and OSW. A particular focus has been given on analyzing existing modeling and simulation techniques, assessing integration and control strategies, and comparing technologies based on water depth. Furthermore, this study provides important insights into the readiness levels of both technologies by highlighting ongoing international projects. By addressing these issues, this review will give researchers and industry stakeholders an outline for assessing the maturity of OSW and MHK systems and facilitating their transition to large-scale, sustainable deployment.

16 - TIDAL AND WAVE POWER↗
DOE OSTI2022

Using a Grid-Forming Inverter to Stabilize a Low-Inertia Power System - Maui Hawaiian Island

As power systems around the world integrate greater amounts of wind and solar photovoltaic power, periods of very high instantaneous power shares of inverters, the primary interfacing technology for these generation sources, are complicating system stability and control. The contemporary, primary mode of inverter operation, grid-following, which explicitly assumes the presence of a local, stable voltage waveform, yields operational inadequacy at high instantaneous power shares potentially leading to instability due to the low-inertia conditions, as well as the correlated reduced voltage forming capacity on the respective system. Parallel connected grid-forming inverters, which directly regulate the local voltage, are a solution that is expected to bolster system stability and mitigate the shortcomings of the grid-following technology. In this paper, grid-following and two types of grid-forming inverter control, the traditional linear droop and the recently introduced nonlinear exponential droop (Droop-e), are simulated on a low-inertia, H = 0.48s, high inverter-based resource scenario, 97%, with a validated electromagentic transient domain model of the Hawaiian island of Maui power system. The benefit of a single grid-forming device over its grid-following counterpart is significant, both in terms frequency deviation and voltage stability. Further, the superiority of the Droop-e and the associated secondary power sharing control over both grid-following and linear droop grid-forming technologies is displayed, with improved nadir and rate of change of frequency over the linear droop control.

droop control↗
DOE OSTI2021

The Polarized Image of a Synchrotron-emitting Ring of Gas Orbiting a Black Hole

Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov and conservation of the Walker–Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The effect of a transducer’s spatial averaging on an elastodynamic guided wave’s wavenumber spectrum

Elastodynamic guided waves propagate in an elastic solid which makes it difficult, if not impractical, to place a receiving transducer in the direct path of the propagating wave (as one would for an acoustic wave in a fluid medium). Thus, to account for this, receiving transducers are often placed on the surface of the solid waveguide such that the transducer surface is parallel to the wave propagation direction. This transducer orientation introduces spatial averaging, which causes the received signal to have an altered signal amplitude and mode bias. We investigate both of these effects and present a simple model from which we derive a scaling-ratio expression that describes the effects of spatial averaging. Furthermore, we then test its performance using finite-element simulations that incorporate “real-world” assumptions (e.g., transient waves, piezoelectric effects, etc.). The results from the simulations demonstrate that the scaling-ratio can characterize the effects resulting from spatial averaging. The scaling-ratio expression will be particularly useful when designing experiments involving high frequency (small wavelength) guided waves. Lastly, the proposed scaling-ratio expression could be applied to other sensing methods, like Laser Doppler Vibrometry (LDV) or piezoelectric waveguides with inter-digitated transducers, because of its generality. However, the authors are careful to note that the scaling-ratio expression is not intended as a replacement of multi-physics analysis or good experimental design, and the effects of spatial averaging should be avoided whenever possible.

2DFFT↗
DOE OSTI2025

Interaction of Planar Brushes with End-Functionalized Carboxylic or Amine Groups: Water Structure, Hydrogen Bonds, and Electrostatic Correlations

In this work we present a detailed analysis of two parallel plates grafted with polyethylene glycol (PEG) terminated with carboxylic or amine groups, and discuss the cases of opposite and same charge in the presence of electrolytes. We analyze the screening of charged groups by electrolytes, the role of hydrogen bonds in enhancing attraction or repulsion, hydration forces, competitive binding with counterions, and the hydrogen bond network, crucial to stabilizing the system. Our simulations reveal that, for large plate separation, the system is well described as a “Stern layer”, characterized by strong counterion binding and a subsequent weak diffuse layer. For short plate separation, the system is characterized by strong correlations among end groups and counterions. We also discuss the relevance of charged carboxylic/amine groups in different problems with a particular emphasis on nanoparticle assembly.

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

Polarization-resolved core exciton dynamics in LiF using attosecond transient absorption spectroscopy

The ability to control absorption by modifying the polarization of light presents an exciting opportunity to experimentally determine the orbital alignment of absorption features. Here, attosecond extreme ultraviolet (XUV) transient absorption spectroscopy is used to investigate the polarization dependence of core exciton dynamics in LiF thin films at the Li + K edge. XUV pulses excite electrons from the Li 1s core level into the conduction band, allowing for the formation of a p-orbital-like core exciton, aligned along the XUV light polarization axis. A sub-5 fs near-infrared (NIR) probe pulse then arrives at variable time delays, perturbing the XUV-excited states and allowing the coherence decay of the core exciton to be mapped. The coherence lifetimes are found to be ≈ 2.4 ± 0.4 fs, which is attributed to a phonon-mediated dephasing mechanism as in previous core exciton studies. The differential absorption features are also shown to be sensitive to the relative polarization of the XUV and NIR fields. The parallel NIR probe induces couplings between the initial XUV-excited p-like bright exciton and s-like dark excitons. When crossed pump and probe polarizations are used, the coupling between the bright and dark states is no longer dipole-allowed, and the transient absorption signal associated with the coupling is suppressed by approximately 90%. This interpretation is supported by simulations of a few-level model system, as well as analysis of the calculated band structure. Here, the results indicate that laser polarization can serve as a powerful experimental tool for exploring the orbital alignment of core excitonic states in solid-state materials.

Attosecond laser spectroscopy↗
DOE OSTI2021

The discrete Green's function paradigm for two-way coupled Euler–Lagrange simulation

We outline a methodology for the simulation of two-way coupled particle-laden flows. The drag force that couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. In addition to hydrodynamics, the method can be extended to other physics present in particle-laden flows such as heat transfer and electromagnetism. The discrete Green's functions for the Navier–Stokes equations are obtained at low particle Reynolds number in a two-plane channel geometry. We perform verification at different Reynolds numbers for a particle settling under gravity parallel to a plane wall, for different wall-normal separations. Compared with other point-particle schemes, the Stokesian discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations. To account for degradation in accuracy away from the wall at finite Reynolds number, we extend the present methodology to an Oseen-like discrete Green's function. The extended discrete Green's function method is found to be accurate within 6% at all wall-normal separations for particle Reynolds numbers up to 24. Furthermore, the discrete Green's function approach is well suited to dilute systems with significant mass loading and this is highlighted by comparison against other Euler–Lagrange as well as particle-resolved simulations of gas–solid turbulent channel flow. Strong particle–turbulence coupling is observed in the form of turbulence modification and turbophoresis suppression, and these observations are placed in context of the different methods.

42 ENGINEERING↗
DOE OSTI2022

Anisotropic Thermally Conductive Composites Enabled by Acoustophoresis and Stereolithography

Opportunities to improve thermal management in electronic devices are currently hindered by processing constraints that limit thermal conductivity in polymer-matrix composites. Active patterning of filler particles is a promising route to improve conductivity while retaining processability by improving particle contact density and directing heat along optimized pathways. Here, we employ acoustic patterning to align and compact filler particles into stripes during stereolithographic 3D printing. This approach produces polymer-based composite materials with highly efficient embedded heat transport pathways which reach 95 vol% particle utilization (relative to the parallel conduction upper limit). These composites exhibit anisotropic thermal conductivity up to 300% higher than unpatterned composites, with in-plane anisotropy ratios of up to 350%. Combining this high conductivity with 3D printing enables materials with engineered heat networks that optimize transport from hot spots to heatsinks while maintaining low viscosity for fast particle patterning and for infiltration around electronic components. Finally, numerical simulations of acoustic assembly of particles with varied geometry, when compared to experimentally characterized particle packing, illuminate pathways for further improving conductivity by optimizing particle geometry for alignment and stacking of particles with maximum contact surface area.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The structure of 3-D collisional magnetized bow shocks in pulsed-power-driven plasma flows

Here we investigate three-dimensional (3-D) bow shocks in a highly collisional magnetized aluminium plasma, generated during the ablation phase of an exploding wire array on the MAGPIE facility (1.4 MA, 240 ns). Ablation of plasma from the wire array generates radially diverging, supersonic ($M_S$ $\sim 7$), super-Alfvénic ($M_A>1$) magnetized flows with frozen-in magnetic flux $(R_M\gg1)$. These flows collide with an inductive probe placed in the flow, which serves both as the obstacle that generates the magnetized bow shock, and as a diagnostic of the advected magnetic field. Laser interferometry along two orthogonal lines of sight is used to measure the line-integrated electron density. A detached bow shock forms ahead of the probe, with a larger opening angle in the plane parallel to the magnetic field than in the plane normal to it. Since the resistive diffusion length of the plasma is comparable to the probe size, the magnetic field decouples from the ion fluid at the shock front and generates a hydrodynamic shock, whose structure is determined by the sonic Mach number, rather than the magnetosonic Mach number of the flow. The 3-D simulations performed using the resistive magnetohydrodynamic (MHD) code Gorgon confirm this picture, but under-predict the anisotropy observed in the shape of the experimental bow shock, suggesting that non-MHD mechanisms may be important for modifying the shock structure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Integrating Advanced Modeling and Accelerated Testing for a Modernized Fuel Qualification Paradigm

With the increasing interest in sodium fast reactor technology, as seen by applications to the U.S. Nuclear Regulatory Commission for the OKLO Aurora plant, fuel testing for the TerraPower Traveling Wave Reactor, and the impending construction and startup of the versatile test reactor (VTR), a modernized, accelerated approach to fuel qualification is needed. To guide this effort, a Phenomena Identification Ranking Table–styled analysis was performed for a U-Pu-Zr sodium-free annular fuel system. This analysis evaluated a series of fuel design properties and parameters against their contributions to key fuel performance phenomena. The resulting priority parameters were then reviewed against existing modeling and experimental capabilities to support investigation of the highest-priority parameters. A pathway for qualification was then established using highthroughput, high-volume experiments from MiniFuel and FAST in parallel with advanced physics-based model development. This effort outlines how the first stages of qualification can be reduced from the typical 20+-year development cycle to 5 to 7 years by deploying accelerated irradiation testing platforms. As with any accelerated test, these methods are prototypic in some aspects and less so in others; however, by coupling with advanced fuel performance modeling and simulation capabilities, the larger space of irradiation parameters and material response provided offers advantages for the validation of physics-based models supporting the deployment of novel fuel designs. As a test case, this paper utilizes a proposed Mark II fuel system for the upcoming VTR. Thus, an accelerated qualification method can be tested for the development of MARK II driver fuel so that by the time of VTR startup, lead test assemblies for a Mark II fuel can be initiated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Multifunctional Heat Pump with Energy Storage Capability

Replacing gas furnaces with heat pumps on a large scale can pose challenges to the electrical grid. It will significantly increase the demand, strain the electrical grid, especially during peak usage times. To ensure that heat pumps contribute to reducing emissions, it's essential to increase the share of renewable energy in the grid. However, the supply of renewable energy is unstable. Large scale of electrification must address the concern of grid resilience to meet necessary heating/cooling demands and shift peak electric load. All these require future heat pumps having energy storage capability. These systems can store excess energy during low-demand periods and release it during high-demand periods, reducing stress on the grid. This paper will introduce an innovative configuration, control, and laboratory investigations of a multifunctional packaged heat pump. In addition to the original indoor and outdoor air coils, the heat pump has a parallel hydronic loop with a brazed plate heat exchanger (BHP). The hydronic loop circulates hot/chilled water to store heating/cooling energy in phase change material (PCM). It is capable of space cooling, space heating, cooling energy storage/defrosting, water heating/heating energy storage with outdoor air source or indoor air source, and cooling/heating energy free discharge. The laboratory studies and breadboard unit are presented. Additionally, building energy simulation demonstrates utility cost reduction potential via charging PCM panels/ceilings during off-peak hours and releasing the energy during peak hours to mitigate the peak power consumption.

Shen, Bo↗
DOE OSTI2023

Scalable DPG multigrid solver for Helmholtz problems: A study on convergence

This paper presents a scalable multigrid preconditioner targeting large-scale systems arising from discontinuous Petrov–Galerkin (DPG) discretizations of high-frequency wave operators. This work is built on previously developed multigrid preconditioning techniques of Petrides and Demkowicz (Comput. Math. Appl. 87 (2021) pp. 12–26) and extends the convergence results from $\mathscr{O}$(10 7 ) degrees of freedom (DOFs) to $\mathscr{O}$(10 9 ) DOFs using a new scalable parallel MPI/OpenMP implementation. Novel contributions of this paper include an alternative definition of coarse-grid systems based on restriction of fine-grid operators, yielding superior convergence results. In the uniform refinement setting, a detailed convergence study is provided, demonstrating h and p robust convergence and linear scaling with respect to the wave frequency. Finally, the paper concludes with numerical results on hp -adaptive simulations including a large-scale seismic modeling benchmark problem with high material contrast.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Multimodal Bayesian registration of noisy functions using Hamiltonian Monte Carlo

Functional data registration is a necessary processing step for many applications. The observed data can be inherently noisy, often due to measurement error or natural process uncertainty; which most functional alignment methods cannot handle. A pair of functions can also have multiple optimal alignment solutions, which is not addressed in current literature. In this paper, a flexible Bayesian approach to functional alignment is presented, which appropriately accounts for noise in the data without any pre-smoothing required. Additionally, by running parallel MCMC chains, the method can account for multiple optimal alignments via the multi-modal posterior distribution of the warping functions. To most efficiently sample the warping functions, the approach relies on a modification of the standard Hamiltonian Monte Carlo to be well-defined on the infinite-dimensional Hilbert space. In this work, this flexible Bayesian alignment method is applied to both simulated data and real data sets to show its efficiency in handling noisy functions and successfully accounting for multiple optimal alignments in the posterior; characterizing the uncertainty surrounding the warping functions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Broadband unidirectional visible imaging using wafer-scale nano-fabrication of multi-layer diffractive optical processors

We present a broadband and polarization-insensitive unidirectional imager that operates at the visible part of the spectrum, where image formation occurs in one direction, while in the opposite direction, it is blocked. This approach is enabled by deep learning-driven diffractive optical design with wafer-scale nano-fabrication using high-purity fused silica to ensure optical transparency and thermal stability. Our design achieves unidirectional imaging across three visible wavelengths (covering red, green, and blue parts of the spectrum), and we experimentally validated this broadband unidirectional imager by creating high-fidelity images in the forward direction and generating weak, distorted output patterns in the backward direction, in alignment with our numerical simulations. This work demonstrates wafer-scale production of diffractive optical processors, featuring 16 levels of nanoscale phase features distributed across two axially aligned diffractive layers for visible unidirectional imaging. This approach facilitates mass-scale production of ~0.5 billion nanoscale phase features per wafer, supporting high-throughput manufacturing of hundreds to thousands of multi-layer diffractive processors suitable for large apertures and parallel processing of multiple tasks. Beyond broadband unidirectional imaging in the visible spectrum, this study establishes a pathway for artificial-intelligence-enabled diffractive optics with versatile applications, signaling a new era in optical device functionality with industrial-level, massively scalable fabrication.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phased plan for the implementation of the time-resolving magnetic recoil spectrometer on the National Ignition Facility (NIF)

The time-resolving magnetic recoil spectrometer (MRSt) is a transformative diagnostic that will be used to measure the time-resolved neutron spectrum from an inertial confinement fusion implosion at the National Ignition Facility (NIF). It uses a CD foil on the outside of the hohlraum to convert fusion neutrons to recoil deuterons. An ion-optical system positioned outside the NIF target chamber energy-disperses and focuses forward-scattered deuterons. A pulse-dilation drift tube (PDDT) subsequently dilates, un-skews, and detects the signal. While the foil and ion-optical system have been designed, the PDDT requires more development before it can be implemented. Therefore, a phased plan is presented that first uses the foil and ion-optical systems with detectors that can be implemented immediately—namely CR-39 and hDISC streak cameras. These detectors will allow the MRSt to be commissioned in an intermediate stage and begin collecting data on a reduced timescale, while the PDDT is developed in parallel. A CR-39 detector will be used in phase 1 for the measurement of the time-integrated neutron spectra with excellent energy-resolution, necessary for the energy calibration of the system. Streak cameras will be used in phase 2 for measurement of the time-resolved spectrum with limited spectral coverage, which is sufficient to diagnose the time-resolved ion temperature. Simulations are presented that predict the performance of the streak camera detector, indicating that it will achieve excellent burn history measurements at current yields, and good time-resolved ion-temperature measurements at yields above 3 × 10 17 . The PDDT will be used for optimal efficiency and resolution in phase 3.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Impact of collisionality on turbulence in the edge of tokamak plasma using 3D global simulations

Collisionality is one of the key parameters in determining turbulent transport in the plasma edge, regulating phenomena such as 'shoulder formation', separation of scale lengths in the scrape-off layer (SOL), turbulence damping and zonal flow dynamics. Understanding its role is therefore of primary importance for future reactors like ITER. Getting reliable predictions and a better characterization of plasma flow properties when varying collisionality remains, however, a critical challenge for the simulations. This paper focuses on the impact of varying collisionality in a non-isothermal three-dimensional fluid model of the plasma edge. A high field side limited configuration encompassing open and closed magnetic field lines with parameters typical of a medium-sized tokamak is considered. The present model can consistently account for the variations of collisionality and its impact on both the parallel resistivity η ∥ and the ion and electron parallel thermal conductivities χ ∥e,i . Details on mean flow and turbulence properties are given. Changing collisionality leads to significant changes in the flow properties both on the mean and fluctuating quantities. In particular, lowering collisionality decreases the size of coherent structures, the fluctuation levels of turbulence, and steepens the density and temperature equilibrium profiles around the separatrix leading to a global reduction of the turbulent transport. The SOL width is observed to increase with collisionality, eventually resulting in the disappearance of the scale lengths separation between near and far SOL, consistently with previous experimental observations. At low collisionality, where the presence of narrow feature is well-established, a contribution of heat conduction increases up to compete with heat convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Cryogenic light detectors with thermal signal amplification for 0 νββ search experiments

As a step towards the realization of cryogenic-detector experiments to search for neutrinoless double-beta decay (such as CROSS, BINGO, and CUPID), we investigated a batch of 10 Ge light detectors (LDs) assisted by Neganov-Trofimov-Luke (NTL) signal amplification. Each LD was assembled with a large cubic light-emitting crystal (45 mm side) using the recently developed CROSS mechanical structure. The detector array was operated at milli-Kelvin temperatures in a pulse-tube cryostat at the Canfranc underground laboratory in Spain. We achieved good performance with scintillating bolometers from CROSS, made of Li 2 100 MoO 4 crystals and used as reference detectors of the setup, and with all LDs tested (except for a single device that encountered an electronics issue). No leakage current was observed for 8 LDs with an electrode bias up to 100 V. Operating the LDs at an 80 V electrode bias applied in parallel, we obtained a gain of around 9 in the signal-to-noise ratio of these devices, allowing us to achieve a baseline noise RMS of O(10 eV). Thanks to the strong current polarization of the temperature sensors, the time response of the devices was reduced to around half a millisecond in rise time. The achieved performance of the LDs was extrapolated via simulations of pile-up rejection capability for several configurations of the CUPID detector structure. Despite the sub-optimal noise conditions of the LDs (particularly at high frequencies), we demonstrated that the NTL technology provides a viable solution for background reduction in CUPID.

47 OTHER INSTRUMENTATION↗