Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mechanical robustness”

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 1,225 records · Page 68

Multi-physics Modeling of Radiative Heat Transfer and Fluid Flow for the Reactor Cavity Cooling System

High-temperature gas-cooled reactors (HTGRs) are notable for their high thermal efficiency and potential for combined heat and power applications. These reactors are particularly appealing due to their advanced passive safety features. HTGRs utilize passive safety systems that function without requiring active components like pumps or compressors during emergencies. These reactor designs depend on a Reactor Cavity Cooling System (RCCS) to manage decay heat removal from the reactor pressure vessel (RPV) during accident conditions. The RCCS consists of vertical rectangular channels known as "risers" or riser ducts positioned around the RPV. These risers receive heat from the RPV through both convective and radiative heat transfer mechanisms. Understanding the interplay of multiple physical phenomena, such as fluid dynamics, heat transfer, and neutron interactions, is essential for the effective design and operation of nuclear reactors, particularly for systems like the RCCS. Multi-physics simulations provide a comprehensive approach to studying these interactions, offering detailed insights and enhancing accuracy. They are especially important in RCCS designs, where the interaction between radiative and convective heat transfer can significantly impact system performance. By leveraging multi-physics simulations, complex reactor behaviors can be modeled without compromising the fidelity of the underlying physical processes. This work aims to establish a robust methodology for coupling multiple physical processes in an air-cooled RCCS. By focusing on validating this multi-physics approach, the study involves designing test cases that simulate various conditions to verify the numerical models employed. The outcomes of this research will provide critical insights for accurately modeling and optimizing complex nuclear systems like the RCCS.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Consistency-Enhanced Evolution for Variable Selection Can Identify Key Chemical Information from Spectroscopic Data

In the last few decades, spectroscopic techniques such as near-infrared (NIR) spectroscopy have gained wide applications in several industries, such as the pharmaceutical, agricultural, oil, and gas industries. As a result, various soft sensors have been developed to predict sample properties from spectroscopic readings. Because the spectroscopic readings at different wavelengths, especially at the adjacent wavelengths, are highly correlated, it has been shown that variable selection could significantly improve a soft sensor’s prediction performance while reducing the model complexity. To improve the prediction performance, most variable selection methods focus on identifying the variables (i.e., wavelengths or wavelength segments) that are strongly correlated with the dependent variable. Although many successful applications have been reported, these variable selection methods do have their limitations. Specifically, the selected wavelengths sometimes show little connection to the chemical bounds or functional groups presenting in the sample. In addition, the selected variables can be quite sensitive to the choice of the training samples. In this work, we address these limitations from a different perspective: if a variable selection algorithm can identify the truly relevant input variables, it should consistently identify the same subset of variables regardless of the choice of the training samples. Therefore, we propose a variable selection method that aims to improve the consistency of variable selection resulting from different training samples. Furthermore, the new algorithm is termed consistency-enhanced evolution for variable selection (CEEVS). To demonstrate the performance and robustness of CEEVS, we compare the proposed method with three representative variable selection methods using five published NIR data sets. These case studies clearly demonstrate that by improving the variable selection consistency, we can not only achieve improved prediction performance, but also identify key chemical information from spectroscopic data.

42 ENGINEERING↗

Detection of Outliers in LiDAR Data Acquired by Multiple Platforms over Sorghum and Maize

High-resolution point cloud data acquired with a laser scanner from any platform contain random noise and outliers. Therefore, outlier detection in LiDAR data is often necessary prior to analysis. Applications in agriculture are particularly challenging, as there is typically no prior knowledge of the statistical distribution of points, plant complexity, and local point densities, which are crop-dependent. The goals of this study were first to investigate approaches to minimize the impact of outliers on LiDAR acquired over agricultural row crops, and specifically for sorghum and maize breeding experiments, by an unmanned aerial vehicle (UAV) and a wheel-based ground platform; second, to evaluate the impact of existing outliers in the datasets on leaf area index (LAI) prediction using LiDAR data. Two methods were investigated to detect and remove the outliers from the plant datasets. The first was based on surface fitting to noisy point cloud data via normal and curvature estimation in a local neighborhood. The second utilized the PointCleanNet deep learning framework. Both methods were applied to individual plants and field-based datasets. To evaluate the method, an F-score was calculated for synthetic data in the controlled conditions, and LAI, the variable being predicted, was computed both before and after outlier removal for both scenarios. Results indicate that the deep learning method for outlier detection is more robust than the geometric approach to changes in point densities, level of noise, and shapes. The prediction of LAI was also improved for the wheel-based vehicle data based on the coefficient of determination (R2) and the root mean squared error (RMSE) of the residuals before and after the removal of outliers.

