Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equations of motion”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Modelling and simulation of brinicle formation

Below the Arctic sea ice, under the right conditions, a flux of icy brine flows down into the sea. The icy brine has a much lower fusion point and is denser than normal seawater. As a result, it sinks while freezing everything around it, forming an ice channel called a brinicle (also known as ice stalactite). In this paper, we develop a mathematical model for this phenomenon, assuming cylindrical symmetry. The fluid is considered to be viscous and quasi-stationary. The heat and salt transport are weakly coupled to the fluid motion and are modelled with the corresponding conservation equations, accounting for diffusive and convective effects. Finite-element discretization is employed to solve the coupled system of partial differential equations. We find that the model can capture the general behaviour of the physical system and generate brinicle-like structures while also recovering dendrite composition, which is a physically expected feature aligned with previous experimental results. This represents, to our knowledge, the first complete model proposed that captures the global structure of the physical phenomenon even though it has some discrepancies, such as brine accumulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Revisiting Theoretical Limits for One-Degree-of-Freedom Wave Energy Converters

This work revisits the theoretical limits of one-degree-of-freedom wave energy converters. This paper considers the floating sphere used in the Ocean Energy Systems Task 10 Wave Energy Converter modeling and verification effort for analysis. Analytical equations are derived to determine bounds on the motion amplitude, time-averaged power, and power-take-off (PTO) force. A unique result was found that shows the time-averaged power absorbed by a wave energy converter can be defined solely by the inertial properties and radiation hydrodynamic coefficients. In addition, a unique expression for the PTO force amplitude was derived that has provided upper and lower bounds when resistive control is used to maximize power generation. For complex conjugate control, this same expression can only provide a lower bound, as there is theoretically no upper bound. These bounds are used to compare the performance of a floating sphere if it were to extract energy using surge or heave motion. The analysis shows that because of the differences in hydrodynamic coefficients of each oscillating mode, there will be different frequency ranges that provide better power capture efficiency. The influence of a motion constraint on power absorption while also utilizing a nonideal power take-off is examined and found to reduce the losses associated with bidirectional energy flow. The expression to calculate the time-averaged power with a nonideal PTO is modified by the mechanical-to-electrical efficiency and the ratio of the PTO spring and damping coefficients. The PTO spring and damping coefficients were separated in the expression, which allows for limits to be set on the possible values of PTO coefficients to ensure a net flow of power to the grid.

