Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “characteristic equations”

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 235 records · Page 13

Observation of critical scaling in spin glasses below T c using thermoremanent magnetization

Time-dependent thermoremanent magnetization (TRM) studies have been instrumental in probing energy dynamics within the spin glass phase. In this paper, we review the evolution of the TRM experiment over the last half century and discuss some aspects related to how it has been used in the understanding of spin glasses. We also report on recent experiments using high-resolution DC SQUID magnetometry to probe the TRM at temperatures less than but near to the transition temperature T c . These experiments have been performed as a function of waiting time, temperature, and five different magnetic fields. We find that as the transition temperature is approached from below, the characteristic time scale of TRM is suppressed up to several orders of magnitude in time. In the highest-temperature region, we find that the waiting time effect subsides, and a waiting time-independent crossover line is reached. We also find that increasing the magnetic field further suppresses the crossover line. Using a first-principles energy argument across the crossover line, we derive an equation that is an excellent fit to the crossover lines for all magnetic fields probed. The data show strong evidence for critical slowing down and an H = 0 Oe phase transition.

36 MATERIALS SCIENCE↗

Analysis of the weighted shifted boundary method for the Poisson and Stokes problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. Accuracy is maintained by properly shifting the location and values of the boundary conditions. This avoids integration over cut cells and the associated implementation issues. Recently, the Weighted SBM (WSBM) was proposed for the Navier-Stokes equations with free surfaces and the Stokes flow with moving boundaries. The attribute “weighted” in the name WSBM stems from the fact that its variational form is weighted with the elemental volume fraction of active fluid. The motivation for the development of the WSBM was the preservation of the volume of active fluid to a higher degree of accuracy, which in turn resulted in improved stability and robustness characteristics in moving-boundary, time-dependent simulations. In this article, we present the numerical analysis of the WSBM formulations for the Poisson and Stokes problems. We give mathematical conditions under which the bilinear forms defining the discrete variational formulations are uniformly coercive (Poisson problem) or inf-sup stable (Stokes problem). By these results, stability and optimal convergence is proven in the natural norm; L2-error estimates can also be derived.

Approximate domain boundaries↗

Statistical description of coalescing magnetic islands via magnetic reconnection

The physical picture of interacting magnetic islands provides a useful paradigm for certain plasma dynamics in a variety of physical environments, such as the solar corona, the heliosheath and the Earth's magnetosphere. In this work, we derive an island kinetic equation to describe the evolution of the island distribution function (in area and in flux of islands) subject to a collisional integral designed to account for the role of magnetic reconnection during island mergers. This equation is used to study the inverse transfer of magnetic energy through the coalescence of magnetic islands in two dimensions. We solve our island kinetic equation numerically for three different types of initial distribution: Dirac delta, Gaussian and power-law distributions. The time evolution of several key quantities is found to agree well with our analytical predictions: magnetic energy decays as $\tilde {t}^{-1}$ , the number of islands decreases as $\tilde {t}^{-1}$ and the averaged area of islands grows as $\tilde {t}$ , where $\tilde {t}$ is the time normalised to the characteristic reconnection time scale of islands. General properties of the distribution function and the magnetic energy spectrum are also studied. Finally, we discuss the underlying connection of our island-merger models to the (self-similar) decay of magnetohydrodynamic turbulence.

Physics↗

Electrode-molecule energy level offsets in a gold-benzene diamine-gold single molecule tunnel junction

One means for describing electron transport across single molecule tunnel junctions (MTJs) is to use density functional theory (DFT) in conjunction with a nonequilibrium Green’s function formalism. This description relies on interpreting solutions to the Kohn–Sham (KS) equations used to solve the DFT problem as quasiparticle (QP) states. Many practical DFT implementations suffer from electron self-interaction errors and an inability to treat charge image potentials for molecules near metal surfaces. For MTJs, the overall effect of these errors is typically manifested as an overestimation of electronic currents. Correcting KS energies for self-interaction and image potential errors results in MTJ current–voltage characteristics in close agreement with measured currents. An alternative transport approach foregoes a QP picture and solves for a many-electron wavefunction on the MTJ subject to open system boundary conditions. It is demonstrated that this many-electron method provides similar results to the corrected QP picture for electronic current. The analysis of these two distinct approaches is related through corrections to a junction’s electronic structure beyond the KS energies for the case of a benzene diamine molecule bonded between two gold electrodes.

