Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamical System Scaling”

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 703 records · Page 39

A reduced-order model of concentration polarization in reverse osmosis systems with feed spacers

Feed spacers in reverse osmosis systems generate complex fluid flows that limit computational fluid dynamics (CFD) simulations to small length and time scales. That limits our ability to simulate mineral scaling and other membrane fouling phenomena, which occur over longer length and time scales. Thus motivated, we develop a reduced model that replaces the CFD simulation of the velocity field with an analytical model that mimics spacers. This focuses the remaining numerical effort on simulating the advection–diffusion equation governing solute transport. We motivate and validate the model with CFD simulations and bench-scale experiments of spacer filaments in three different arrangements, including cases of unsteady vortex shedding. We show that the model produces a roughly 10,000-fold speedup compared to CFD, and accurately reproduces CFD predictions of not only the average and maximum concentrations, but also the local concentration distribution along the membrane. Here, we also demonstrate the model for simulating a feed channel with a length-to-height ratio of 200. The model provides a simple testbed for exploratory studies of multispecies transport, precipitation, and membrane fouling phenomena for which simulating spacers is often prohibitive.

Computational fluid dynamics↗

Phonon-enhanced nonlinearities in hexagonal boron nitride

Abstract Polar crystals can be driven into collective oscillations by optical fields tuned to precise resonance frequencies. As the amplitude of the excited phonon modes increases, novel processes scaling non-linearly with the applied fields begin to contribute to the dynamics of the atomic system. Here we show two such optical nonlinearities that are induced and enhanced by the strong phonon resonance in the van der Waals crystal hexagonal boron nitride (hBN). We predict and observe large sub-picosecond duration signals due to four-wave mixing (FWM) during resonant excitation. The resulting FWM signal allows for time-resolved observation of the crystal motion. In addition, we observe enhancements of third-harmonic generation with resonant pumping at the hBN transverse optical phonon. Phonon-induced nonlinear enhancements are also predicted to yield large increases in high-harmonic efficiencies beyond the third.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Simulations of the IR and Raman spectra of water confined in amorphous silica slit pores

Water in nano-scale confining environments is a key element in many biological, material, and geological systems. The structure and dynamics of the liquid can be dramatically modified under these conditions. Probing these changes can be challenging, but vibrational spectroscopy has emerged as a powerful tool for investigating their behavior. A critical, evolving component of this approach is a detailed understanding of the connection between spectroscopic features and molecular-level details. In this paper, this issue is addressed by using molecular dynamics simulations to simulate the linear infrared (IR) and Raman spectra for isotopically dilute HOD in D 2 O confined in hydroxylated amorphous silica slit pores. The effect of slit-pore width and hydroxyl density on the silica surface on the vibrational spectra is also investigated. The primary effect of confinement is a blueshift in the frequency of OH groups donating a hydrogen bond to the silica surface. This appears as a slight shift in the total (measurable) spectra but is clearly seen in the distance-based IR and Raman spectra. Analysis indicates that these changes upon confinement are associated with the weaker hydrogen-bond accepting properties of silica oxygens compared to water molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The origin of the Stokes–Einstein relation in simple dense liquids

Here, we investigate the origin of the universal relation between structural relaxation and diffusion in simple dense liquids, known as the Stokes–Einstein (SE) relation. The fact that this relation, originally derived from a hydrodynamic model of a macroscopic particle in a viscous medium, can describe the microscopic-scale liquid dynamics still eludes understanding. We introduce a new universal measure of structural relaxation in a system of N identical particles based on an explicit decomposition of the configuration space into N! congruent convex polyhedra. This measure makes it possible to quantify the correlation between two distinct particle configurations in terms of their minimal Euclidean distance, optimized with respect to particle permutations. Using this measure alongside a model of independent random walkers under the single-occupancy constraint, we derive a master equation that quantifies the SE relation. It allows us to demonstrate that the universal relation between structural relaxation and diffusion in simple dense liquids is caused by two conditions: (a) the confinement of the dominant density fluctuations to the first coordination shell, manifested by de Gennes narrowing, and (b) Gaussianity of the diffusion process; the former is shown to be violated in low-density fluids, and the latter is known to be violated in supercooled liquids.

Physics - Condensed matter physics↗

Mixed stochastic-deterministic density functional theoretic decomposition of Kubo–Greenwood conductivities in the projector augmented wave formalism