50 EE - Wind and Water Power Program - Water (EE-4↗
DOE OSTI2025

Real-Time Bayesian Inference at Extreme Scale: A Digital Twin for Tsunami Early Warning Applied to the Cascadia Subduction Zone

We present a Bayesian inversion-based digital twin that employs acoustic pressure data from seafloor sensors, along with 3D coupled acoustic–gravity wave equations, to infer earthquake-induced spatiotemporal seafloor motion in real time and forecast tsunami propagation toward coastlines for early warning with quantified uncertainties. Our target is the Cascadia subduction zone, with one billion parameters. Computing the posterior mean alone would require 50 years on a 512 GPU machine. Instead, exploiting the shift invariance of the parameter-to-observable map and devising novel parallel algorithms, we induce a fast offline–online decomposition. The offline component requires just one adjoint wave propagation per sensor; using MFEM, we scale this part of the computation to the full El Capitan system (43,520 GPUs) with 92% weak parallel efficiency. Moreover, given real-time data, the online component exactly solves the Bayesian inverse and forecasting problems in 0.2 seconds on a modest GPU system, a ten-billion-fold speedup.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Effective Strength of Sinusoidally-Varying Focussing

Particle beam focussing and other confined physical systems often obey the equation ¨x = kx, where k < 0 implies stable motion. A constant negative value of k yields simple harmonic motion (constant focussing), while positive values give unstable (defocussing) motion. Alternating gradient focussing is the surprising principle that if k(t) is varied periodically, the overall motion can sometimes be stable even if the average value of k is zero or negative. It is used where focussing in one direction implies defocussing in one or more other directions, such as with electromagnetic fields in free space that have to satisfy ∇ · E = 0 and similar.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Revisiting Theoretical Limits for One-Degree-of-Freedom Wave Energy Converters: Preprint

This work revisits the theoretical limits of one-degree of freedom wave energy converters (WECs). This paper considers the floating sphere used in the Ocean Energy Systems (OES) Task 10 Wave Energy Converter modeling and verification effort for analysis. Analytical equations are derived to determine bounds on the motion amplitude, time-averaged power, and power take-off (PTO) force. A unique result was found that shows the time-averaged power absorbed by a WEC can be defined solely by the inertial properties and the radiation hydrodynamic coefficients. In addition, a unique expression for the PTO force amplitude was derived which has provided upper and lower bounds when resistive control is used to maximize power generation. For complex conjugate control, this same expression can only provide a lower bound as there is theoretically no upper bound. These bounds are used to compare the performance of a floating sphere if it were to extract energy using surge or heave motion. The analysis shows that because of the differences in hydrodynamic coefficients of each oscillating mode there will be different frequency ranges that provide better power capture efficiency. The influence of a motion constraint on power absorption while also utilizing a nonideal power take-off is examined and found to reduce the losses associated with bidirectional energy flow. The expression to calculate the time-averaged power with a nonideal PTO is modified by the mechanical-to-electrical efficiency and the ratio of the PTO spring and damping coefficients. These were separated in the expression which allows for limits to be set on the possible values of PTO coefficients to ensure a net flow of power to the grid.

50 EE - Wind and Water Power Program - Water (EE-4↗
DOE OSTI2022

Effect of Solvent Motion on Ion Transport in Electrolytes

We use concentrated solution theory to derive an equation governing solvent velocity in a binary electrolyte when a current passes through it. This equation, in combination with the material balance equation, enables the prediction of electrolyte concentration profiles and species velocities as a function of space and time. This framework is used to predict ion velocities in Li-Li symmetric cells containing a mixture of lithium bis(trifluoromethanesulfonyl)imide and poly(ethylene oxide) (LiTFSI/PEO), for which the cation transference number relative to the solvent velocity, ${t}_{+}^{0},$ can be either positive or negative, depending on salt concentration. Accounting for the solvent motion is increasingly important at higher concentrations. Especially for negative ${t}_{+}^{0},$ if solvent velocity is set to zero, the cation velocity, based on the electrode-electrolyte interface reference frame, is pointed opposite to the current flow. However, when solvent motion is taken into account, the cation velocity, based on the same reference frame, is in the same direction as the current. This analysis demonstrates the importance of accounting for solvent velocity rigorously in seemingly simple systems such as symmetric cells.

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

On the relationship between eastern China aerosols and western North Pacific tropical cyclone activity

How aerosols affect the weather and climate system has received increasing attention. This study finds a significant negative correlation between March–May eastern China aerosol optical depth (AOD) and July–November western North Pacific (WNP) tropical cyclone (TC) frequency during 2003–2020. This time period spans when several aerosol reanalyses are available and both Terra and Aqua Moderate Resolution Imaging Spectroradiometer AOD retrievals are assimilated therein. Composite analyses and budget analyses of dynamical genesis potential indices indicate the importance of large-scale environmental factors, especially vertical velocity and vertical wind shear, associated with changes in AOD that in turn modulate changes in WNP TC frequency. Increased eastern China AOD may facilitate negative Pacific meridional mode development via modulation of the westerly jet, which then forces an anticyclonic circulation over the WNP basin. Increased AOD can also directly decrease the inter-hemispheric temperature differential and increases the intra-hemispheric temperature differential between the equator and the mid-latitudes, thus weakening ascending motion and enhancing vertical wind shear over the WNP, especially the southeastern portion of the basin. Further, all of these large-scale environment changes induced by increased eastern China AOD tend to suppress WNP TCs. This study highlights the potential influence of eastern China aerosol loadings on WNP TCs, thus improving our understanding of TC climate variability over the WNP.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

EZ: An efficient, charge conserving current deposition algorithm for electromagnetic particle-in-cell simulations

Here, we present EZ, a novel current deposition algorithm for particle-in-cell (PIC) simulations. EZ calculates the current density on the electromagnetic grid due to macro-particle motion within a time step by solving the continuity equation of electrodynamics. Being a charge conserving hybridization of Esirkepov's method and ZigZag, we refer to it as “EZ” as shorthand for “Esirkepov meets ZigZag”. Simulations of a warm, relativistic plasma with PIConGPU show that EZ achieves the same level of charge conservation as the commonly used method by Esirkepov, yet reaches higher performance for macro-particle assignment-functions up to third-order. In addition to a detailed description of the functioning of EZ, reasons for the expected and observed performance increase are given, and guidelines for its implementation aiming at highest performance on GPUs are provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Goal-oriented real-time Bayesian inference for linear autonomous dynamical systems with application to digital twins for tsunami early warning

We present a goal-oriented framework for constructing digital twins with the following properties: (1) they employ discretizations of high-fidelity partial differential equation (PDE) models governed by autonomous dynamical systems, leading to large-scale forward problems; (2) they solve a linear inverse problem to assimilate observational data to infer uncertain model components followed by a forward prediction of the evolving dynamics; and (3) the entire end-to-end, data-to-inference-to-prediction computation is carried out without approximation and in real time through a Bayesian framework that rigorously accounts for uncertainties. Several challenges must be overcome to realize this framework, including the large scale of the forward problem, the high dimensionality of the parameter space, and for a class of problems including those we target, the slow decay of the singular values of the parameter-to-observable map. Here we introduce a methodology to overcome these challenges by exploiting the autonomous structure of the forward model to decompose the solution of the inverse problem into a one-time-only offline phase in which the PDE model is solved a limited number of times (equal to the number of sensors), and an online phase that maps well onto GPUs and computes the parameter inference and prediction of quantities of interest in real time, given observational data. Our ultimate goal is to apply this framework to construct digital twins for subduction zones, including Cascadia, to provide early warning for tsunamis generated by megathrust earthquakes. To this end, we demonstrate how our methodology can be used to employ seafloor pressure observations, along with the coupled acoustic–gravity wave equations, to infer the earthquake-induced spatiotemporal seafloor motion (discretized with $\mathscr{O}$ (10 9 ) parameters) and forward predict the tsunami propagation. We present results of an end-to-end inference, prediction, and uncertainty quantification for a representative test problem with $\mathscr{O}$ (10 8 ) inversion parameters for which goal-oriented Bayesian inference is accomplished exactly and in real time, that is, in a matter of seconds.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

The Pantheon+ Analysis: Evaluating Peculiar Velocity Corrections in Cosmological Analyses with Nearby Type Ia Supernovae

Separating the components of redshift due to expansion and peculiar motion in the nearby universe (z < 0.1) is critical for using Type Ia Supernovae (SNe Ia) to measure the Hubble constant (H 0 ) and the equation-of-state parameter of dark energy (w). Here, we study the two dominant "motions" contributing to nearby peculiar velocities: large-scale, coherent-flow (CF) motions and small-scale motions due to gravitationally associated galaxies deemed to be in a galaxy group. We use a set of 584 low-z SNe from the Pantheon+ sample, and evaluate the efficacy of corrections to these motions by measuring the improvement of SN distance residuals. We study multiple methods for modeling the large and small-scale motions and show that, while group assignments and CF corrections individually contribute to small improvements in Hubble residual scatter, the greatest improvement comes from the combination of the two (relative standard deviation of the Hubble residuals, Rel. SD, improves from 0.167 to 0.157 mag). We find the optimal flow corrections derived from various local density maps significantly reduce Hubble residuals while raising H 0 by ~0.4 km s -1 Mpc -1 as compared to using CMB redshifts, disfavoring the hypothesis that unrecognized local structure could resolve the Hubble tension. We estimate that the systematic uncertainties in cosmological parameters after optimally correcting redshifts are 0.06–0.11 km s -1 Mpc -1 in H 0 and 0.02–0.03 in w which are smaller than the statistical uncertainties for these measurements: 1.5 km s -1 Mpc -1 for H 0 and 0.04 for w.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Heat transfer enhancement due to cold cap motion from bubbling in a waste glass melter

In this study, a computational fluid dynamics (CFD) model was developed to model the motion of a solid cold cap in a waste glass melter. Forced convection bubblers at the base of the melter release air into the molten glass, which forms large bubbles that travel upwards to the cold cap and augment heat transfer from the glass to the cold cap. The CFD model employs the Navier- Stokes equations to solve for the fluctuating flowfield using a rigid body motion dynamic fluid body interaction module. This allows for movement of the floating body in response to the bubbling forces calculated at each time step. Further, the heat flux delivered to the cold cap by the convective bubbling is studied as a function of the bubbling rate. Results for the moving cold cap are compared with the computed heat flux trends for a stationary cold cap. The heat flux delivered to the cold cap from the glass is 25% higher for the case with the moving cold cap. The heat flux was found to be proportional to v0.6 as opposed to v0.9 (where v is the normalized bubbling rate) for the stationary cold cap.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Multi-fidelity modeling to predict the rheological properties of a suspension of fibers using neural networks and Gaussian processes

Unveiling the rheological properties of fiber suspensions is of paramount interest to many industrial applications. There are multiple factors, such as fiber aspect ratio and volume fraction, that play a significant role in altering the rheological behavior of suspensions. Three-dimensional (3D) numerical simulations of coupled differential equations of the suspension of fibers are computationally expensive and time-consuming. Machine learning algorithms can be trained on the available data and make predictions for the cases where no numerical data are available. However, some widely used machine learning surrogates, such as neural networks, require a relatively large training dataset to produce accurate predictions. Multi-fidelity models, which combine high-fidelity data from numerical simulations and less expensive lower fidelity data from resources such as simplified constitutive equations, can pave the way for more accurate predictions. Here, we focus on neural networks and the Gaussian processes with two levels of fidelity, i.e., high and low fidelity networks, to predict the steady-state rheological properties, and compare them to the single-fidelity network. High-fidelity data are obtained from direct numerical simulations based on an immersed boundary method to couple the fluid and solid motion. The low-fidelity data are produced by using constitutive equations. Multiple neural networks and the Gaussian process structures are used for the hyperparameter tuning purpose. Results indicate that with the best choice of hyperparameters, both the multi-fidelity Gaussian processes and neural networks are capable of making predictions with a high level of accuracy with neural networks demonstrating marginally better performance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Milky Way mass with K giants and BHB stars using LAMOST, SDSS/SEGUE, and Gaia : 3D spherical Jeans equation and tracer mass estimator

ABSTRACT We measure the enclosed Milky Way mass profile to Galactocentric distances of ∼70 and ∼50 kpc using the smooth, diffuse stellar halo samples of Bird et al. The samples are Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) and Sloan Digital Sky Survey/Sloan Extension for Galactic Understanding and Exploration (SDSS/SEGUE) K giants (KG) and SDSS/SEGUE blue horizontal branch (BHB) stars with accurate metallicities. The 3D kinematics are available through LAMOST and SDSS/SEGUE distances and radial velocities and Gaia DR2 proper motions. Two methods are used to estimate the enclosed mass: 3D spherical Jeans equation and Evans et al. tracer mass estimator (TME). We remove substructure via the Xue et al. method based on integrals of motion. We evaluate the uncertainties on our estimates due to random sampling noise, systematic distance errors, the adopted density profile, and non-virialization and non-spherical effects of the halo. The tracer density profile remains a limiting systematic in our mass estimates, although within these limits we find reasonable agreement across the different samples and the methods applied. Out to ∼70 and ∼50 kpc, the Jeans method yields total enclosed masses of 4.3 ± 0.95 (random) ±0.6 (systematic) × 1011 M⊙ and 4.1 ± 1.2 (random) ±0.6 (systematic) × 1011 M⊙ for the KG and BHB stars, respectively. For the KG and BHB samples, we find a dark matter virial mass of $M_{200}=0.55^{+0.15}_{-0.11}$ (random) ±0.083 (systematic) × 1012 M⊙ and $M_{200}=1.00^{+0.67}_{-0.33}$ (random) ±0.15 (systematic) × 1012 M⊙, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Progress towards a phenomenological picture and theoretical understanding of glassy dynamics and vitrification near interfaces and under nanoconfinement

The nature of alterations to dynamics and vitrification in the nanoscale vicinity of interfaces-commonly referred to as "nanoconfinement" effects on the glass transition-has been an open question for a quarter century. We first analyze experimental and simulation results over the last decade to construct an overall phenomenological picture. Key features include the following: after a metrology- and chemistry-dependent onset, near-interface relaxation times obey a fractional power law decoupling relation with bulk relaxation; relaxation times vary in a double-exponential manner with distance from the interface, with an intrinsic dynamical length scale appearing to saturate at low temperatures; the activation barrier and vitrification temperature T-g approach bulk behavior in a spatially exponential manner; and all these behaviors depend quantitatively on the nature of the interface. We demonstrate that the thickness dependence of film-averaged T-g for individual systems provides a poor basis for discrimination between different theories, and thus we assess their merits based on the above dynamical gradient properties. Entropy-based theories appear to exhibit significant inconsistencies with the phenomenology. Diverse free-volume-motivated theories vary in their agreement with observations, with approaches invoking cooperative motion exhibiting the most promise. The elastically cooperative nonlinear Langevin equation theory appears to capture the largest portion of the phenomenology, although important aspects remain to be addressed. A full theoretical understanding requires improved confrontation with simulations and experiments that probe spatially heterogeneous dynamics within the accessible 1-ps to 1-year time window, minimal use of adjustable parameters, and recognition of the rich quantitative dependence on chemistry and interface. Published under license by AIP Publishing.

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

Theory, analysis, and testing of an angular resonator for wave energy generation

This article describes the theory, analysis, and initial bench-top testing of a minimally invasive, rotational resonator designed to produce small amounts of electrical energy for use in oceanic observation buoys. Here, this work details the systems of equations that govern such a resonator, its potential power production, and its predicted effects on the modified motion of the buoy. Finally, a bench-top test apparatus is designed and experimented upon to identify the system and verify the system of equations empirically.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Collective mode Brownian dynamics: A method for fast relaxation of statistical ensembles

Sampling equilibrium configurations of correlated systems of particles with long relaxation times (e.g., polymeric solutions) using conventional molecular dynamics and Monte Carlo methods can be challenging. This is especially true for systems with complicated, extended bond network topologies and other interactions that make the use and design of specialized relaxation protocols infeasible. We introduce a method based on Brownian dynamics simulations that can reduce the computational time it takes to reach equilibrium and draw decorrelated samples. Importantly, the method is completely agnostic to the particle configuration and the specifics of interparticle forces. In particular, we develop a mobility matrix that excites non-local, collective motion of N particles and can be computed efficiently in $\mathcal{O}(N)$ time. Particle motion in this scheme is computed by integrating the overdamped Langevin equation with an Euler–Maruyama scheme, in which Brownian displacements are drawn efficiently using a low-rank representation of the mobility matrix in position and wave space. We demonstrate the efficacy of the method with various examples from the realm of soft condensed matter and release a massively parallel implementation of the code as a plugin for the open-source package HOOMD-blue [J. A. Anderson et al., J. Comput. Phys. 227, 5342 (2008) and J. Glaser et al., Comput. Phys. Commun. 192, 97 (2015)] which runs on graphics processing units

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

T-BMT Spin Resonance Tracker Code for He3 with Six Snakes

Polarization lifetime for He3 using two and six snakes are studied using the T-BMT Spin Resonance Tracker code. This code integrates a reduced spinor form of the T-BMT equation including only several spin resonances and the kinematics of synchrotron motion. It was previously benchmarked against RHIC polarization lifetime under the two snake system *. * Phys. Rev.Accel. Beams 22 (2019) 9, 091001

Ranjbar, V. H.↗
DOE OSTI2022

A classical equation that accounts for observations of non-Arrhenius and cryogenic grain boundary migration

Observations of microstructural coarsening at cryogenic temperatures, as well as numerous simulations of grain boundary motion that show faster migration at low temperature than at high temperature, have been troubling because they do not follow the expected Arrhenius behavior. This work demonstrates that classical equations, that are not simplified, account for all these oddities and demonstrate that non-Arrhenius behavior can emerge from thermally activated processes. According to this classical model, this occurs when the intrinsic barrier energies of the processes become small, allowing activation at cryogenic temperatures. Additional thermal energy then allows the low energy process to proceed in reverse, so increasing temperature only serves to frustrate the forward motion. This classical form is shown to reconcile and describe a variety of diverse grain boundary migration observations.

36 MATERIALS SCIENCE↗