36 MATERIALS SCIENCE↗

Performance measurement of emerging 3- and 4-terminal tandem solar cells

Tandem solar cells are not limited to the conventional two-terminal (2-T) configuration. Multi-terminal designs like three-terminal (3-T) and four-terminal (4-T) devices have gained increasing attention in the PV community due to their relaxed current-matching requirement between subcells and their potential for enhanced energy yield. However, reliable and standardized methods for evaluating the performance of multi-terminal tandems remain underdeveloped. This work addresses this gap by providing comprehensive measurement guidelines tailored to these advanced configurations. We examine key coupling mechanisms between subcells, including the shared electrical load in 3-T devices and optical luminescent coupling in both 3-T and 4-T devices, to enable accurate and consistent performance evaluation. Furthermore, we propose two stabilized measurement methods for emerging 3-T tandem cells incorporating perovskite subcells: (1) a two-dimensional maximum-power-point tracking (MPPT) approach that continuously tracks both subcells' maximum power points (P MAX ) until convergence to stabilized outputs, and (2) a hybrid approach that combines MPPT for one subcell with stabilized current recording under fixed voltage biases near the P MAX of the other, allowing robust extraction of the overall stabilized P MAX (termed “MPPT + asymptotic P MA X scan” method). These methods directly address the dynamic current responses inherent to perovskite-containing tandems, providing a foundation for meaningful and reproducible performance comparisons.

14 SOLAR ENERGY↗

Structure and Function of REP34 Implicates Carboxypeptidase Activity in Francisella tularensis Host Cell Invasion

Francisella tularensis is the etiological agent of tularemia, or rabbit fever. Although F. tularensis is a recognized biothreat agent with broad and expanding geographical range, its mechanism of infection and environmental persistence remain poorly understood. Previously, we identified seven F. tularensis proteins that induce a rapid encystment phenotype (REP) in the free-living amoeba, Acanthamoeba castellanii. Encystment is essential to the pathogen's long term intracellular survival in the amoeba. Here, we characterize the cellular and molecular function of REP34, a REP protein with a mass of 34 kDa. A REP34 knock-out strain of F. tularensis has a reduced ability to both induce encystment in A. castellanii and invade human macrophages. We determined the crystal structure of REP34 to 2.05-Å resolution and demonstrate robust carboxypeptidase B-like activity for the enzyme. REP34 is a zinc-containing monomeric protein with close structural homology to the metallocarboxypeptidase family of peptidases. REP34 possesses a novel topology and substrate binding pocket that deviates from the canonical funnelin structure of carboxypeptidases, putatively resulting in a catalytic role for a conserved tyrosine and distinct S1' recognition site. Taken together, these results identify REP34 as an active carboxypeptidase, implicate the enzyme as a potential key F. tularensis effector protein, and may help elucidate a mechanistic understanding of F. tularensis infection of phagocytic cells.

59 BASIC BIOLOGICAL SCIENCES↗

Robust Molecular Predictive Methods for Novel Polymer Discovery and Applications

