Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

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

At least 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Sensitivity of time-dependent density functional theory to initial conditions

Time-dependent density-functional theory is mathematically formulated through nonlinear coupled time-dependent three-dimensional partial differential equations, and it is natural to expect a strong sensitivity of its solutions to variations of the initial conditions, akin to the butterfly effect ubiquitous in classical dynamics. Since the Schrödinger equation for an interacting many-body system is, however, linear and mathematically the exact equations of the density-functional theory reproduce the corresponding one-body properties, it would follow that the Lyapunov exponents are also vanishing within a density-functional theory framework. Whether for realistic implementations of the time-dependent density-functional theory the question of the absence of the butterfly effect and whether the dynamics provided is indeed a predictable theory was never discussed. At the same time, since the time-dependent density-functional theory is a unique tool allowing us to study the nonequilibrium dynamics of strongly interacting many-fermion systems, the question of predictability of this theoretical framework is of paramount importance. Here our analysis, for a number of quantum superfluid many-body systems (unitary Fermi gas, nuclear fission, and heavy-ion collisions) with a classical equivalent number of degrees of freedom O(10 10 ) and larger, suggests that its maximum Lyapunov exponents are negligible for all practical purposes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Transient Efficiency Flexibility and Reliability Optimization of Coal-Fired Power Plants: Model-Predictive Control Library Development (Report)

This document pertains to the reporting requirements of DOE contract FE-0031767. The document covers the development of a model predictive control (MPC) library for implementing MPC for a general dynamic system. The library is implemented in a standardized manner in Matlab/Simulink, where core functions on model prediction, linearization and formulation and solution of a quadratic programming (QP) optimization problem is done in the core library - independent of the specific application. The user can provide the application-specific dynamic model in continuous and discrete time, to rapidly implement and test the MPC performance in a desktop simulation. The MPC optimization objective and constraints are also easily configured via an Excel file to allow iterative refinement as needed. Finally, the MPC library enables a rapid deployment to a target environment through auto C-code generation and containerization. The MPC library works seamlessly with the model based estimation (MBE) library to obtain the overall output feedback control solution. In this program, the reduced order model (ROM) of a coal-fired power plant (CFPP) is used to implement and test the MPC solution.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

Effect of Sulfonation Level on the Percolated Morphology and Proton Conductivity of Hydrated Fluorine-Free Copolymers: Experiments and Simulations

Using all-atom molecular dynamics simulations and a variety of experimental methods, we previously reported on a linear polyethylene with pendant phenyl sulfonated groups precisely on every fifth carbon along the backbone. With increasing relative humidity this fluorine-free polymer self-assembled to form nanoscale water channels and exhibited exceptional proton conductivity. Expanding upon those findings, here we explore partially sulfonated random copolymers, referred to as p 5PhSH-Y. Using either acetyl sulfate or sulfuric acid, a wide range of sulfonation levels were prepared ( Y = 34−98%) corresponding to ion-exchange capacities (IEC) of 2.0−4.4 mmol/g. Combining experimental techniques and all-atom molecular dynamics simulations, we study the effect of Y on water uptake, nanoscale morphology, and the proton/water transport properties of p5PhSH- Y . The proton conductivity of p 5PhSH- Y increases with relative humidity and with Y and achieves values in excess of 0.1 S/cm. These high conductivities are attributed to high IEC and welldeveloped nanoscale percolated hydrophilic domains made possible by the flexible backbone. We quantitatively describe the nature of the water channels using the characteristic distance, channel width distribution, the area per sulfonate group at the hydrophilic/ hydrophobic interface, and the fractal dimension. Notably, the channel widths and the areas per sulfonate group are nominally independent of the level of sulfonation, while depending significantly on the level of hydration. The fractal dimension of the water channels correlates strongly with the water diffusion coefficients calculated from the molecular dynamics (MD) simulations. These findings demonstrate that the p 5PhSH- Y hydrocarbon copolymers can be modified to tune properties, particularly proton conductivity.

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

Defect-Limited Mobility and Defect Thermochemistry in Mixed A-Cation Tin Perovskites: (CH3NH3)1-xCsxSnBr3

