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

A numerical solver for investigating the space charge effect on the electric field in liquid argon time projection chambers

This paper reports the development of a numerical solver aimed to simulate the interaction between the space charge (i.e. ions) distribution and the electric field in liquid argon time projection chamber (LArTPC) detectors. The ion transport equation is solved by a time-accurate, cell-centered finite volume method and the electric potential equation by a continuous finite element method. The electric potential equation updates the electric field which provides the drift velocity to the ion transport equation. The ion transport equation updates the space charge density distribution which appears as the source term in the electric potential equation. The interaction between the space charge distribution and the electric field is numerically simulated within each physical time step. The convective velocity in the ion transport equation can include the background flow velocity in addition to the electric drift velocity. The numerical solver has been parallelized using the Message Passing Interface (MPI) library. Numerical tests show and verify the capability and accuracy of the current numerical solver. It is planned that the developed numerical solver, together with a Computational Fluid Dynamics (CFD) package which provides the flow velocity field, can be used to investigate the space charge effect on the electric field in large-scale particle detectors.

42 ENGINEERING↗
DOE OSTI2022

Electronic energy loss and ion velocity correlation effects in track production in swift-ion-irradiated LiNbO3: A quantitative assessment between structural damage morphology and energy deposition

The primary motivation for studying how irradiation modifies the structures and properties of solid materials involves the understanding of undesirable phenomena, including irradiation-induced degradation of components in nuclear reactors and space exploration, and beneficial applications, including material performance tailoring through ion beam modification and defect engineering. In this work, the formation mechanism of latent tracks with different damage morphologies in LiNbO 3 crystals under 0.09–6.17 MeV/u ion irradiation with an electronic energy loss from 2.6–13.2 keV/nm is analyzed by experimental characterizations and numerical calculations. Irradiation-induced damage is preliminarily evaluated via the prism coupling technique to analyze the correlation between the dark-mode spectra and energy loss profiles of irradiated regions. Under the irradiation conditions of different ion velocities and electronic energy losses, different damage morphologies, from individual spherical defects to discontinuous and continuous tracks, are experimentally characterized. During ion penetration process, the ion velocity determines the spatiotemporal distribution of deposited irradiation energy induced by electronic energy loss, meaning that the two essential factors including electronic energy loss and ion velocity co-affect the track damage. The inelastic thermal spike model is used to numerically calculate the spatiotemporal evolutions of energy deposition and the corresponding atomic temperature under different irradiation conditions, and a quantitative relationship is proposed by comparison with corresponding experimentally observed track damage morphologies. Additionally, the obtained quantitative relationship between irradiation conditions and track damage provides deep insight and guidance for understanding the damage behavior of crystal materials in extreme radiation environments and selecting irradiation parameters, including ion species and energies, for ion beam technique application in atomic-level defect manipulation, material modification, and micro/nanofabrication.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Extreme Scale Infrasound Inversion and Prediction for Weather Characterization and Acute Event Detection

Accurate and timely weather predictions are critical to many aspects of society with a profound impact on our economy, general well-being, and national security. In particular, our ability to forecast severe weather systems is necessary to avoid injuries and fatalities, but also important to minimize infrastructure damage and maximize mitigation strategies. The weather community has developed a range of sophisticated numerical models that are executed at various spatial and temporal scales in an attempt to issue global, regional, and local forecasts in pseudo real time. The accuracy however depends on the time period of the forecast, the nonlinearities of the dynamics, and the target spatial resolution. Significant uncertainties plague these predictions including errors in initial conditions, material properties, data, and model approximations. To address these shortcomings, a continuous data collection occurs at an effort level that is even larger than the modeling process. It has been demonstrated that the accuracy of the predictions depends on the quality of the data and is independent to a certain extent on the sophistication of the numerical models. Data assimilation has become one of the more critical steps in the overall weather prediction business and consequently substantial improvements in the quality of the data would have transformational benefits. This paper describes the use of infrasound inversion technology, enabled through exascale computing, that could potentially achieve orders of magnitude improvement in data quality and therefore transform weather predictions with significant impact on many aspects of our society.

58 GEOSCIENCES↗
DOE OSTI2020

Printing 2D Conjugated Polymer Monolayers and Their Distinct Electronic Properties

