Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Characterizing Variability in Lignocellulosic Biomass: A Review

Feedstock variability is a significant barrier to the scale-up and commercialization of lignocellulosic biofuel technologies. Variability in feedstock characteristics and behavior creates numerous challenges to the biorefining industry by affecting continuous operation and biofuels yields. Recently, feedstock variability is understood and explained largely on the basis of chemical composition. Physical and mechanical properties and behavior of lignocellulosic feedstock in various unit operations, studied through advanced analytical methods, can further explain variability. Such studies will enable us in developing processes and designing equipment to improve operation and conversion performance. In this perspective, we review several advanced analytical methods that measure density, moisture content, thermal properties, flowability, grindability, rheology properties, and micromorphological characteristics. We also discuss the correlations and interactions among these properties that reflect the complexity of lignocellulosic biomass as a feedstock and the associated quality metrics and logistics of supplying consistent quality feedstock to a biorefinery. We further examine methods that have not traditionally been used to characterize lignocellulosic feedstocks but have the potential to bridge the gap in our explanation of feedstock variability.

09 BIOMASS FUELS↗
DOE OSTI2022

Real-time laser ultrasonic monitoring of laser-induced thermal processes

Abstract Intra- and inter-layer integrity of components fabricated with advanced manufacturing techniques, such as laser powder bed fusion, is dependent upon rapid heating, melting, and solidification processes. There is a need for new techniques to provide in situ feedback of these processes. Here a laser-based ultrasonic technique to probe thermal effects induced by a high-power continuous wave laser in titanium samples is described. Numerical simulations were performed to show that, for a spatially uniform heating beam, laser-induced surface acoustic waves are strongly influenced by surface heating conditions, are dispersive in the case of rapid heating, and that an abrupt velocity reduction happens upon the onset of surface melting. Furthermore, laser-based ultrasound experimental results which monitor the transient change of surface wave travel time associated with high power laser surface heating are provided. A pulsed laser is used to generate high frequency surface acoustic waves that propagate through the laser-heated region and are detected using a photorefractive crystal-based interferometer. Qualitative agreement is observed between theory and experiment with both showing a rapid reduction in the surface wave velocity at the onset of illumination and further decrease in surface wave velocity associated with melting. It is demonstrated that changes in the surface wave velocity can be used to track local heating and detect the onset of surface melting in real time.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Disentangling the chemistry of Australian plant exudates from a unique historical collection

For thousands of years, the unique physicochemical properties of plant exudates have defined uses in material culture and practical applications. Native Australian plant exudates, including resins, kinos, and gums, have been used and continue to be used by Aboriginal Australians for numerous technical and cultural purposes. A historic collection of well-preserved native Australian plant exudates, assembled a century ago by plant naturalists, gives a rare window into the history and chemical composition of these materials. Here we report the full hierarchical characterization of four genera from this collection, Xanthorrhoea, Callitris, Eucalyptus, and Acacia, from the local elemental speciation, to functional groups and main molecular markers. We use high-resolution X-ray Raman spectroscopy (XRS) to achieve bulk-sensitive chemical speciation of these plant exudates, including insoluble, amorphous, and cross-linked fractions, without the limitation of invasive and/or surface specific methods. Combinatorial testing of the XRS data allows direct classification of these complex natural species as terpenoid, aromatic, phenolic, and polysaccharide materials. Differences in intragenera chemistry was evidenced by detailed interpretation of the XRS spectral features. We complement XRS with Fourier-transform infrared (FT-IR) spectroscopy, gas chromatography–mass spectrometry (GC-MS), and pyrolysis–GC-MS (Py-GC-MS). This multimodal approach provides a fundamental understanding of the chemistry of these natural materials long used by Aboriginal Australian peoples.

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

Rheology of concentrated fiber suspensions with a load-dependent friction coefficient