Hybrid organic-inorganic semiconductors crystallizing in the perovskite structure present a significant opportunity for realizing defect-tolerant semiconductors. In this work, we examine the solid solution, (CH3NH3)1-xCsxSnBr3, and identify the thermochemistry dictating the intrinsic carrier concentrations and how local structural distortions influences this electronic behavior. This family of compounds exhibits the expected systematic trend in decreasing optical gap with the cesium to methylammonium A-site mixing ratio, x, in the visible region. However, the carrier mobility, as determined from time-resolved microwave conductivity measurements, trends opposite to that expected from first-principles calculations combined with Boltzmann scattering theory calculations of the carrier mobility. We propose that this is a result of increasing carrier scattering with x in (CH3NH3)1-xCsxSnBr3 due to a significant increase in the carrier density with x. By examining the dependence of the carrier density as a function of x, we infer the compositional-dependence of the average enthalpy and nonconfigurational entropy per defect. While diffraction reveals a cubic aristotypic perovskite structure as a function of x at room temperature, the pair distribution functions obtained from synchrotron X-ray total scattering from these materials are better described by symmetry-adapted displacement modes of the Pm3m crystal structure, which we attribute to a large degree of anharmonic dynamics of the atom positions. This analysis shows that CH3NH3+-rich compositions retain mostly linear Sn-Br-Sn bonding environments through the displacements, while Cs+-rich compositions lead to more significantly bent Sn-Br-Sn environments. We propose that these bent bonding arrangements in Cs-rich compositions yield a higher propensity for defect formation. This also provides a rationale for carrier trapping that gives rise to anomalous microwave transients. Together, these results provide insight into the structure-dynamics-properties relationships in this highly anharmonic system with high amplitude atomic motions and low defect formation energies.

defect thermochemistry↗
DOE OSTI2022

Dynamics of pectic homogalacturonan in cellular morphogenesis and adhesion, wall integrity sensing and plant development

Homogalacturonan (HG) is the most abundant pectin subtype in plant cell walls. Although it is a linear homopolymer, its modification states allow for complex molecular encoding. HG metabolism affects its structure, chemical properties, mobility and binding capacity, allowing it to interact dynamically with other polymers during wall assembly and remodelling and to facilitate anisotropic cell growth, cell adhesion and separation, and organ morphogenesis. HGs have also recently been found to function as signalling molecules that transmit information about wall integrity to the cell. Here we highlight recent advances in our understanding of the dual functions of HG as a dynamic structural component of the cell wall and an initiator of intrinsic and environmental signalling. Here, we also predict how HG might interconnect the cell wall, plasma membrane and intracellular components with transcriptional networks to regulate plant growth and development.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Extended x-ray emission times of clusters in intense x-ray pulses

Novel experimental spectroscopic and imaging techniques have been developed to exploit x-ray free-electron laser (XFEL) for studying structure and dynamic of non-periodic entities at atomic spatial resolution and femtosecond timescales.Recently, intensity correlation of x-ray fluorescence, based on the principle introduced by Hanbury Brown and Twiss [R. Hanbury Brown and R. Q. Twiss, Nature (London) 177, 27 (1956)], has been proposed for high-resolution imaging of a three-dimensional arrangement of atoms. To explore the applicability of this fluorescence approach, we theoretically investigate fluorescence dynamics of nonperiodic systems subject to x-ray free electron laser (XFEL) pulses over a range of fluences from the linear to nonlinear x-ray absorption regimes as a function of system size, from a single atom to a cluster of 149 171 atoms. Fluorescence dynamics in the nonlinear x-ray regime differs from that in the weak x-ray field in that intense x-ray pulses interrupt the fluorescence dynamics by multiphoton absorption creating a dense electron environment within the sample on a femtosecond timescale. In large systems, the presence of both recombination and photoionization pathways gives rise to an enhanced Kα and Kα h emission yield and an extended emission time beyond the lifetime of the core-excited states. Our analysis suggests that, in an intense x-ray pulse, the Kα h emission line can be a good candidate for fluorescence imaging as it has a higher yield and an emission time that is short relative to the x-ray-induced lattice distortion time.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

QLBT: a linear Boltzmann transport model for heavy quarks in a quark-gluon plasma of quasi-particles