Abstract Recently, 2D monolayer films of conjugated polymers have gained increasing attention owing to the preeminence of 2D inorganic films that exhibit unique optoelectronic and mechanical properties compared to their bulk analogs. Despite numerous efforts, crystallization of semiconducting polymers into highly ordered 2D monolayer films still remains challenging. Herein, a dynamic‐template‐assisted meniscus‐guided coating is utilized to fabricate continuous, highly ordered 2D monolayer films of conjugated polymers over a centimeter scale with enhanced backbone π–π stacking. In contrast, monolayer films printed on solid substrates confer upon the 1D fiber networks strong alkyl side‐chain stacking at the expense of backbone packing. From single‐layers to multilayers, the polymer π‐stacks change from edge‐on to bimodal orientation as the film thickness reaches ≈20 nm. Spectroscopic and cyclic voltammetry analysis reveals an abrupt increase in J‐aggregation and absorption coefficient and a decrease in bandgap and highest occupied molecular orbital level until critical thickness, possibly arising from the straightened polymer backbone. This is corroborated by an abrupt increase in hole mobility with film thickness, reaching a maximum of 0.7 cm 2 V −1 s −1 near the critical thickness. Finally, fabrication of chemical sensors incorporating polymer films of various thicknesses is demonstrated, and an ultrahigh sensitivity of the ≈7 nm thick ultrathin film (bilayers) to 1 ppb ammonia is shown.

Kafle, Prapti↗
DOE OSTI2024

Observation and Modeling of Dynamic Fracture Behaviors of Battery Cell Under Impact Loading Using Enhanced Representative Volume Element Concept

The burgeoning electric automobile industry has increased interest in battery safety. Battery cells experience significant mechanical stress during operation, including the impact of accidents and vibrations from driving. The potential for thermal runaway reactions in battery cells raises safety concerns. Although numerous researchers have defined the dynamic behavior of battery cells and proposed numerical models to describe it, few studies have focused on the high-strain rate mechanical impact phase correlated with the onset of fracture. In this study, we describe the dynamic behavior of pouch battery cells and propose a modeling method to study their mechanical failure under impact situations. Impact tests are conducted at various velocities and heights. To overcome numerical issues commonly encountered under rapid deformation scenarios, a new finite element model is developed based on the representative volume element model. The proposed approach efficiently simulates continuous crack propagation and brittleness behavior during impact by permitting the individual behavior of the cell components. Therefore, engineers can reliably design safer electric vehicle battery cells by measuring the properties of the cell components.

ENERGY STORAGE↗
DOE OSTI2022

Neuromorphic Computing for Scientific Applications

Neuromorphic computing technology continues to make strides in the development of new algorithms, devices, and materials. In addition, applications have begun to emerge where neuromorphic computing shows promising results. However, numerous barriers to further development and application remain. In this work, we identify several science areas where neuromorphic computing can either make an immediate impact (within 1 to 3 years) or the societal impact would be extremely high if the technological barriers can be addressed. We identify both opportunities and hurdles to the development of neuromorphic computing technology for these areas. Finally, we discuss future directions that need to be addressed to expand both the development and application of neuromorphic computing.

Patton, Robert↗
DOE OSTI2022

Machine Learning-Assisted High-Temperature Reservoir Thermal Energy Storage Optimization: Numerical Modeling and Machine Learning Input and Output Files

This data set includes the numerical modeling input files and output files used to synthesize data, and the reduced-order machine learning models trained from the synthesized data for reservoir thermal energy storage site identification. In this study, a machine-learning-assisted computational framework is presented to identify High-Temperature Reservoir Thermal Energy Storage (HT-RTES) site with optimal performance metrics by combining physics-based simulation with stochastic hydrogeologic formation and thermal energy storage operation parameters, artificial neural network regression of the simulation data, and genetic algorithm-enabled multi-objective optimization. A doublet well configuration with a layered (aquitard-aquifer-aquitard) generic reservoir is simulated for cases of continuous operation and seasonal-cycle operation scenarios. Neural network-based surrogate models are developed for the two scenarios and applied to generate the Pareto fronts of the HT-RTES performance for four potential HT-RTES sites. The developed Pareto optimal solutions indicate the performance of HT-RTES is operation-scenario (i.e., fluid cycle) and reservoir-site dependent, and the performance metrics have competing effects for a given site and a given fluid cycle. The developed neural network models can be applied to identify suitable sites for HT-RTES, and the proposed framework sheds light on the design of resilient HT-RTES systems. All the simulations and the neural network model were done by Idaho National Laboratory. A detailed description of the work was reported in publication linked below.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Tomographic X-ray particle tracking velocimetry