Here we investigate the effects of fiber aspect ratio, roughness, flexibility, and volume fraction on the rheology of concentrated suspensions in a steady shear flow using direct numerical simulations. We model the fibers as inextensible continuous flexible slender bodies with the Euler-Bernoulli beam equation governing their dynamics suspended in an incompressible Newtonian fluid. The fiber dynamics and fluid flow coupling is achieved using the immersed boundary method. In addition, the fiber surface roughness might lead to interfiber contacts, resulting in normal and tangential forces between the fibers, which follow Coulomb’s law of friction. The surface roughness is modeled as hemispherical protrusions on the fiber surfaces. Their deformation results in a normal load-dependent friction coefficient. Our simulations accurately predict the experimentally observed shear thinning in fiber suspensions. Furthermore, we find that the suspension viscosity η increases with increasing the volume fraction, roughness, fiber rigidity, and aspect ratio. The increase in η is the macroscopic manifestation of a similar increase in the microscopic contact contribution to the total stress with these parameters. In addition, we observe positive and negative first N 1 and second N 2 normal stress differences, respectively, with |N 2 | < |N 1 |, in agreement with previous experiments. Last, we propose a modified Maron-Pierce law to quantify the reduction in the jamming volume fraction by increasing the fiber aspect ratio and roughness. Our results and analysis establish the use of fiber surface tribology to tune the suspension flow behavior.

42 ENGINEERING↗
DOE OSTI2022

Laplacian-level meta-generalized gradient approximation for solid and liquid metals

In this work, we derive and motivate a Laplacian-level, orbital-free meta-generalized-gradient approximation (LL-MGGA) for the exchange-correlation energy, targeting accurate ground-state properties of sp and sd metallic condensed matter, in which the density functional for the exchange-correlation energy is only weakly nonlocal due to perfect long-range screening. Our model for the orbital-free kinetic energy density restores the fourth-order gradient expansion for exchange to the r 2 SCAN meta-GGA [Furness et al., J. Phys. Chem. Lett. 11, 8208 (2020)], yielding a LL-MGGA we call OFR2. OFR2 matches the accuracy of SCAN for prediction of common lattice constants and improves the equilibrium properties of alkali metals, transition metals, and intermetallics that were degraded relative to the PBE GGA values by both SCAN and r 2 SCAN. We compare OFR2 to the r 2 SCAN -L LL-MGGA [Mejia-Rodriguez and Trickey, Phys. Rev. B 102, 121109(R) (2020)] and show that OFR2 tends to outperform r 2 SCAN -L for the equilibrium properties of solids, but r 2 SCAN -L much better describes the atomization energies of molecules than OFR2 does. For best accuracy in molecules and nonmetallic condensed matter, we continue to recommend SCAN and r 2 SCAN. Numerical performance is discussed in detail, and our paper provides an outlook to machine learning.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Topological Langmuir-cyclotron wave

A theory is developed to describe the topological Langmuir-cyclotron wave (TLCW), a topological excitation in magnetized plasmas recently identified by numerical simulations. As a topological wave in a continuous medium, the TLCW propagates unidirectionally without scattering in complex boundaries and can be explored as an effective mechanism to energize particles. We show that, because momentum space in continuous media is contractible in general, the topology of the wave bundles is trivial over momentum space that contains no degeneracy points. This is in stark contrast to condensed matters with periodic lattice structures that impose nontrivial topology on momentum space. In continuous media without lattice structures, nontrivial topology of the eigenmode bundles manifests over phase space, and it is the nontrivial topology over phase space that underpins the topological excitations, such as the TLCW. It is shown that the TLCW can be faithfully modeled by a generic tilted Dirac cone in phase space, whose entire spectrum, including the spectral flow, is given.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Design of quantum optical experiments with logic artificial intelligence