Szepieniec, Mark S.↗

Simultaneous velocity and density measurements of fully developed Rayleigh-Taylor mixing

Here, the dynamics of molecular mixing and the energy transfer process in the Rayleigh-Taylor instability (RTI) are studied through the collection of simultaneous density-velocity measurements. These experiments provide simultaneous density-velocity field measurements, in contrast to previous point measurements. Statistically stationary experiments are performed in a “convective-type” gas tunnel facility, with density contrast achieved through the injection of helium into the bottom stream. Three experiments at Atwood number ≈ 0.1 are captured at three outer scale Reynolds numbers Re = 520, 2260, and 4050. Particle image velocimetry and laser induced fluorescence are employed simultaneously. Statistics of the density and velocity show self-similar collapse of RTI profiles at large Reynolds number Re >2000. Flat velocity profiles indicate homogeneous turbulence characteristics in the core of the mixing region. Significant anisotropy develops in the flow, with horizontal velocity fluctuations being only 60% of the vertical velocity fluctuations. The turbulent mass flux is found to be asymmetric about the centerline, with increased peak towards the spike. Measurements of the molecular mixing show that mixing is maximized at the core of the flow and increases with increased Reynolds number. The transport equation of density-specific-volume correlation b shows that it is mostly produced in the core of the mixing region, and that the spatial evolution of its profile is the result of transport by bulk motion of the bubble and spike. Energy transfer from gravitational potential energy to turbulent kinetic energy and viscous dissipation is observed to occur in the experiment with a ratio of dissipated energy to potential energy released of 38%. The analysis of the turbulent kinetic energy transport equation budget reveals that production is the primary mechanism towards the growth of turbulent kinetic energy in the core of the flow, and is asymmetrically slightly skewed towards the spike. However, it is through the transport that the strong advection at the edges of the mixing region is maintained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Formation of turing patterns in strongly magnetized electric discharges

Abstract Pattern formation and self-organization in many biological and non-biological systems can be explained through Turing’s activator-inhibitor model. Here we show how this model can be employed to describe the formation of filamentary structures in a low-pressure electric discharge exposed to a strong magnetic field. Theoretical investigation reveals that the fluid equations describing a magnetized plasma can be rearranged to take the mathematical form of Turing’s activator-inhibitor model. Numerical simulations based on the equations derived from this approach could reproduce the various patterns observed in the experiments. Also, it is shown that a density imbalance between electrons and ions exists in the bulk of the magnetized plasma that generates an electric field structure transverse to the applied magnetic field. This electric field is responsible for the stability of the filamentary patterns in the magnetized plasma over time scales much longer than the characteristic time scales of the electric discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Interpretation of wide-angle x-ray diffraction patterns from mesophase pitch-based carbon fibers – a simulation and experimental study

Carbon fiber is a critical material in a wide range of industries, where it is highly valued for its high specific strength/stiffness, excellent wear resistance, efficient electrical and thermal transport properties, chemical resistance, and low coefficient of thermal expansion. The properties of a specific carbon fiber are closely tied to its structural characteristics at all length scales. Here, in this work, we applied wide-angle x-ray diffraction to a set of heat-treated mesophase pitch-based carbon fibers, with the goal of elucidating the crystalline structures as a function of fiber orientation. To assist with analysis and interpretation of the experimental data, we employed diffraction pattern simulations using the scalar and vector forms of the Debye scattering equation to determine the influence of basal plane orientation, crystalline ordering (turbostratic-graphitic), and basal plane asymmetry on the diffraction patterns. The results presented here suggest that growth of the transverse crystallites in mesophase pitch-based carbon fiber is fixed until graphitization temperatures are reached. The work completed here provides a framework for the analysis of carbon fiber and other oriented carbon-based materials via diffraction.

36 MATERIALS SCIENCE↗