Abstract We develop a new heavy quark transport model, QLBT, to simulate the dynamical propagation of heavy quarks inside the quark-gluon plasma (QGP) created in relativistic heavy-ion collisions. Our QLBT model is based on the linear Boltzmann transport (LBT) model with the ideal QGP replaced by a collection of quasi-particles to account for the non-perturbative interactions among quarks and gluons of the hot QGP. The thermal masses of quasi-particles are fitted to the equation of state from lattice QCD simulations using the Bayesian statistical analysis method. Combining QLBT with our advanced hybrid fragmentation-coalescence hadronization approach, we calculate the nuclear modification factor $$R_\mathrm {AA}$$ R AA and the elliptic flow $$v_2$$ v 2 of D mesons at the Relativistic Heavy-Ion Collider and the Large Hadron Collider. By comparing our QLBT calculation to the experimental data on the D meson $$R_\mathrm {AA}$$ R AA and $$v_2$$ v 2 , we extract the heavy quark transport parameter $$\hat{q}$$ q ^ and diffusion coefficient $$D_\mathrm {s}$$ D s in the temperature range of $$1-4~T_\mathrm {c}$$ 1 - 4 T c , and compare them with the lattice QCD results and other phenomenological studies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Numerical Analysis of Offshore Wind-Farm-Induced Drag Effects on Coastal Upwelling Dynamics

Wind farms extract momentum from the atmospheric flow, generating wind-speed deficits both within the plant, and extending downstream. When located offshore, these deficits modulate air-sea coupling, potentially impacting coastal upwelling in sensitive regions. We investigate impacts of wind farms on coastal upwelling using kilometre-scale, three-way-coupled simulations with the coupled ocean-atmosphere-wave-sediment transport system for the US West Coast. Wind-farm effects are represented by a generalised turbine drag formulation, an idealised, height-dependent body force whose magnitude is systematically varied. This approach isolates the leading-order fluid-dynamical response in a realistic coastal configuration. The atmospheric adjustment exhibits an approximately linear relation between drag force and wind-speed deficit, with wakes that expand downstream and increase in magnitude as drag increases. An empirical orthogonal function analysis of sea-surface-temperature anomalies reveals the emergence of a canonical dipole pattern under strong drag forcing. Subsurface diagnostics show consistent shoaling of the mixed layer and suppressed upward velocities in areas near wind-farm region, accompanied by compensating enhancements of shoaling closer to the coast. These results identify turbine drag as a control parameter in assessing interactions between wind-farm wake and coastal upwelling and provide scaling relationships for understanding offshore wind-farm effects on the coastal circulation dynamics.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Relaxation dynamics of deformed polymer nanocomposites as revealed by small-angle scattering and rheology

Here, the relaxation dynamics of polystyrene (PS)/silica nanocomposites after a large step deformation are studied by a combination of small-angle scattering techniques and rheology. Small-angle X-ray scattering measurements and rheology show clear signatures of nanoparticle aggregation that enhances the mechanical properties of the polymer nanocomposites (PNCs) in the linear viscoelastic regime and during the initial phase of stress relaxation along with accelerated relaxation dynamics. Small-angle neutron scattering experiments under the zero-average-contrast condition reveal, however, smaller structural anisotropy in the PNCs than that in the neat polymer matrix, as well as accelerated anisotropy relaxation. In addition, the degrees of anisotropy reduction and relaxation dynamics acceleration increase with increasing nanoparticle loading. These results are in sharp contrast to the prevailing viewpoint of enhanced molecular deformation as the main mechanism for the mechanical enhancement in PNCs. Furthermore, the observed acceleration of stress relaxation and reduction in structural anisotropy point to two types of nonlinear effects in the relaxation dynamics of PNCs at large deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

LaSDI: Parametric Latent Space Dynamics Identification

Enabling fast and accurate physical simulations with data has become an important area of computational physics to aid in inverse problems, design-optimization, uncertainty quantification, and other various decision-making applications. This paper presents a data-driven framework for parametric latent space dynamics identification procedure that enables fast and accurate simulations. The parametric model is achieved by building a set of local latent space model and designing an interaction among them. An individual local latent space dynamics model achieves accurate solution in a trust region. By letting the set of trust region to cover the whole parameter space, our model shows an increase in accuracy with an increase in training data. Herein we introduce two different types of interaction mechanisms, i.e., point-wise and region-based approach. Both linear and nonlinear data compression techniques are used. We illustrate the framework of Latent Space Dynamics Identification (LaSDI) enable a fast and accurate solution process on various partial differential equations, i.e., Burgers’ equations, radial advection problem, and nonlinear heat conduction problem, achieving 0 (100)x speed-up and 0 (1)% relative error with respect to the corresponding full order models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Threading–Unthreading Transition of Linear-Ring Polymer Blends in Extensional Flow