Polymeric materials are ubiquitous in modern society and they play an instrumental role in almost all industries, undoubtedly including the energy and environment sectors. Increased demand of energy and awareness to sustainability both necessitates the development of novel polymers with enhanced properties. Unfortunately, their structural and behavioral complexity render such discovery challenging and impeded. To address this problem, scientists are developing various computational modeling techniques and leveraging their power to depict the relationship between structural characteristics of polymers and their properties (such as rheological behaviors), and use such prediction to guide the design and syntheses of novel polymeric materials with enhanced performances. Unfortunately, predicting the relationships between polymer structure and composition with rheological properties via atomistic modeling is still a major challenge because of the extended time and length scales involved. Studying dynamic shear viscosity and linear viscoelasticity using molecular models requires capabilities that have been elusive, including representation of large molecular weight chains with an effective internal scale capable of describing entanglement, shear-rates that are in the s-1 scale with accurate quantitative stresses, and chemically-realistic combinations of both homogeneous and heterogeneous systems. Motivated by these unmet challenges, the overall technical objective of this DOE-STTR Phase II project is to develop robust molecular predictive methods for advanced polymer discovery and applications and especially for designing and demonstrating the “smart” polymer-based waterflooding enhanced oil recovery (EOR) process. In particular, we apply state-of-the-art molecular modeling methods developed by our academic partner, Materials Stimulation Center (MSC) at California Institute of Technology (Caltech), to facilitate and accelerate the experimental discovery processes. During the Phase I of this project, we had focused on development and demonstration of the molecular modeling methods to describe rheological properties of non-Newtonian polymer fluids, and to improve our fundamental understandings of shear-thickening mechanism and kinetics. In Phase II, we further apply the theoretical models to guide our experimental programs to improve our design of smart rheology modifier (SRM) polymers and their optimization for EOR. Specifically, we have three objectives in the Phase II study: (1) to further improve out computational modeling methods, coupling with the advanced machine learning algorithms; (2) to develop cost-effective and efficient SRM-flooding process suitable for EOR applications under typical reservoir conditions; and (3) to further explore the application of our molecular predictive models for innovative material discovery in other industrial applications. The recent development of our multiscale predictive framework allows the successful prediction of rheological properties from the chemical structure for polymers of experimentally relevant molecular weights, and provides an in-silico machine learning engine for screening novel compositions and structures with optimized non-Newtonian response, required for both shear-thinning and shear-thickening applications. Our framework provides: (1) procedures and tools for systematic coarsening from atomistic models and reverse mapping of coarse-grain models to atomistic, (2) unique ab initio methods to characterize the atomistic origin of colloidal and interfacial interactions and phenomena, (3) systematic structure and composition builders based on practical descriptors that drive rheological changes in polymer melts and diluted polymer mixtures, (4) a rheological properties engine capable of predicting viscosity in the zero-shear limit and under realistic dynamic conditions (for shear-rates commensurate with experiments) for large heterogeneous systems, (5) coarse-grain force fields with improved non-bond descriptions based on accurate quantum mechanics, (6) an in-silico screening machine learning engine that feeds from the systematic model builders to cover the descriptors search space, computes the rheological properties from converged trajectories spanning sub-milliseconds and ranks them for each structure/composition using an automated viscosity-vs-shear rate fitness function that can be tuned for shear-thickening, shear-thinning and other rheological responses.

02 PETROLEUM↗

Molten sodium batteries: advances in chemistries, electrolytes, and interfaces

The need for clean, renewable energy has driven the expansion of renewable energy generators, such as wind and solar. However, to achieve a robust and responsive electrical grid based on such inherently intermittent renewable energy sources, grid-scale energy storage is essential. The unmet need for this critical component has motivated extensive grid-scale battery research, especially exploring chemistries “beyond Li-ion”. Among others, molten sodium (Na) batteries, which date back to the 1960s with Na-S, have seen a strong revival, owing mostly to raw material abundance and the excellent electrochemical properties of Na metal. Recently, many groups have demonstrated important advances in battery chemistries, electrolytes, and interfaces to lower material and operating costs, enhance cyclability, and understand key mechanisms that drive failure in molten Na batteries. For widespread implementation of molten Na batteries, though, further optimization, cost reduction, and mechanistic insight is necessary. In this light, this work provides a brief history of mature molten Na technologies, a comprehensive review of recent progress, and explores possibilities for future advancements.

25 ENERGY STORAGE↗

Advancing the HERO WEC Through Integrated Modeling, Testing, and Field Deployments

The Hydraulic and Electric Reverse Osmosis Wave Energy Converter (HERO WEC) was developed by the National Laboratory of the Rockies as a modular platform for advancing wave-powered desalination technologies. Designed to operate in either a hydraulic or electric configuration, the system enables direct comparison of mechanical-to-water and electrical-to-water conversion pathways using a single hardware architecture. Deployments in 2022 and 2024 at the Jennette's Pier test site in North Carolina demonstrated freshwater production between 60 and 300 gallons per day, including successful operation in wave heights as low as 0.25 m. A structured evaluation approach combining numerical modeling of hydrodynamic and power take-off response, controlled laboratory testing, and in-ocean field deployments has been used to characterize and refine system performance. Specifically, insights from these efforts are now informing the second-generation HERO WEC (V2) to increase robustness and overall survivability. Changes from the original prototype include improved drivetrain components, refined geometry, enhanced control systems, and design. The HERO WEC platform continues to serve as an open-access research asset supporting broader marine energy and desalination development.