Logic Artificial Intelligence (AI) is a subfield of AI where variables can take two defined arguments, True or False, and are arranged in clauses that follow the rules of formal logic. Several problems that span from physical systems to mathematical conjectures can be encoded into these clauses and solved by checking their satisfiability (SAT). In contrast to machine learning approaches where the results can be approximations or local minima, Logic AI delivers formal and mathematically exact solutions to those problems. In this work, we propose the use of logic AI for the design of optical quantum experiments. We show how to map into a SAT problem the experimental preparation of an arbitrary quantum state and propose a logic-based algorithm, called Klaus, to find an interpretable representation of the photonic setup that generates it. We compare the performance of Klaus with the state-of-the-art algorithm for this purpose based on continuous optimization. We also combine both logic and numeric strategies to find that the use of logic AI significantly improves the resolution of this problem, paving the path to developing more formal-based approaches in the context of quantum physics experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Experimental Investigation of Low-Frequency Distributed Acoustic Sensor Responses to Two Parallel Propagating Fractures

Low-frequency distributed acoustic sensing (LF-DAS) is a diagnostic tool for hydraulic fracture propagation with far-field monitoring using fiber optic sensors. LF-DAS senses strain rate variation caused by stress field change due to fracture propagation. Fiber optic sensors are installed in the monitoring wells in the vicinity of a fractured well. From the strain responses, fracture propagation can be evaluated. To understand subsurface conditions with multiple propagating fractures, a laboratory-scale hydraulic fracture experiment was performed simulating the LF-DAS response to fracture propagation with embedded distributed optical fiber strain sensors under these conditions. The experiment was performed using a transparent cube of epoxy with two parallel radial initial flaws centered in the cube. Fluid was injected into the sample to generate fractures along the initial flaws. The experiment used distributed high-definition fiber optic strain sensors with tight spatial resolutions. The sensors were embedded at two different locations on opposite sides of the initial flaws, serving as observation/monitoring locations. We also employed finite element modeling to numerically solve the linear elastic equations of equilibrium continuity and stress–strain relationships. The measured strains from the experiment were compared to simulation results from the finite element model. The experimentally derived strain and strain-rate waterfall plots from this study show the responses to both fractures propagating, while the fracture at the lower position took most of the fluid during the experiment. Interestingly, a fracture first began propagating from the upper flaw of the two flaws, but once the lower fracture was initiated, it grew much faster than the upper fracture. Both fibers were intercepted by the lower fracture, further verifying the strain signature as a fracture is approaching and intersecting an offset fiber.

Chemistry↗
DOE OSTI2020

Microwave Radiometer: Meteorological Temperature Profiler (raw)

The microwave radiometer profiler (MWRP) provides vertical profiles of temperature, humidity, and cloud liquid water content as a function of height or pressure at approximately 5-minute intervals for nearly all weather conditions. The profiles are derived from measurements of absolute microwave radiances (expressed as "brightness temperatures") obtained at 12 frequencies in the range of 22-30 GHz and 51-59 GHz. The data are useful for input to numerical weather forecast models and others that require continuous, high-temporal-resolution profiles. The microwave radiance measurements are useful for testing models of microwave radiation transfer used to derive the profiles.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

High order interpolation of magnetic fields with vector potential reconstruction for particle simulations

We propose a method for interpolating divergence-free continuous magnetic fields via vector potential reconstruction using Hermite interpolation, which ensures high-order continuity for applications requiring adaptive, high-order ordinary differential equation (ODE) integrators, such as the Dormand-Prince method. The method provides C(m) continuity and achieves high-order accuracy, making it particularly suited for particle trajectory integration and Poincaré section analysis under optimal integration order and timestep adjustments. Through numerical experiments, we demonstrate that the Hermite interpolation method preserves volume and continuity, which are critical for conserving toroidal canonical momentum and magnetic moment in guiding center simulations, especially over long-term trajectory integration. Furthermore, we analyze the impact of insufficient derivative continuity on Runge-Kutta schemes and show how it degrades accuracy at low error tolerances, introducing discontinuity-induced truncation errors. Lastly, we demonstrate performant Poincaré section analysis in two relevant settings of field data collocated from finite element meshes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Exact quantization and analytic continuation