Empirical radius formulas for canonical neutron stars from bidirectionally selecting features of equations of state in extended Bayesian analyses of observational data

Significant advancement in Bayesian inference of nuclear equation of state (EOS) from gravitational wave and x-ray observations of neutron stars (NSs) has been made by the nuclear astrophysics community especially since GW170817. By extending the traditional Bayesian analysis which normally ends at presenting the marginalized posterior probability distribution functions (PDFs) of individual EOS parameters and their correlations (or sometimes only the Pearson correlation coefficients which are only reliably useful when the variables are linearly correlated while they are actually often not), we search for a data-driven and robust empirical formula for the radius 𝑅 1.4 of canonical NSs in terms of the characteristic EOS parameters (features). We also identify the single most important but currently poorly known EOS parameter for determining the 𝑅 1.4 . Using three regression-model-building methodologies: bidirectional stepwise feature selection, least absolute shrinkage selection operator (LASSO) regression, and neural network regression on a large set of posterior EOSs and the corresponding 𝑅 1.4 values inferred from earlier comprehensive Bayesian analyses of NS observational data, we systematically and rigorously develop the most probable 𝑅 1.4 formulas with varying statistical accuracy and technical complexity. Here, the most important EOS parameters for determining 𝑅 1.4 are found consistently in each of the feature selection processes to be (in order of decreasing importance): curvature 𝐾 sym , slope 𝐿, skewness 𝐽 sym of nuclear symmetry energy, skewness 𝐽 0 , incompressibility 𝐾 0 of symmetric nuclear matter, and the magnitude 𝐸 sym ⁡(𝜌 0 ) of symmetry energy at the saturation density 𝜌 0 of nuclear matter.

Bayesian methods↗

Experimental observation and numerical investigation of filamentary structures in magnetized plasmas

Filamentary structures in low-pressure, low-temperature plasmas are produced when strong magnetic fields are applied parallel to the electric field defined by parallel electrodes. Filamentary structures are regions within the plasma that have distinct properties such as optical brightness and extend along the magnetic field lines. In our experiments, an argon, radio frequency discharge is exposed to a strong background magnetic field in the magnetized dusty plasma experiment at Auburn University. Different forms of filamentary structures emerge in the plasma such as columns, target, or spiral-like structures. To investigate the origin and the characteristics of these patterns, we have developed a three-dimensional fluid model that can reproduce the experimental observations which enables us to investigate the underlying mechanisms of the filamentation process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fluctuating hydrodynamics of chiral active fluids

Active materials are characterized by continuous injection of energy at the microscopic level and typically cannot be adequately described by equilibrium thermodynamics. Here we study a class of active fluids in which equilibrium-like properties emerge when fluctuating and activated degrees of freedom are statistically decoupled, such that their mutual information is negligible. We analyse three paradigmatic systems: chiral active fluids composed of spinning frictional particles that are free to translate, oscillating granular gases and active Brownian rollers. In all of these systems, a single effective temperature generated by activity parameterizes both the equation of state and the emergent Boltzmann statistics. The same effective temperature, renormalized by velocity correlations, relates viscosities to steady-state stress fluctuations via a Green–Kubo relation. To rationalize these observations, we develop a theory for the fluctuating hydrodynamics of these non-equilibrium fluids and validate it through large-scale molecular dynamics simulations. Our work sheds light on the microscopic origin of odd viscosities and stress fluctuations characteristic of parity-violating fluids, in which mirror symmetry and detailed balance are broken.

74 ATOMIC AND MOLECULAR PHYSICS↗

Low-Memory, Discrete Ordinates, Discontinuous Galerkin Methods for Radiative Transport

The discrete ordinates discontinuous Galerkin (SN-DG) method is a well-established and practical approach for solving the radiative transport equation. In this paper, we study a low-memory variation of the upwind SN-DG method. The proposed method uses a smaller finite element space that is constructed by coupling spatial unknowns across collocation angles, thereby yielding an approximation with fewer degrees of freedom than the standard method. Like the original SN-DG method, the low-memory variation still preserves the asymptotic diffusion limit and maintains the characteristic structure needed for mesh sweeping algorithms. While we observe second-order convergence in the scattering dominated, diffusive regime, the low-memory method is in general only first-order accurate. To address this issue, we use upwind reconstruction to recover second-order accuracy. Finally, for both methods, numerical procedures based on upwind sweeps are proposed to reduce the system dimension in the underlying Krylov solver strategy.