We investigate the feasibility of in-laboratory tomographic X-ray particle tracking velocimetry (TXPTV) and consider creeping flows with nearly density matched flow tracers. Specifically, in these proof-of-concept experiments we examined a Poiseuille flow, flow through porous media and a multiphase flow with a Taylor bubble. For a full 360° computed tomography (CT) scan we show that the specially selected 60 micron tracer particles could be imaged in less than 3 seconds with a signal-to-noise ratio between the tracers and the fluid of 2.5, sufficient to achieve proper volumetric segmentation at each time step. In the pipe flow, continuous Lagrangian particle trajectories were obtained, after which all the standard techniques used for PTV or PIV (taken at visible wave lengths) could also be employed for TXPTV data. And, with TXPTV we can examine flows inaccessible with visible wave lengths due to opaque media or numerous refractive interfaces. In the case of opaque porous media we were able to observe material accumulation and pore clogging, and for flow with Taylor bubble we can trace the particles and hence obtain velocities in the liquid film between the wall and bubble, with thickness of liquid film itself also simultaneously obtained from the volumetric reconstruction after segmentation. While improvements in scan speed are anticipated due to continuing improvements in CT system components, we show that for the flows examined even the presently available CT systems could yield quantitative flow data with the primary limitation being the quality of available flow tracers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

On the Equilibrium Electrostatic Potential and Light-Induced Charge Redistribution in Halide Perovskite Structures

Lead halide perovskites are semiconductor materials which are employed as nonintentionally doped absorbers inserted between two selective carrier transport layers (SCTL), realizing a p-i-n or n-i-p heterojunction. In our study, we have developed and investigated a lateral device, based on methylammonium lead iodide (MAPbI3) in which the p-i-n heterojunction develops in the horizontal direction. Our research suggests that the effective doping level in the MAPbI3 film should be very low, below 1012 cm-3. Along the vertical direction, this doping level is not enough to screen the electric field of the buried heterojunction with the SCTL. The perovskite work function is therefore affected by the work function of the SCTL underneath. From drift-diffusion simulations, we show that intrinsic perovskite-SCTL structures develop mV range surface photovoltages (SPVs) under continuous illumination. However, perovskite-SCTL structures can develop SPVs of hundreds of mV, as confirmed by our measurements. We therefore analyzed the compatibility between low doping and low defect densities in the perovskite layer and such high SPV values using numerical modeling. It is shown that these high SPV values could originate from electronic processes due to large band offsets in the buried perovskite-SCTL heterojunctions, or at the SCTL-transparent conductive oxide (TCO) buried heterojunction. However, such electronic processes can hardly explain the long SPV persistence after switching off the illumination.

defects↗
DOE OSTI2021

Coupling Surface Flow with High-performance Subsurface Reactive Flow and Transport Code PFLOTRAN

Water exchange between the surface and subsurface is important for both water resource management and environmental protection. In this paper, we develop coupled surface and subsurface flow simulation capability in a parallel subsurface flow and reactive transport code PFLOTRAN. We sequentially couple the diffusion wave-based surface flow with the subsurface flow governedby the Richards equation in PFLOTRAN. These two flow domains are linked with a boundary condition switching method that ensures continuity of pressure and flux at the surface-subsurface interface. We verify the coupled code against other existing hydrologic models and observation data using a number of numerical experiments. The coupled hydrological model exhibits good performance in strong parallel scaling tests. The new coupled surface and subsurface simulator significantly advance community simulation capability towards improving integrated hydrologic and biogeochemical understanding of complex systems such as watersheds and river corridors. Keywords: Surface flow, Integrated hydrological modeling, Boundary condition switching, Parallel computing

Wu, Runjian↗
DOE OSTI2026

Magnon-Magnon Interaction Induced by Dynamic Coupling in a Hybrid Magnonic Crystal

We report a combined experimental and numerical investigation of spin-wave dynamics in a hybrid magnonic crystal consisting of a CoFeB artificial spin ice (ASI) of stadium-shaped nanoelements patterned atop a continuous NiFe film separated by a 5 nm Al 2 O 3 spacer. Using Brillouin light scattering spectroscopy, we probe the frequency dependence of thermal spin waves as functions of applied magnetic field and wavevector, revealing the decisive role of interlayer dipolar coupling in the magnetization dynamics. Micromagnetic simulations complement the experiments, showing a strong interplay between ASI edge modes and backward volume modes in the NiFe film. The contrast in saturation magnetization between CoFeB and NiFe enhances this coupling, leading to a pronounced hybridization manifested as a triplet of peaks in the BLS spectra predicted by simulations and observed experimentally. This magnon−magnon coupling persists over a wide magnetic field range, shaping both the spin-wave dispersion at fixed fields and the full frequency-field response throughout the magnetic hysteresis loop. Our findings establish how ASI geometry can selectively enhance specific spin-wave wavelengths in the underlying film, thereby boosting their amplitude and identifying them as preferential channels for spin wave transmission and manipulation.