In this paper we give a streamlined derivation of the exact quantization condition (EQC) on the quantum periods of the Schrödinger problem in one dimension with a general polynomial potential, based on Wronskian relations. We further generalize the EQC to potentials with a regular singularity, describing spherical symmetric quantum mechanical systems in a given angular momentum sector. We show that the thermodynamic Bethe ansatz (TBA) equations that govern the quantum periods undergo nontrivial monodromies as the angular momentum is analytically continued between integer values in the complex plane. The TBA equations together with the EQC are checked numerically against Hamiltonian truncation at real angular momenta and couplings, and are used to explore the analytic continuation of the spectrum on the complex angular momentum plane in examples.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Automatic speech recognition predicts contemporaneous earthquake fault displacement

Abstract Significant progress has been made in probing the state of an earthquake fault by applying machine learning to continuous seismic waveforms. The breakthroughs were originally obtained from laboratory shear experiments and numerical simulations of fault shear, then successfully extended to slow-slipping faults. Here we apply the Wav2Vec-2.0 self-supervised framework for automatic speech recognition to continuous seismic signals emanating from a sequence of moderate magnitude earthquakes during the 2018 caldera collapse at the Kīlauea volcano on the island of Hawai’i. We pre-train the Wav2Vec-2.0 model using caldera seismic waveforms and augment the model architecture to predict contemporaneous surface displacement during the caldera collapse sequence, a proxy for fault displacement. We find the model displacement predictions to be excellent. The model is adapted for near-future prediction information and found hints of prediction capability, but the results are not robust. The results demonstrate that earthquake faults emit seismic signatures in a similar manner to laboratory and numerical simulation faults, and artificial intelligence models developed for encoding audio of speech may have important applications in studying active fault zones.

58 GEOSCIENCES↗
DOE OSTI2023

A general non-Fourier Stefan problem formulation that accounts for memory effects

The Stefan problem is the classical model of a melting phase change. In heterogeneous systems, such phase changes can exhibit non-Fourier (anomalous) behaviors, where the advance of the melt interface does not follow the expected time scaling. These situations can be modeled by replacing the derivatives, in the governing partial differential equations, with fractional order derivatives. In particular, replacing the time derivatives leads to non-Fourier models that account for memory effects in the system. In this work, by using appropriate time convolution integrals, a general thermodynamic balance statement for melting phase problems, explicitly accounting for memory effects, is developed. From this balance, a general model formulation applicable to problems involving melting over a temperature range (i.e., a mushy region) is derived. A key component in this model is the representation of memory effects through the use of fractional derivative based constitutive models of the enthalpy and heat flux. Further, on shrinking the mushy region to a single isotherm, a general sharp interface melting model is obtained. Here, in contrast to the classic Stefan problem, the fractional derivatives induce a natural regularization, such that the constitutive models for enthalpy and heat flux are continuous at the melt interface; a result confirmed through numerical simulation. To further support the theoretical findings, a physical example of a non-Fourier Stefan problem is presented. Overall the development and results in this paper underscore the importance of explicitly relating the development of fractional calculus models to the appropriate thermodynamic balance statements.

42 ENGINEERING↗
DOE OSTI2020

Model predictive control for voltage restoration in microgrids using temporal logic specifications

Power system operation will encounter numerous voltage variabilities as the proliferation of renewable energy continues. Real-time monitoring and communication technologies can potentially improve voltage stability by enabling the rapid detection of voltage deviations and the implementation of corrective actions. These corrective actions will only be effective in restoring stability if they are chosen in a timely and scalable manner. This paper considers the problem of power systems’ load voltage control in order to simultaneously address both magnitude and time voltage specifications. In order to comply with grid codes and avoid unnecessary relay protection actions, a model predictive control (MPC)-based control strategy employing temporal logic specifications (TLSs) is proposed. The TLSs strategy is introduced as a formalism to control the voltage variation at a critical load bus against operational bounds over time. The control objective is to schedule optimal control input signals from the available supportive energy storage systems, which provide reactive power injections, leading to satisfying the specified finite-time restoration described by the TLSs with minimal control efforts. The simulation results display that the TLSs strategy for power systems’ voltage control synthesis is extremely significant.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Circuit Quantum Electrodynamics in Hyperbolic Space: From Photon Bound States to Frustrated Spin Models