97 MATHEMATICS AND COMPUTING↗

Equations of Motion for the Vertical Rigid-Body Rotor: Linear and Nonlinear Cases

Centuries ago, the prolific mathematician Leonhard Euler (1707–1783) wrote down the equations of motion (EOM) for the heavy symmetrical top with one point fixed. The resulting set of equations turned out to be nonlinear and had a limited number of closed-form solutions.Today, tools such as transfer matrix and finite elements enable the calculation of the rotordynamic properties for rotor-bearing systems. Some of these tools rely on the “linearized” version of the EOM to calculate the eigenvalues, unbalance response, or transients in these systems.In fact, industry standards mandate that rotors be precisely balanced to have safe operational characteristics. However, in some cases, the nonlinear aspect of the EOM should be considered.The purpose of this chapter is to show examples of how the linear vs. nonlinear formulations differ. This chapter also shows how excessive unbalance is capable of dramatically altering the behavior of the system and can produce chaotic motions associated with the “jump” phenomenon.

Laos, Hector↗

An integrated approach to optimizing concentration shock wave electrodialysis using 2D multicell simulation and response surface models

Shock wave electrodialysis (SWED) is a highly promising technique for energy-efficient ion separation in the context of a circular economy. This paper presents a approach way of modeling and improving SWED using a two-dimensional multicell model combined with the COMSOL program and response surface methodology. The model integrates the Nernst-Planck equation, Darcy's law, and first-order electroosmosis to examine the local concentration, flux of ionic species, distribution of current, and velocity of flow in SWED cells under various operating conditions. We first illustrate the clear depiction of concentration, velocity, and electric potential distribution through contours which aids in identifying optimal operating conditions and designing scalable SWED systems. The results emphasize the significance of surface charge density and voltage in influencing the features of shock waves for obtaining effective ion separation while optimizing energy consumption and improving current efficiency by controlling the retention time of feed flow. Here, this study defines two crucial characteristics of shock waves, namely the length of the flat depletion zone of a fully developed shock wave (shock wave height) and the distance of shock wave propagation (shock wave length). These properties significantly impact separation performance, as determined by the simulation results. Additionally, the response surface methodology is incorporated with the COMSOL models to develop predictive models and graph responses, enabling a more comprehensive understanding of the interactions between parameters and performance indicators, such as removal ratio, energy consumption, and water recovery. Finally, this work suggests design tactics for expanding SWED processes and outlines potential areas for further research. This research provides valuable insights into the prospective applications, design optimization, and scalability of SWED in the field of electrokinetic separation technologies for green chemistry and a circular economy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accelerating Biomimetic Solar - Energy Harvesting: Mapping the Interaction Landscape of Plasmonic-Excitonic Hybrid Nanosystems (Final Report)