16 TIDAL AND WAVE POWER↗

Extended Barton–Bandis model for rock joints under cyclic loading: Formulation and implicit algorithm

In this paper, the Barton–Bandis model for rock joints is extended to cyclic loading conditions, without any new material parameter. Also developed herein is an implicit solution algorithm for the extended Barton–Bandis model, which can also be used for the original Barton–Bandis model for which an implicit algorithm has been unavailable. To this end, we first cast the Barton–Bandis model into an incremental elasto-plastic framework, deriving an expression for the elastic shear stiffness being consistent with the original model formulation. We then extend the model formulation to cyclic loading conditions, incorporating the dependence of shear stress and dilation on the joint position and the shearing direction. The extension is achieved by introducing a few state-dependent variables which can be calculated with the existing material parameters. For robust and accurate utilization of the model, we also develop an implicit algorithm based on return mapping, which is unconditionally stable and guarantees the satisfaction of the strength criterion. We verify that the proposed model formulation and algorithm produce virtually the same results as the original Barton–Bandis model under monotonic shearing conditions. Here, we then validate the extended Barton–Bandis model against experimental data on natural rock joints under cycling loading conditions. The present work thus enables the Barton–Bandis model, which has been exceptionally popular in research and practice, to be applicable to a wider range of problems in rock mechanics and rock engineering.

58 GEOSCIENCES↗

A quantitative framework reveals ecological drivers of grassland microbial community assembly in response to warming

Abstract Unraveling the drivers controlling community assembly is a central issue in ecology. Although it is generally accepted that selection, dispersal, diversification and drift are major community assembly processes, defining their relative importance is very challenging. Here, we present a framework to quantitatively infer community assembly mechanisms by phylogenetic bin-based null model analysis (iCAMP). iCAMP shows high accuracy (0.93–0.99), precision (0.80–0.94), sensitivity (0.82–0.94), and specificity (0.95–0.98) on simulated communities, which are 10–160% higher than those from the entire community-based approach. Application of iCAMP to grassland microbial communities in response to experimental warming reveals dominant roles of homogeneous selection (38%) and ‘drift’ (59%). Interestingly, warming decreases ‘drift’ over time, and enhances homogeneous selection which is primarily imposed on Bacillales. In addition, homogeneous selection has higher correlations with drought and plant productivity under warming than control. iCAMP provides an effective and robust tool to quantify microbial assembly processes, and should also be useful for plant and animal ecology.

59 BASIC BIOLOGICAL SCIENCES↗

Distributed quantum sensing enhanced by continuous-variable error correction

Abstract A distributed sensing protocol uses a network of local sensing nodes to estimate a global feature of the network, such as a weighted average of locally detectable parameters. In the noiseless case, continuous-variable (CV) multipartite entanglement shared by the nodes can improve the precision of parameter estimation relative to the precision attainable by a network without shared entanglement; for an entangled protocol, the root mean square estimation error scales like 1/ M with the number M of sensing nodes, the so-called Heisenberg scaling, while for protocols without entanglement, the error scales like 1 / M . However, in the presence of loss and other noise sources, although multipartite entanglement still has some advantages for sensing displacements and phases, the scaling of the precision with M is less favorable. In this paper, we show that using CV error correction codes can enhance the robustness of sensing protocols against imperfections and reinstate Heisenberg scaling up to moderate values of M . Furthermore, while previous distributed sensing protocols could measure only a single quadrature, we construct a protocol in which both quadratures can be sensed simultaneously. Our work demonstrates the value of CV error correction codes in realistic sensing scenarios.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Data-Driven Framework for Direct Local Tensile Property Prediction of Laser Powder Bed Fusion Parts

This article proposes a generalizable, data-driven framework for qualifying laser powder bed fusion additively manufactured parts using part-specific in situ data, including powder bed imaging, machine health sensors, and laser scan paths. To achieve part qualification without relying solely on statistical processes or feedstock control, a sequence of machine learning models was trained on 6299 tensile specimens to locally predict the tensile properties of stainless-steel parts based on fused multi-modal in situ sensor data and a priori information. A cyberphysical infrastructure enabled the robust spatial tracking of individual specimens, and computer vision techniques registered the ground truth tensile measurements to the in situ data. The co-registered 230 GB dataset used in this work has been publicly released and is available as a set of HDF5 files. The extensive training data requirements and wide range of size scales were addressed by combining deep learning, machine learning, and feature engineering algorithms in a relay. The trained models demonstrated a 61% error reduction in ultimate tensile strength predictions relative to estimates made without any in situ information. Lessons learned and potential improvements to the sensors and mechanical testing procedure are discussed.