Adding small amounts of ring polymers to a matrix of their linear counterparts is known to increase the zero-shear-rate viscosity because of linear-ring threading. Uniaxial extensional rheology measurements show that, unlike its pure linear and ring constituents, the blend exhibits an overshoot in the stress growth coefficient. By combining these measurements with ex-situ small-angle neutron scattering and nonequilibrium molecular dynamics simulations, this overshoot is shown here to be driven by a transient threading–unthreading transition of rings embedded within the linear entanglement network. Prior to unthreading, embedded rings deform affinely with the linear entanglement network and produce a measurably stronger elongation of the linear chains in the blend compared to the pure linear melt. Thus, rings uniquely alter the mechanisms of transient elongation in linear polymers.

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

Modeling dynamics of acute HIV infection incorporating density-dependent cell death and multiplicity of infection

Understanding the dynamics of acute HIV infection can offer valuable insights into the early stages of viral behavior, potentially helping uncover various aspects of HIV pathogenesis. The standard viral dynamics model explains HIV viral dynamics during acute infection reasonably well. However, the model makes simplifying assumptions, neglecting some aspects of HIV infection. For instance, in the standard model, target cells are infected by a single HIV virion. Yet, cellular multiplicity of infection (MOI) may have considerable effects in pathogenesis and viral evolution. Further, when using the standard model, we take constant infected cell death rates, simplifying the dynamic immune responses. Here, we use four models—1) the standard viral dynamics model, 2) an alternate model incorporating cellular MOI, 3) a model assuming density-dependent death rate of infected cells and 4) a model combining (2) and (3)—to investigate acute infection dynamics in 43 people living with HIV very early after HIV exposure. We find that all models qualitatively describe the data, but none of the tested models is by itself the best to capture different kinds of heterogeneity. Instead, different models describe differing features of the dynamics more accurately. For example, while the standard viral dynamics model may be the most parsimonious across study participants by the corrected Akaike Information Criterion (AICc), we find that viral peaks are better explained by a model allowing for cellular MOI, using a linear regression analysis as analyzed by R 2 . These results suggest that heterogeneity in within-host viral dynamics cannot be captured by a single model. Depending on the specific aspect of interest, a corresponding model should be employed.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Penalized ensemble Kalman filters for high dimensional non-linear systems

The ensemble Kalman filter (EnKF) is a data assimilation technique that uses an ensemble of models, updated with data, to track the time evolution of a usually non-linear system. It does so by using an empirical approximation to the well-known Kalman filter. However, its performance can suffer when the ensemble size is smaller than the state space, as is often necessary for computationally burdensome models. This scenario means that the empirical estimate of the state covariance is not full rank and possibly quite noisy. To solve this problem in this high dimensional regime, we propose a computationally fast and easy to implement algorithm called the penalized ensemble Kalman filter (PEnKF). Under certain conditions, it can be theoretically proven that the PEnKF will be accurate (the estimation error will converge to zero) despite having fewer ensemble members than state dimensions. Further, as contrasted to localization methods, the proposed approach learns the covariance structure associated with the dynamical system. These theoretical results are supported with simulations of several non-linear and high dimensional systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Ultrafast electron dynamics of graphene quantum dots: High harmonic generation

We study theoretically nonlinear optical properties of graphene quantum dots placed in a field of a short and strong linearly polarized optical pulse. We address the problem of high harmonic generation in quantum dots and how such nonlinear effect is affected by dephasing processes in a quantum dot. The dephasing makes the ultrafast electron dynamics more irreversible with a large residual population of the excited quantum dot levels. In relation to the high-harmonic spectrum, with increasing the dephasing time, the intensities of the low-frequency harmonics increase while the cutoff energy decreases. The dependence of the cutoff energy on the amplitude of the optical pulse is also sensitive to the frequency of the pulse. When the frequency of the optical pulse is much less than the quantum dot band gap, this dependence is almost linear, but when the frequency of the pulse is comparable to the band gap, the cutoff energy shows saturation behavior at large field amplitude, >0.4 V/Å.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Tuning the charge transfer character of the multiexciton state in singlet fission