In general, excitonic and plasmonic nanoscale materials in close proximity show high potential for significant breakthroughs in energy related materials research. The interactions between these two kinds of materials result in coupled optical transitions (plexcitons), distinct from those of both the individual exciton and plasmon as well as from those of the sum of their constituents (synergistic effects). By linking together materials-research and physical-research approaches, this project contributes to a concerted approach on nanomaterials energy research. The project’s overall goal is to accelerate the development of well-defined plexcitonic model systems consisting of carefully engineered plasmonic and excitonic nanomaterial— essential for both gaining a fundamental understanding of plexcitonic nanomaterials and the development of novel design principles for biomimetic solar energy harvesting. During the 3-year project period and the terminal renewal with limited support for a 12-month period, we successfully synthesized and characterized (1) a robust excitonic nanomaterial and (2) a library of plasmonic nanoparticles as well as developed (3) a microfluidic platform for homogenous nanosynthesis as summarized below: (1) Robust Excitonic Nanomaterial. Supramolecular assemblies are Nature’s most successful material system for solar energy harvesting. However, photovoltaic devices based on artificial supramolecular assemblies continue to be stymied by disappointing efficiencies and poor stability. The conceptual failure may lie in current solar cell architectures, which rely on solidifying supramolecular assemblies as an ensemble into a solid matrix, neglecting the intrinsic fragility of the assemblies’ internal structure, thus disrupting or even destroying their delicate optoelectronic properties, that is, delicate Frenkel excitonic properties. Supramolecular assemblies may finally serve as usable light harvesting material systems for solar energy conversion technologies, only if they meet the following criteria: (a) Stability, that is, the fragile structure including its delicate Frenkel excitonic character needs to be stable, (b) Robustness, that is, resistant against elevated and fluctuating temperatures, and (c) Viability for device integration, that is, capable of being immobilized onto solid substrates. Here, by developing a nanocomposite via a tunable, cage-like scaffold design, we successfully provided stable supramolecular nanocomposites, that inhabit robust Frenkel excitons despite harming environmental conditions such as extreme heat stress. (2) Library of Plasmonic Nanoparticles. Naturally, current models describing plasmonic hybrid quantum states—plasmonic hybridizations—parallel those developed for molecular orbitals, equating individual plasmonic nanostructures with “atoms” and the plasmonic nanoassemblies with “molecules.” In analogy to organic synthesis, a suitably robust fabrication method would allow for “atom-like” manipulation of “molecule-like” plasmonic nanoassemblies; of high value for next-generation energy nanotechnologies. Despite this frequent comparison, current plasmonic nanoassembly fabrication methods favor top-down templating over wet-chemical synthesis, however, achieving precise control over nanostructure’s geometry and surface characteristics remain an art and a scientific challenge. The conceptual failure may lie in the current wet-chemical synthesis paradigm, as it relies on the accessibility of a multi-dimensional synthesis parameter space through limited, rather one-dimensional synthesis procedures by employing step-by-step approaches. Solution-based nanoarchitectonics for rational design of precisely built plasmonic nanoassemblies via solution-based fabrication may finally be possible only if multi-dimensional syntheses approaches are available that allow for comprehensive control over the plasmonic nanomaterials’ (a) Structural Properties and (b) Surface Properties. Here, by developing an innovative multidimensional 1,3-propanediol based polyol synthesis, we successfully provided control over the plasmonic building-block’s geometry (size and shape) together with its surface characteristics. Our results present a critical step toward the vision of a “periodic table-like” system for plasmonic materials based on straightforward wet-chemical syntheses for energy nanotechnologies. Developing deliberate modifications on this synthesis, we generated a library of plasmonic nanostructures covering the vast parameter space—opening the door for fundamental investigation of plexcitonic model systems. (3) Microfluidic Platform for Homogenous Nanosynthesis. Control over structural properties of plexcitonic nanocomposites remains a challenge due to current limitations in nanosynthesis techniques. Slight variations in nanostructure’s geometry impact their optoelectronic properties, demanding precise synthesis beyond the capabilities of solution-based (batch) synthesis processes. In contrast, the small, confined liquid volumes used in microfluidics—a reaction technique where the manipulation of fluids takes place in channels with dimensions of tens of micrometers—allows for homogenous synthesis conditions, providing excellent control of the reaction and, as a result, of the materials’ geopmetry and composition. However, thus far, the majority of plexcitonic systems has been developed via batch synthesis. Here, by successfully developing a two-channel microreactor, our microfluidic-supported synthesis approach combines the advantages of both microfluidics and batch platforms, allowing for precise spatio-temporal control over all synthesis parameters opening the possibility for homogenous nanosythnesis of well-defined plexcitonic model systems.

14 SOLAR ENERGY↗

Numerical study of proppant transport and settling processes in fractures