36 MATERIALS SCIENCE↗

Hierarchical Defect Engineering for LiCoO 2 through Low-Solubility Trace Element Doping

Real-world industry-relevant battery composite electrodes are hierarchically structured. Their structural and chemical complexity is featured by ubiquitous multi-scale porosity and cracks, solid-solid and solid-liquid interfaces, compositional and redox heterogeneity, as well as lattice disordering and deformation. In particular for the active cathode particles, which are the fundamental building blocks for the energy reservoir, it is a consensus that these structural and chemical defects could have a profound impact on the battery performance. An in-depth understanding of the underlying mechanisms could critically inform the cathode material engineering, which would have a tremendous potential but remains a daunting challenge at present. In this work, we tackle this question by studying LiCoO 2 (LCO) with trace doping of Ti, which exhibits a low solubility in the LCO layered lattice. Additionally, we observed the spontaneous and heterogeneous segregation of the dopant (Ti) across a wide range of length scales. In addition to the modification of the particle surface and the buried grain boundaries within the particle, we reveal that the Ti doping has induced a significant amount of lattice distortions, which, in turn, promotes the robustness of the LCO lattice at high state of charge (above 4.5V). Our result formulates a multi-scale defect engineering strategy that could be applicable to the synthesis of a broad range of energy materials for applications in batteries and beyond.

36 MATERIALS SCIENCE↗

Origins of Conformational Heterogeneity in Peptoid Helices Formed by Chiral N -1-Phenylethyl Sidechains

N-substituted glycines (polypeptoids) containing chiral hydrophobic sidechains are known to fold into biomimetic alpha helices. These helix formers often produce conformationally heterogeneous structures and are difficult to characterize at a sub-nanometer resolution. Previously, peptoid N-1-phenylethyl (S)-enantiomer sidechains (Nspe) were inferred from various experiments to form right-handed helices and (R)-enantiomers (Nrpe), left-handed helices. Prior computational work for N(s/r)pe oligomers has struggled to reproduce this trend. Herein, quantum mechanics calculations and molecular dynamics simulations are used to understand the origins of this discrepancy. Results from DFT and molecular mechanics calculations on a variety of Nspe and Nrpe oligomers as a function of chain length are in agreement, showing that Nspe and Nrpe prefer left- and right-handed helices, respectively. Additional metadynamics simulations are used to study Nrpe and Nspe oligomers folding in water. These results show that the free-energy driving forces for assembly into a helical backbone configuration are very small (within ~k B T). Lastly, we compare DFT calculations for other experimentally characterized peptoid sidechains, N(r/s)sb, N(r/s)tbe, and N(r/s)npe. In this analysis, we show that peptoid sidechains determined to be more robust experimentally (tbe and npe) have helical preferences opposite the trend seen in less robust assemblies formed by N(r/s)pe and N(r/s)sb chemistries. Furthermore, the more robust tbe and nnpe favor the (S)-enantiomer to right-handed and the (R)-enantiomers to left-handed helices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Human Factors and Technologies Design to Improve User Acceptance of Pooled Rideshare for Increasing Transportation System Energy Efficiency