Pairing the accuracy of the Kohn–Sham density-functional framework with the efficiency of a stochastic algorithmic approach, mixed stochastic-deterministic density functional theory (mDFT) achieves a favorable computational scaling with system sizes and electronic temperatures. We employ the recently developed mDFT formalism to investigate the dynamic charge-transport properties of systems in the warm dense matter regime. The optical conductivity spectra are computed for single- and multi- component mixtures of carbon, hydrogen, and beryllium using two complementary approaches: Kubo–Greenwood in the mDFT picture and real-time time-dependent mDFT. We further devise a decomposition of the Onsager coefficients leading up to the Kubo–Greenwood spectra to exhibit contributions from the deterministic, stochastic, and mixed electronic state transitions at different incident photon energies.

36 MATERIALS SCIENCE↗

Inertia-induced mixing and reaction maximization in laminar porous media flows

Solute transport and biogeochemical reactions in porous and fractured media flows are controlled by mixing, as are subsurface engineering operations such as contaminant remediation, geothermal energy production, and carbon sequestration. Porous media flows are generally regarded as slow, so the effects of fluid inertia on mixing and reaction are typically ignored. Here, we demonstrate through microfluidic experiments and numerical simulations of mixing-induced reaction that inertial recirculating flows readily emerge in laminar porous media flows and dramatically alter mixing and reaction dynamics. An optimal Reynolds number that maximizes the reaction rate is observed for individual pore throats of different sizes. This reaction maximization is attributed to the effects of recirculation flows on reactant availability, mixing, and reaction completion, which depend on the topology of recirculation relative to the boundary of the reactants or mixing interface. Recirculation enhances mixing and reactant availability, but a further increase in flow velocity reduces the residence time in recirculation, leading to a decrease in reaction rate. The reaction maximization is also confirmed in a flow channel with grain inclusions and randomized porous media. Interestingly, the domain-wide reaction rate shows a dramatic increase with increasing Re in the randomized porous media case. This is because fluid inertia induces complex three-dimensional flows in randomized porous media, which significantly increases transverse spreading and mixing. This study shows how inertial flows control reaction dynamics at the pore scale and beyond, thus having major implications for a wide range of environmental systems.

Chen, Michael A. (ORCID:0000000294614710)↗

Enhancing Graphene Plasmonic Device Performance via its Dielectric Environment

Graphene plasmons provide a compelling avenue toward chip-scale dynamic tuning of infrared light. Dynamic tunability emerges through controlled alterations in the optical properties of the system defining graphene’s plasmonic dispersion. Typically, electrostatic induced alterations of the carrier concentration in graphene working in conjunction with mobility have been considered the primary factors dictating plasmonic tunability. We find here that the surrounding dielectric environment also plays a primary role, dictating not just the energy of the graphene plasmon but so too the magnitude of its tuning and spectral width. To arrive at this conclusion, poles in the imaginary component of the reflection coefficient are used to efficiently survey the effect of the surrounding dielectric on the tuning of the graphene plasmon. By investigating several common polar materials, optical phonons (i.e., the Reststrahlen band) of the dielectric substrate are shown to appreciably affect not only the plasmon’s spectral location but its tunability, and its resonance shape as well. In particular, tunability is maximized when the resonances are spectrally distant from the Reststrahlen band, whereas sharp resonances (i.e., high-Q) are achievable at the band’s edge. Overall, these observations both underscore the necessity of viewing the dielectric environment in aggregate when considering the plasmonic response derived from two-dimensional materials and provide heuristics to design dynamically tunable graphene-based infrared devices.

36 MATERIALS SCIENCE↗

Compiler and Runtime Approaches to Enable Large-Scale Irregular Programs. Final report, July 2013 - July 2019

While regular algorithms, characterized by operations on dense matrices and arrays, have long been the mainstay of scientific, high-performance computing, irregular algorithms, which feature unpredictable accesses to pointer-based data structures, are becoming increasingly common in high performance computing, arising in graph analysis, data mining and visualization, among other domains. Unfortunately, the defining characteristics of irregular applications, their dynamic, unpredictable, data-dependent access patterns and data layouts, make achieving high performance on large scale systems difficult. Scaling applications to peta- and exa-scale requires carefully controlling communication and data movement and placement, an inherently difficult task when access patterns and data layouts are unpredictable! Most irregular applications that attain high performance must be painstakingly hand-written and hand-tuned, with few common principles or paradigms uniting various implementations and easing future development. Despite the increasing importance of irregular applications, there is little programmer knowledge, and even less compiler ability, devoted to optimizing them. This project aims to solve these problems. By allowing programmers to write irregular applications in high level forms, with at most a few annotations highlighting key structural properties, programmers can focus on developing their algorithms and methods. The compiler and run-time system can take on the tedious task of optimizing the application for execution at large scales, and can automatically provide efficient implementations. This will provide portability and ease maintenance for existing irregular applications, but, more importantly, open up whole new domains of computational science to large-scale, high-performance simulation codes.