Reservoir stimulation by creating hydraulically conductive fractures is the key step for enabling enhanced geothermal systems (EGS). The effectiveness of stimulation is significantly influenced by the deposition of proppant inside induced fractures. The transportation and settling of proppant in a propagating fracture is controlled by a multitude of operational and physical parameters, including the fracturing fluid rheology, injection rate, proppant concentration, fracture length/aperture evolution, proppant size/density/shape, etc. A numerical tool that robustly and efficiently accounts for all important attributes can facilitate the design and optimization of reservoir stimulation. This study presents the novel computational tool ELK (ELectrical fracKing) developed for the numerical simulation of proppant-fluid mixture circulation in a fractured geothermal reservoir. We enriched the MOOSE-based PorousFlow module with a suite of equations to consider the fluid-proppant mixture with particle-particle/fluid interactions, which include gravitational settling, particle convection, particle hampering, and strong density and viscosity contrasts. The computational tool is validated by comparing the predicted proppant bed evolution against two different laboratory scale experiments of proppant transport in a fixed aperture channel. Further parameter studies were performed, and the modeling results show that the proppant deposition is determined by the mixing characteristics and settling of the particles from the slurry. Concentration-dependent density and viscosity lead to an inhomogeneous distribution of the proppant, particle collision, and enhanced settling at the bottom of the fractures. Preliminary coupling with dynamic fracture propagation shows promising results and will be further developed to simulate hydraulic stimulation at high fidelity.

15 GEOTHERMAL ENERGY↗

Computational and numerical analysis of AC optimal power flow formulations on large-scale power grids

We report that alternating current optimal power flow (AC-OPF) is a fundamental tool in electric utilities to determine optimal operation of the various resources. Typically, the AC-OPF problem uses power balance formulation containing voltages and power equations. Yet, there is no comprehensive comparison of the different AC-OPF formulations, especially for large-scale networks. This paper presents a detailed comparative evaluation of different formulations of the AC-OPF problem on networks ranging from 9-bus to 25,000 buses. Three different formulations: 1) power balance with polar voltages, 2) power balance with Cartesian voltages, and 3) current balance with Cartesian voltages are discussed in detail by comparing their characteristics, and numerical and computational performance.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Equations of Motion for the Vertical Rigid-Body Rotor: Linear and Nonlinear Cases

Centuries ago, the prolific mathematician Leonhard Euler (1707–1783) wrote down the equations of motion (EOM) for the heavy symmetrical top with one point fixed. The resulting set of equations turned out to be nonlinear and had a limited number of closed-form solutions. Today, tools such as transfer matrix and finite elements enable the calculation of the rotor dynamic properties for rotor-bearing systems. Some of these tools rely on the “linearized” version of the EOM to calculate the eigenvalues, unbalance response, or transients in these systems. In fact, industry standards mandate that rotors be precisely balanced to have safe operational characteristics. However, in some cases, the nonlinear aspect of the EOM should be considered. The purpose of this paper is to show examples of how the linear vs. nonlinear formulations differ. This paper will also show how excessive unbalance is capable of dramatically altering the behavior of the system and can produce chaotic motions associated with the “jump” phenomenon.

42 ENGINEERING↗

Adaptive Data-Driven Deep-Learning Surrogate Model for Frontal Polymerization in Dicyclopentadiene

Frontal polymerization (FP) is a self-sustaining curing process that enables rapid and energy-efficient manufacturing of thermoset polymers and composites. Computational methods conventionally used to simulate the FP process are time-consuming, and repeating simulations are required for sensitivity analysis, uncertainty quantification, or optimization of the manufacturing process. Here, in this work, we develop an adaptive surrogate deep-learning model for FP of dicyclopentadiene (DCPD), which predicts the evolution of temperature and degree of cure orders of magnitude faster than the finite-element method (FEM). The adaptive algorithm provides a strategy to select training samples efficiently and save computational costs by reducing the redundancy of FEM-based training samples. The adaptive algorithm calculates the residual error of the FP governing equations using automatic differentiation of the deep neural network. A probability density function expressed in terms of the residual error is used to select training samples from the Sobol sequence space. The temperature and degree of cure evolution of each training sample are obtained by a 2D FEM simulation. The adaptive method is more efficient and has a better prediction accuracy than the random sampling method. With the well-trained surrogate neural network, the FP characteristics (front speed, shape, and temperature) can be extracted quickly from the predicted temperature and degree-of-cure fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