Brillouin light scattering spectroscopy↗
DOE OSTI2026

Benchmark for two-dimensional large scale coherent structures in partially magnetized E × B plasmas—community collaboration & lessons learned

Low-temperature plasmas (LTPs) are essential to both fundamental scientific research and critical industrial applications. As in many areas of science, numerical simulations have become a vital tool for uncovering new physical phenomena and guiding technological development. Code benchmarking remains crucial for verifying implementations and evaluating performance. This work continues the Landmark benchmark initiative, a series specifically designed to support the verification of LTP codes. In this study, seventeen simulation codes from a collaborative community of nineteen international institutions modeled a partially magnetized E × B Penning discharge. The emergence of large scale coherent structures, or rotating plasma spokes, endows this configuration with an enormous range of time scales, making it particularly challenging to simulate. The codes showed excellent agreement on the rotation frequency of the spoke as well as key plasma properties, including time-averaged ion density, plasma potential, and electron temperature profiles. Achieving this level of agreement came with challenges, and we share lessons learned on how to conduct future benchmarking campaigns. Comparing code implementations, computational hardware, and simulation runtimes also revealed interesting trends, which are summarized with the aim of guiding future plasma simulation software development.

benchmarking↗
DOE OSTI2022

Modeling and Market Design Considerations for Conventional and Decarbonized Resources

The following paper reviews various market design issues that will need to be reconsidered due to anticipated changes in the resources that supply energy in wholesale electricity markets. As the energy supply continues towards decarbonization, the change in resource technologies will affect the fundamentals of production scheduling and the policies to address supply variability and uncertainty. We show through simple numerical examples that the existing market design approach based on fuel costs will result in $0/MWh prices with intermittent price spikes during reserve shortages, but that incorporating more granular reserve pricing, energy storage participation, and price-responsive demand participation can restore efficient market clearing with reasonable pricing outcomes. The market design approach fundamentally shifts from fuel based to opportunity cost based. We briefly review alternative market design frameworks that can support this shift, including intraday markets, decentralized markets, flexibility options, and swing contracts. Market designs may also be required to accommodate or support various out-of-market policies and agreements; ideally, these external factors can be integrated into the market design to facilitate efficient exchanges across longer time scales, between other markets, and in support of public policy goals. The paper concludes by discussing how decarbonization trends may affect the design and use of production cost models for short term operations and long term planning studies.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Electronic Interactions Between the Receptor-Binding Domain of Omicron Variants and Angiotensin-Converting Enzyme 2: A Novel Amino Acid–Amino Acid Bond Pair Concept

SARS-CoV-2 remains a severe threat to worldwide public health, particularly as the virus continues to evolve and diversify into variants of concern (VOCs). Among these VOCs, Omicron variants exhibit unique phenotypic traits, such as immune evasion, transmissibility, and severity, due to numerous spike protein mutations and the rapid subvariant evolution. These Omicron subvariants have more than 15 mutations in the receptor-binding domain (RBD), a region of the SARS-CoV-2 spike protein that is important for recognition and binding with the angiotensin-converting enzyme 2 (ACE2) human receptor. To address the impact of these high numbers of Omicron mutations on the binding process, we have developed a novel method to precisely quantify amino acid interactions via the amino acid–amino acid bond pair (AABP). We applied this concept to investigate the interface interactions of the RBD–ACE2 complex in four Omicron Variants (BA.1, BA.2, BA.5, and XBB.1.16) with its Wild Type counterpart. Based on the AABP analysis, we have identified all the sites that are affected by mutation and have provided evidence that unmutated sites are also impacted by mutation. We have calculated that the binding between RBD and ACE2 is strongest in OV BA.1, followed by OV BA.2, WT, OV BA.5, and OV XBB.1.16. We also present the partial charge values for all 311 residues across these five models. Our analysis provides a detailed understanding of changes caused by mutation in each Omicron interface complex.

Biochemistry & Molecular Biology↗
DOE OSTI2023

An Efficient High-Order Solver for Diffusion Equations with Strong Anisotropy on Non-Anisotropy-Aligned Meshes