97 MATHEMATICS AND COMPUTING↗

Atomic Ejecta Source Optical Probe (AESOP) (L3 Milestone Report FY2021 Work Package)

As part of its mission to ensure the reliability of the nation's nuclear weapons stockpile, the NNSA has funded a broad range of projects to characterize the mass of material ejected from shocked metal surfaces with special emphasis placed on determining the size distribution of macroscopic "ejecta" particles. Substantial work has been done for particles with diameters down to roughly 1 μm, but little data is available for smaller particles, and no data exist for the amount of atomic-scale ejecta produced in such dynamic environments. Such data are important due to the implications for weapons systems behavior but also for distinguishing between different potential ejecta production mechanisms, such as Rayleight-Taylor instability or the shallow-bubble-collapse phenomenon recently proposed by G. Maskally.

36 MATERIALS SCIENCE↗

MAGNET Scaling and Methodology

The purpose of this study was to analyze the heat transfer of the Microreactor Agile Non-nuclear Experimental Testbed (MAGNET) within the Dynamic Energy Transport and Integration Laboratory (DETAIL) and develop scaling equations and models to couple with other systems. Hierarchical Two-Tiered Scaling (H2TS) methodologies were applied to DETAIL’s MAGNET facility to scale and project data sets while conserving the observed behavior based on first principles. The MAGNET system was successfully scaled using H2TS and multiple system parameters were determined or calculated from experimental data including steady state and transient data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Mental workload in decision and control

This paper briefly reviews the problems of defining and measuring the 'mental workload' of aircraft pilots and other human operators of complex dynamic systems. Of the alternative approaches the author indicates a clear preference for the use of subjective scaling. Some recent experiments from MIT and elsewhere are described which utilize subjective mental workload scales in conjunction with human decision and control tasks in the laboratory. Finally a new three-dimensional mental workload rating scale, under current development for use by IFR aircraft pilots, is presented.

Sheridan, T. B.↗

A numerical experiment on the equilibrium and stability of a rotating galactic bar

A self-consistent, three-dimensional numerical experiment is performed on an N-body system whose initial state is a realization of a certain theoretical model of a rotating triaxial galaxy. The model is a stellar-dynamical counterpart of a uniformly rotating polytrope of index equal to 0.5. The aim of the experiment is to study the equilibrium of the system and, in particular, to test its stability. The experimental system behaves in the mean like a realization of the theoretical model for at least seven crossing times. The principal departure of the system from equilibrium is an oscillation which is identified as a radial pulsation. There is no indication in its behavior that the system is unstable with respect to anu mode with an e-folding time shorter than or of the order of two crossing times. Certain changes that occur in the state of the system are interpreted, with the aid of the theoretical model, as secular changes which result from a slight failure of our numerical methods to conserve the mass, energy, and angular momentum of the system; these effects are small enough that they do not vitiate the experiment on a dynamical time scale.

Miller, R. H.↗

Control law parameterization for an aeroelastic wind-tunnel model equipped with an active roll control system and comparison with experiment

Nominal roll control laws were designed, implemented, and tested on an aeroelastically-scaled free-to-roll wind-tunnel model of an advanced fighter configuration. The tests were performed in the NASA Langley Transonic Dynamics Tunnel. A parametric study of the nominal roll control system was conducted. This parametric study determined possible control system gain variations which yielded identical closed-loop stability (roll mode pole location) and identical roll response but different maximum control-surface deflections. Comparison of analytical predictions with wind-tunnel results was generally very good.

Perry, Boyd, III↗

Control law parameterization for an aeroelastic wind-tunnel model equipped with an active roll control system and comparison with experiment

Nominal roll control laws were designed, implemented, and tested on an aeroelastically-scaled free-to-roll wind-tunnel model of an advanced fighter configuration. The tests were performed in the NASA Langley Transonic Dynamics Tunnel. A parametric study of the nominal roll control system was conducted. This parametric study determined possible control system gain variations which yielded identical closed-loop stability (roll mode pole location) and identical roll response but different maximum control-surface deflections. Comparison of analytical predictions with wind-tunnel results was generally very good.