This multi-year project delivered a comprehensive, human-factors-driven framework to understand, model, and improve pooled rideshare (PR) adoption in the United States. Through three large-scale national survey studies involving more than 16,000 participants across multiple cities and demographic groups, the research established one of the most extensive datasets to date on user perceptions, behavioral barriers, and service expectations related to pooled rideshare. These data revealed key human factors barriers of user acceptance of PR and suggested potential actionable experience optimizations that could lead to increased PR usage. This foundational knowledge guided the development of novel human-factors models and behavioral choice models that quantify how psychological, demographic, and trip-level factors influence willingness to pool. Building on these empirical insights, the project developed advanced behavioral modeling tools, including mixed logit and integrated choice and latent variable models, to capture both observable and latent influences on PR adoption. These models significantly improved the ability to predict riders’ acceptance of pooled trips, explaining choice heterogeneity through latent constructs such as safety, service experience, privacy concerns, time sensitivity, and environmental attitudes. Together, these models provide a robust analytical foundation for designing PR systems that more effectively meet user needs. The project translated human-factors insights and behavioral models into actionable technology innovations by extending POLARIS—an agent-based, activity-based travel simulation platform—into a fully functional pooled rideshare simulation environment. New PR modules, acceptance models, and regional scenarios were implemented for Greenville, SC and Austin, TX, enabling high-fidelity validation of algorithmic strategies under realistic demand and traffic conditions. The simulation platform supported the development and evaluation of adaptive discount-based assignment algorithms, enhanced willingness-to-pay formulations, demographic-aware incentive mechanisms, and a proactive joint assignment and repositioning strategy. Simulation results demonstrated substantial gains in pooling uptake, average vehicle occupancy, energy efficiency, and fleet profitability. In Greenville, pooling adoption more than doubled, while reductions in vehicle-miles traveled and energy consumption were significant. In Austin, pooling improvements were achieved with minimal service-quality trade-offs, and profitability increased across all fleet sizes. Through this research, we developed a comprehensive understanding of the human factors barriers that limit user acceptance of pooled rideshare services. These insights enabled the design of human-factors-aware pooled rideshare technologies that more effectively address user concerns and improve adoption rates. By integrating these models into an advanced agent-based simulation framework, we demonstrated that higher adoption of pooled rideshare can lead to measurable improvements in energy efficiency and system performance. Together, these contributions establish a validated pathway from human-centered analysis to technology development and energy-saving outcomes, supporting national goals for more sustainable and efficient mobility systems.

Jia, Yunyi↗

Orbital Chern Insulator and Quantum Phase Diagram of a Kagome Electron System with Half-Filled Flat Bands

In this work, we study the quantum phase diagram of electrons on kagome lattice with half-filled lowest flat bands by considering the antiferromagnetic Heisenberg interaction J, and short-range Coulomb interaction V. In the weak J regime, we identify a fully spin-polarized phase. The presence of finite V drives a spontaneous chiral current, which makes the system an orbital Chern insulator by contributing an orbital magnetization. Such an out-of-plane orbital magnetization allows the presence of a Chern insulating phase independent of the spin orientation in contrast to the spin-orbit coupling induced Chern insulator that disappears with inplane ferromagnetism constrained by symmetry. Such a symmetry difference provides a criterion to distinguish the physical origin of topological responses in kagome systems. The orbital Chern insulator is robust against small coupling J. By further increasing J, we find that the ferromagnetic topological phase is suppressed, which first becomes partially polarized and then enters a nonmagnetic phase with spin and charge nematicity. The frustrated flat band allows the spin and Coulomb interaction to play an essential role in determining the quantum phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Scalable Interpolation on GPUs for Thermal Fluids Applications

The need for efficient interpolation arises in many use cases for spectral element simulations. Towards this end, we present a fast and robust GPU-enable interpolation utility for NekRS, a spectral-element, Navier-Stokes solver. This utility is based on the findpts and findpts eval routines from the gslib library. To demonstrate the functionality, we added support for particle tracking. The GPU implementation of the interpolation allows for tracking 100 times more particles for the same computational cost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analysis into Asymptotic Convergence to Full Nonlinear Solutions and Exploration of the Implication of Numerical Operator Mutation of Differential Systems

A robust, sufficiently accurate and practical hydrodynamic simulation toolset is required as a key component of the modeling and simulation of air-gap electrostatic discharge events. This work was performed to complement these ongoing efforts. In particular, hydrodynamic simulations must be vetted to ensure they are robust and sufficiently accurate over relevant characteristic scales. Verification models were generated in order to cultivate the technical knowledge and expertise needed to properly create, implement and execute numerical simulations. Furthermore, this effort was utilized extensively to educate students on the mathematical and numerical principles underlying hydrodynamic simulations. This education opportunity, provided in a holistic and rigorous manner, has greatly benefited developing scientists and engineers with the necessary understandings and toolsets required to excel at accomplishing the task at hand, and, more generally, it has enabled them to generate key programmatic deliverables. This report articulates several subtilties; specifically, how perturbations, nonlinear behavior, and dissipative mechanisms influence numerical stability, how to properly structure mathematical and numerical solutions, and how to properly generate error estimation/assignment. A more rigorous discussion of the consequences of such topics can be found in the body of this report in Chapters 2 and 3 with qualitative findings discussed in Chapter 4.

97 MATHEMATICS AND COMPUTING↗