This paper concerns numerical solution of the diffusion equation with strong anisotropy on meshes not aligned with the anisotropic vector field. In order to resolve the numerical pollution for simulations on a non-anisotropy-aligned mesh and reduce the associated high computational cost we propose an effective preconditioner, extending our previous work. Similar to the anisotropy-aligned mesh case, we apply the auxiliary space preconditioning framework to design a preconditioner where a continuous finite element space is used as the auxiliary space for the discontinuous finite element space. The key component is an effective line smoother that can mitigate the high-frequency errors perpendicular to the magnetic field. We design a graph-based approach to find such a line smoother that is approximately perpendicular to the vector fields when the mesh does not align with the anisotropy. Finally, numerical experiments for several benchmark problems are presented, demonstrating the effectiveness and robustness of the proposed preconditioner when applied to Krylov iterative methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A least-squares finite element method based on the Helmholtz decomposition for hyperbolic balance laws

In this paper, a least-squares finite element method for scalar nonlinear hyperbolic balance laws is proposed and studied. The approach is based on a formulation that utilizes an appropriate Helmholtz decomposition of the flux vector and is related to the standard notion of a weak solution. This relationship, together with a corresponding connection to negative-norm least-squares, is described in detail. As a consequence, an important numerical conservation theorem is obtained, similar to the famous Lax–Wendroff theorem. The numerical conservation properties of the method in this paper do not fall precisely in the framework introduced by Lax and Wendroff, but they are similar in spirit as they guarantee that when L 2 convergence holds, the resulting approximations approach a weak solution to the hyperbolic problem. The least-squares functional is continuous and coercive in an H -1 -type norm, but not L 2 -coercive. Nevertheless, the L 2 convergence properties of the method are discussed. Convergence can be obtained either by an explicit regularization of the functional, that provides control of the L 2 norm, or by properly choosing the finite element spaces, providing implicit control of the L 2 norm. Numerical results for the inviscid Burgers equation with discontinuous source terms are shown, demonstrating the L 2 convergence of the obtained approximations to the physically admissible solution. The numerical method utilizes a least-squares functional, minimized on finite element spaces, and a Gauss–Newton technique with nested iteration. Finally, we believe that the linear systems encountered with this formulation are amenable to multigrid techniques and combining the method with adaptive mesh refinement would make this approach an efficient tool for solving balance laws (this is the focus of a future study).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Repurposing Inactive Oil and Gas Wells for Energy Storage: Maximizing the Potential Via Optimal Drivetrain Control: Preprint

In recent years, there has been a growing emphasis on utilizing energy storage to enhance grid resilience against disruptive events. While renewable energy supply continues to expand, gravity-based solutions like pumped hydro remain dominant in the commercial space. However, their geographical limitations constrain availability, scalability, and increase costs for co-locating solar and wind energy. An alternative approach proposes repurposing idle oil and gas wells located closer to existing grid infrastructure, offering a promising and cost-effective solution. This paper addresses the optimization and control of a regenerative drive system coupled to an interior permanent magnet synchronous machine in a 300-meter well with a 100-Newton weight. The study employs a dynamic MATLAB/Simulink model to simulate the operation of the electric drivetrain system during storage and discharge operations. The results demonstrate an initial round-trip efficiency of 85.9% for the electrical system alone and identify crucial factors for maximizing efficiency. The optimized operation and control of the electromechanical drivetrain system hold great potential for minimizing the levelized cost of storage while maximizing efficiency and revenue generation.

ENERGY STORAGE,ENGINEERING↗
DOE OSTI2023

Learning subgrid-scale models with neural ordinary differential equations

We propose a new approach to learning the subgrid-scale model when simulating partial differential equations (PDEs) solved by the method of lines and their representation in chaotic ordinary differential equations, based on neural ordinary differential equations (NODEs). Solving systems with fine temporal and spatial grid scales is an ongoing computational challenge, and closure models are generally difficult to tune. Machine learning approaches have increased the accuracy and efficiency of computational fluid dynamics solvers. In this approach neural networks are used to learn the coarse- to fine-grid map, which can be viewed as subgrid-scale parameterization. We propose a strategy that uses the NODE and partial knowledge to learn the source dynamics at a continuous level. Our method inherits the advantages of NODEs and can be used to parameterize subgrid scales, approximate coupling operators, and improve the efficiency of low-order solvers. Numerical results with the two-scale Lorenz 96 ODE, the convection-diffusion PDE, and the viscous Burgers' PDE are used to illustrate this approach.

97 MATHEMATICS AND COMPUTING↗