Intramolecular singlet fission (SF) produces the multiexciton correlated triplet pair state, (T 1 T 1 ), prior to the formation of free triplet excitons. The nature of the multiexciton state is complex, as generation of the (T 1 T 1 ) state may involve a charge transfer (CT) intermediate and has been shown to have both mixed electronic and spin characters. According to transient absorption spectroscopy, a linear terrylene-3,4:11,12-bis(dicarboximide) dimer (TDI 2 ) exhibits solvent-dependent excited-state dynamics. As solvent polarity increases from 1,2,4-trichlorobenzene (ε = 2.2) to chlorobenzene (ε = 5.6) to 1,2-dichlorobenzene (ε = 9.9), the SF rate in TDI2 increases and the multiexciton state, which can be thought of as a linear combination of the 1 (S 1 S 0 ), CT, and (T 1 T 1 ) states, gains more CT character. Eventually, the CT state becomes a trap state as indicated by symmetry-breaking charge separation in TDI2 in pyridine (ε = 12.3). The dielectric environment influences not only the SF rate and the relative contributions of the 1 (S 1 S 0 ), CT, and (T 1 T 1 ) states to the overall multiexciton state but also the rate at which the state mixing evolves, with faster dynamics in higher polarity solvents. More importantly, the tunability and presence of strong CT character in the multiexciton state have implications for SF applications since they often rely on electron transfer from the free triplet excitons. This enhanced CT character in the (T 1 T 1 ) state may assist with two-electron transfer directly from the (T 1 T 1 ) state, allowing for facile extraction of charges in intramolecular SF systems whose (T 1 T 1 ) states do not always efficiently dissociate to two triplet excitons.

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

Fine-root dynamics vary with soil depth and precipitation in a lown-utrient tropical forest in the Central Amazonia

A common assumption in tropical ecology is that root systems respond rapidly to climatic cues but that most of that response is limited to the uppermost layer of the soil, with relatively limited changes in deeper layers. However, this assumption has not been tested directly, preventing models from accurately predicting the response of tropical forests to environmental change. We measured seasonal dynamics of fine roots in a Central Amazonia mature forest using minirhizotrons to 90 cm depth, which were calibrated with fine roots extracted from soil cores. Root productivity and mortality in superficial soil layers were positively correlated with precipitation, whereas root standing length was greater during the dry periods at the deeper layers. Contrary to historical assumptions, a large fraction of fine-root standing biomass (46%) and productivity (41%) were found in soil layers deeper than 30 cm. Furthermore, root turnover decreased linearly with soil depth. Our findings demonstrate a strong relationship between fine-root dynamics and precipitation regimes in Central Amazonia. Our results also emphasize the importance of deep roots for accurate estimates of primary productivity and the interaction between roots and carbon, water, and nutrients.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Catastrophic Cooling Instability in Optically Thin Plasmas

The solar corona is the prototypical example of a low-density environment heated to high temperatures by external sources. The plasma cools radiatively, and because it is optically thin to this radiation, it becomes possible to model the density, velocity, and temperature structure of the system by modifying the MHD equations to include an energy source term that approximates the local heating and cooling rates. The solutions can be highly inhomogeneous and even multiphase because the well-known linear instability associated with this source term, thermal instability, leads to a catastrophic heating and cooling of the plasma in the nonlinear regime. Here we show that there is a separate, much simpler linear instability accompanying this source term that can rival thermal instability in dynamical importance. The stability criterion is the isochoric one identified by Parker (1953), and we demonstrate that cooling functions derived from collisional ionization equilibrium are highly prone to violating this criterion. If catastrophic cooling instability can act locally in global simulations, then it is an alternative mechanism for forming condensations, and due to its nonequilibrium character, it may be relevant to explaining a host of phenomena associated with the production of cooler gas in hot, low density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Inertial collapse of a gas bubble in a shear flow near a rigid wall

Despite the extensive research on bubble collapse near rigid walls, the bubble collapse dynamics in the presence of shear flow near a rigid wall is poorly understood. We conduct direct simulations of the Navier–Stokes equations to explore the bubble dynamics and pressures during bubble collapse near a rigid, flat wall under linear shear flow conditions. We examine the dependence of the bubble collapse morphology and wall pressures on the initial bubble location and shear rate. We find that shear distorts the bubble, generating two re-entrant jets – one developing from the side opposite to the mean flow and the other from the far end toward the wall. Upon impact of the jet on the opposite side of the bubble, water-hammer shocks are produced, which propagate outward and interact with the convoluted bubble shape. The shock stretches the bubble towards the wall, resulting in a closer impact location for the jet originating from the far end compared with the case with no shear flow. The water-hammer pressure location can be approximated as the theoretical distance travelled by a particle initialised at the bubble centre with the corresponding constant shear flow velocity. The maximum wall pressures can thus be predicted by considering the distance between the far jet impingement location and the wall along the wall-normal direction. As the shear rate is increased, the maximum wall pressure increases, although only marginally. We determine the critical initial stand-off distance from the wall at which the bubble morphology is shear dominated, i.e. characterised by converging re-entrant jets.

Bhola, Sahil (ORCID:0000000236766302)↗