Circuit quantum electrodynamics is one of the most promising platforms for efficient quantum simulation and computation. In recent groundbreaking experiments, the immense flexibility of superconducting microwave resonators was utilized to realize hyperbolic lattices that emulate quantum physics in negatively curved space. Here we investigate experimentally feasible settings in which a few superconducting qubits are coupled to a bath of photons evolving on the hyperbolic lattice. We compare our numerical results for finite lattices with analytical results for continuous hyperbolic space on the Poincaré disk. We find good agreement between the two descriptions in the long-wavelength regime. We show that photon-qubit bound states have a curvature-limited size. We propose to use a qubit as a local probe of the hyperbolic bath, for example, by measuring the relaxation dynamics of the qubit. We find that, although the boundary effects strongly impact the photonic density of states, the spectral density is well described by the continuum theory. We show that interactions between qubits are mediated by photons propagating along geodesics. We demonstrate that the photonic bath can give rise to geometrically frustrated hyperbolic quantum spin models with finite-range or exponentially decaying interaction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Deep Reinforcement Learning for Distribution System Cyber Attack Defense with DERs

The use of smart inverter capabilities of distributed energy resources (DERs) enhances the grid reliability but in the meanwhile exhibits more vulnerabilities to cyber-attacks. This paper proposes a deep reinforcement learning (DRL)-based defense approach. The defense problem is reformulated as a Markov decision making process to control DERs and minimizing load shedding to address the voltage violations caused by cyber-attacks. The original soft actor-critic (SAC) method for continuous actions has been extended to handle discrete and continuous actions for controlling DERs' setpoints and loadshedding scenarios. Numerical comparison results with other control approaches, such as Volt-VAR and Volt-Watt on the modified IEEE 33-node, show that the proposed method can achieve better voltage regulation and have less power losses in the presence of cyber-attacks.

active distribution systems↗
DOE OSTI2021

BETO 2021 Peer Review - Algal Biofuels Techno-Economic Analysis 1.3.5.200

The objective of NREL's Algal Biofuel Techno-Economic Analysis (TEA) project is to provide process modeling and analysis to support Algae Program activities, utilizing TEA models to relate key process parameters with overall economics for cultivation, processing, and conversion of algal biomass to fuels and coproducts. By quantifying economic implications of key process metrics, TEA models highlight the technical requirements to achieve future program cost goals as well as enabling a means to track progress towards these goals. This project provides high impact and relevance though generation of critical cost data tied to funded research, with our analyses subsequently exercised by BETO to guide program plans, FOA priorities, and other directives. This includes costs for both algal biomass production and downstream conversion, most notably to support BETO's fuel cost targets below $2.5/GGE by 2030. To mitigate a key risk/challenge in constraining our work to academic analyses rooted only in future projections, our work also seeks to provide near-term value to today's algae industry through frequent industry engagement, while maintaining close ties with other BETO collaborators. This project has made numerous accomplishments since the 2019 peer review, including a continued focus on opportunities for value-added products, with related analyses for new pathway opportunities to achieve BETO cost goals and notable State-of-Technology (SOT) improvements over prior cost benchmarks.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Multistage Stabilized Continuation for Indirect Optimal Control of Three-Dimensional Hypersonic Trajectories

This work presents the development of a multistaged stabilized continuation for the three-dimensional unpowered hypersonic trajectory planning problem using indirect optimal control methods. The stabilized continuation method is noniterative and guaranteed to terminate within a finite number of floating point operations, thereby making it well suited for onboard autonomous implementations. Here, we present a multistage formulation of the stabilized continuation scheme that involves starting with a “loose” integration tolerance during the first stage and ramping up toward a “strict” integration tolerance through subsequent stages. An important benefit of this approach is that even when the solution to the underlying optimal control problem is numerically unstable, such as with the hypersonic vehicle footprint generation problem, the stabilized continuation algorithm is shown to be successful in finding a solution while providing some additional insights into the underlying cause of the numerical instability.

33 ADVANCED PROPULSION SYSTEMS↗