Perry, Boyd, III↗

High-speed inlet research program and supporting analysis

The technology challenges faced by the high speed inlet designer are discussed by describing the considerations that went into the design of the Mach 5 research inlet. It is shown that the emerging three dimensional viscous computational fluid dynamics (CFD) flow codes, together with small scale experiments, can be used to guide larger scale full inlet systems research. Then, in turn, the results of the large scale research, if properly instrumented, can be used to validate or at least to calibrate the CFD codes.

Coltrin, Robert E.↗

Extension to the dynamic modeling of the large angle magnetic suspension test fixture

The Large Angle Magnetic Suspension Test Fixture (LAMSTF) is a laboratory scale proof-of-concept system. The configuration is unique in that the electromagnets are mounted in a circular planar array. A mathematical model of the system had previously been developed, but was shown to have inaccuracies. These inaccuracies showed up in the step responses. Eddy currents were found to be the major cause of the modeling errors. In the original system, eddy currents existed in the aluminum baseplate, iron cores, and the sensor support frame. An attempt to include the eddy current dynamics in the system model is presented. The dynamics of a dummy sensor ring were added to the system. Adding the eddy current dynamics to the simulation improves the way it compares to the actual experiment. Also presented is a new method of determining the yaw angle of the suspended element. From the coil currents the yaw angle can be determined and the controller can be updated to suspend at the new current. This method has been used to demonstrate a 360 degree yaw angle rotation.

Foster, Lucas E.↗

Nanoscale Etching and Indentation of Silicon Surfaces with Carbon Nanotubes

The possibility of nanolithography of silicon and germanium surfaces with bare carbon nanotube tips of scanning probe microscopy devices is considered with large scale classical molecular dynamics (MD) simulations employing Tersoff's reactive many-body potential for heteroatomic C/Si/Ge system. Lithography plays a key role in semiconductor manufacturing, and it is expected that future molecular and quantum electronic devices will be fabricated with nanolithographic and nanodeposition techniques. Carbon nanotubes, rolled up sheets of graphene made of carbon, are excellent candidates for use in nanolithography because they are extremely strong along axial direction and yet extremely elastic along radial direction. In the simulations, the interaction of a carbon nanotube tip with silicon surfaces is explored in two regimes. In the first scenario, the nanotubes barely touch the surface, while in the second they are pushed into the surface to make "nano holes". The first - gentle scenario mimics the nanotube-surface chemical reaction induced by the vertical mechanical manipulation of the nanotube. The second -digging - scenario intends to study the indentation profiles. The following results are reported in the two cases. In the first regime, depending on the surface impact site, two major outcomes outcomes are the selective removal of either a single surface atom or a surface dimer off the silicon surface. In the second regime, the indentation of a silicon substrate by the nanotube is observed. Upon the nanotube withdrawal, several surface silicon atoms are adsorbed at the tip of the nanotube causing significant rearrangements of atoms comprising the surface layer of the silicon substrate. The results are explained in terms of relative strength of C-C, C-Si, and Si-Si bonds. The proposed method is very robust and does not require applied voltage between the nanotube tips and the surface. The implications of the reported controllable etching and hole-creating for nanolithography on silicon are discussed in detail.

Dzegilenko, Fedor N.↗

Simulations of Li+ in Ionic Liquids: Structure, Transport, and Electrochemical Windows

Ionic liquids have been proposed as candidate electrolytes for a number of electrochemical applications. The Li+ solvation structure in these liquids is of central importance to electrolyte properties, like ionic conductivity and electrochemical stability. To this point, we employ simulations at three different size scales to better understand various aspects of the interplay between Li+ solvation structure and dynamics. The smallest systems are Li(Anion)n clusters that are treated with high-accuracy density functional theory (DFT) techniques to provide insight into solvation shell structure through energetics and comparisons to experimental IRRaman spectra. Mid-range sized liquid-phase systems (12-24 ion pairs) are treated with DFT molecular dynamics (MD) to provide temperature-dependent insight into Li+ solvation structure, diffusion, and electrochemical window. The largest systems (144-216 ion pairs) are treated with polarizable MD simulations to evaluate the influence of Li-networks on structure and provide size independent values of transport properties. We perform this procedure on three technologically important ionic liquids and comment on property correlations with solvation structure.

batteries↗