Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Analytical”

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

Analytical closure to the spatially-filtered Euler equations for shock-dominated flows

To ensure numerical stability in the vicinity of shocks, a variety of methods have been used, including shock-capturing schemes such as weighted essentially non-oscillatory schemes, as well as the addition of artificial diffusivities to the governing equations. Centered finite difference schemes are often avoided near discontinuities due to the tendency for significant oscillations. However, such schemes have desirable conservation properties compared to many shock-capturing schemes. The objective of this work is to derive all necessary viscous/diffusion terms from first principles and then demonstrate the performance of these analytical terms within a centered differencing framework. The physical Euler equations are spatially-filtered with a Gaussian-like filter. Sub-filter scale (SFS) terms arise in the momentum and energy equations. Analytical closure is provided for each of them by leveraging the jump conditions for a shock. No SFS terms are present in the continuity or species equations. Here, this approach is tested for several problems involving shocks in one and two dimensions. Implemented within a centered difference code, the SFS terms perform well for a range of flow conditions without introducing excessive diffusion.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Analytical workflow dependence of experimental observables for uranium chemistries

In this work, we investigated how the sequencing of laboratory analytical methods used for chemical and morphological characterization influences analytical findings for particulate materials relevant to the nuclear fuel cycle, including UO2, U3O8, studtite (UO2O2·4H2O), and β-UO3, in the context of nuclear forensic analysis. Particles of each chemistry obtained from consistent production batches were exposed to Raman spectroscopy and scanning electron microscopy in varying orders to elucidate how the order in which the techniques are applied influences morphological and chemical observations as a function of particle size. The results indicate that particles from all four chemistries exposed to high-resolution electron imaging before Raman spectral analysis demonstrate optical vibrational spectral changes that reduce accurate interpretation of the underlying chemistry via Raman spectral analysis. We hypothesize that these changes are due to the thermal load of the electron beam imparted to the sample being unable to be dissipated by materials with poor thermal conduction properties. Results from this study will aid in determining best practices for forensic analysis procedures to reduce uncertainty in chemical determination of unknown particulate samples.

Manns, Rebecca [University of Nevada, Las Vegas]↗
DOE OSTI2020

A data analytical approach for assessing the efficacy of Operational Technology active defenses against insider threats

In recent years, the need for Operational Technology (OT) defenses has been recognized, serving as an additional line of defense when Information Technology (IT) defenses are bypassed. This is no longer considered an uncommon possibility when dealing with advanced persistent threat (APT) actors expected to be state-sponsored and receiving insider assistance. In these extreme adversarial situations, OT defenses aim to provide another layer of defense for the system, introduced directly at the physical process level, as described by the sensors data, the system model, and control actions. Just like IT defenses, two schools of thought, i.e., passive and active defenses, have emerged to address this challenge. In active defenses, representing the focus of this paper, known signatures, synthesized based on the system's unique characteristics, are inserted into the system. In contradistinction, passive methods rely solely on observing system behavior in search of patterns of normal behavior with deviations thereof representing abnormal behavior. In their most sophisticated implementations, both passive and active defenses rely on the use of data analytics to identify the patterns and synthesize the observed and/or inserted signatures. Past research has shown that passive defenses may be bypassed by APT actors relying on data analytics and their intimate knowledge of the system to evade detection by respecting the patterns identified by the defenders. Here this manuscript explores the use of active defenses under the assumption that the attacker has privileged access to the system, including access to the system's model and sensors data. Specifically, this manuscript assesses the ability of active defenses to remain invisible to the attackers, and discusses the associated challenges that must be addressed to ensure their resiliency against APT actors.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2021

Simple, general, realistic, robust, analytic tokamak equilibria. Part 2. Pedestals and flow

Part 1 described a wide range of analytic tokamak equilibria modelling smooth limiter surfaces, double- and single-null divertor surfaces, arbitrary aspect ratio, elongation, triangularity and beta. Part 2 generalizes the analysis to further include edge pedestals and toroidal flow. Specifically, edge pedestals are allowed in the pressure, pressure gradient and toroidal current density. Also, an edge-localized contribution to the bootstrap current is treated. In terms of flow, analytic solutions are obtained for two cases: a $\gamma = 2$ adiabatic and a $\gamma = \infty $ incompressible energy conservation relation.

Physics↗
DOE OSTI2026

Analytic model for neutral penetration and plasma fueling

Neutral atoms recycled from wall interaction interact with confined plasma, thereby refuelling it, most strongly in the region closest to the wall. This occurs near the X-point in diverted configurations, or else near the wall itself in limited configurations. A progression of analytic models is developed for neutral density in the vicinity of a planar or linear source in an ionising domain. First-principles neutral transport simulations with DEGAS2 are used throughout to test the validity and limits of the model when using equivalent sources. The model is further generalised for strong plasma gradients or the inclusion of charge exchange. An important part of the problem of neutral fuelling from recycling is thereby isolated and solved with a closed-form analytic model. A key finding is that charge exchange with the confined plasma can be significantly simplified with a reasonable sacrifice of accuracy by treating it as a loss. The several assumptions inherent to the model (and the simulations with which it is compared) can be adapted according to the particular behaviour of neutrals in the divertor and the manner in which they cross the separatrix.

fusion plasma↗
DOE OSTI2025

Ionic Liquids in Analytical Chemistry: Fundamentals, Technological Advances, and Future Outlook

The development of new analytical methods most often focus on novel materials used to impart selectivity or sensitivity to the protocol. Ionic liquids (ILs) are a class of solvents that have been extensively explored as promising materials for various applications and continue to be explored due to their tunable physicochemical properties. These materials possess melting temperatures below 100 °C and can interact with analytes through a multitude of interactions afforded by their readily tunable chemical structure. These interactions include electrostatic, dispersive, hydrogen bonding, π–π, and dipolar interactions and can be modulated or strengthened based on the functional groups present within the chemical structure. ILs consist predominately of organic cations and either inorganic or organic anions, both of which can be functionalized with desired moieties. Common cation and anions found in IL chemical structures are presented in Figure 1. The unique polarity afforded by the ionic structure has also led to their increasing use in areas including sample preparation, chemical separations, electrochemistry, mass spectrometry, and spectroscopy.

Zeger, Victoria R. [Iowa State Univ., Ames, IA (Un↗
DOE OSTI2026

Analytic Nuclear Gradients Including Oriented External Electric Fields in a Molecule-Fixed Frame

Electric-field-assisted chemistry has attracted much attention in recent years, particularly in the context of oriented external electric fields for controlling molecular structure and reactivity. Such fields have been explored in a wide range of applications, including switching materials, nanoparticles, controllable catalysts, medicines, and clinical therapies. However, the determination of fixed fields in the laboratory frame becomes ineffective for flexible molecules, as conformational changes can significantly alter the relative orientation between the applied field and molecular structure. In this work, we propose two molecular reference frames─the principal axis frame and the local reference frame─to define oriented electric fields within the molecular framework. These coordinate systems powerfully eliminate ambiguities in the relative orientation between the applied field and the molecule. Analytic nuclear gradients in the presence of external electric fields are derived and implemented, with an initial application to field-dependent geometry optimizations of cis - and trans -formanilide. Analysis of the resulting field-induced equilibrium structures reveals distinct structural responses, validating the accuracy and robustness of the proposed formalism. The analytic gradient framework enables systematic investigations of molecular properties and reactivity under arbitrarily oriented electric fields, opening new opportunities for computational modeling and rational design in electric-field-controlled chemistry.

electric fields↗
DOE OSTI2021

Numerical and Analytical Modeling of Flow Partitioning in Partially Saturated Fracture Networks

Infiltration processes in fractured-porous media remain a crucial, yet not very well understood component of recharge and vulnerability assessment. Under partially-saturated conditions flows in fractures, percolating fracture networks, and fault zones contribute to the fastest spectrum of infiltration velocities via preferential pathways. Specifically, the partitioning dynamics at fracture intersections determine the magnitude of flow fragmentation into vertical and horizontal components and hence the bulk flow velocity and dispersion of fracture networks. In this work, we derive an approximate analytical solution for the partitioning process and validate it using smoothed particle hydrodynamics simulations. The developed transfer function is conceptually based on simulation results and laboratory experiments carried out in previous works. It allows to efficiently simulate flow through fracture networks with simple cubic structure and arbitrary number of fractures and aperture sizes via linear response theory and convolution of a given input signal. We derive a non-dimensional bulk flow velocity ($\widetilde{v}$) and dispersion coefficient ($\widetilde{D}$) to characterize fracture networks in terms of dimensionless horizontal and vertical time scales $\tau_m$ and $\tau_0$. The dispersion coefficient is shown to strongly depend on the horizontal time scale and converges towards a constant value of $0.08$ within reasonable fluid and geometrical parameter ranges, while the non-dimensional velocity exhibits a characteristic $\widetilde{v} \sim \tau_m^{-1/2}$ scaling. Given that hydraulic information is often only available at limited places within (fractured-porous) aquifer system, such as boreholes or springs, our study intends to provide a rudimentary analytical concept to potentially reconstruct internal fracture network geometries from external boundary information, e.g., the dispersive properties of discharge (groundwater level fluctuations).

42 ENGINEERING↗
DOE OSTI2022

Leakage From Coexisting Geologic Forcing and Injection‐Induced Pressurization: A Semi‐Analytical Solution for Multilayered Aquifers With Multiple Wells

Abstract Abnormal fluid pressures (above or below hydrostatic pressure) can develop and persist in sedimentary basins. The common occurrence of abnormal pressures may cause challenges for project permitting of geological carbon sequestration (GCS), particularly in reservoirs with pre‐injection overpressure. The leaky wells that may exist in some sedimentary basins can provide flow paths between deep brine aquifers and shallower freshwater aquifers. Pre‐injection relative overpressures can cause brine leakage through leaky wells even before any injection occurs. The tendency for flow through leaky wells is coupled with the process of pressure dissipation that occurs through aquitards. Specifically, with non‐zero permeability, aquitards can dissipate pressure over large areal extents and thereby reduce leakage rates through leaky wells. This study presents development of a semi‐analytical solution for hydraulic head and brine leakage in multilayered aquifer–aquitard systems with geologic pressure forcing. The geologic forcing that causes abnormal pressures in the multilayered system can coexist with any number of injection, extraction, and leaky wells that also affect fluid pressure. The semi‐analytical model is applied to explore how leakage through leaky wells varies as functions of pressurization rate, along with aquitard and leaky well properties in an overpressured multilayered system. The results show that although injection‐induced pressures can dissipate rapidly through suitably permeable aquitards, coexisting geologic forcing may create sustained rates of brine leakage into freshwater aquifers through leaky wells. In GCS, a very low‐permeability aquitard with high capillary entry pressure to free‐phase CO 2 is desired to serve as the caprock to prevent leakage of CO 2 from the storage reservoir. Nevertheless, the results from this study show that the brine leakage impact to shallow freshwater aquifers through leaky wells might substantially decrease with increasing aquitard permeability values, as long as small aquitard permeability and high capillary entry pressure serve to prevent CO 2 leakage.

Cihan, A.↗
DOE OSTI2023

One-Parameter Analytical Derivation in Modified Budyko Framework for Unsteady-State Streamflow Elasticity in Humid Catchments

The changes in climate and catchment properties have altered the hydrological processes significantly at different spatiotemporal scales around the world. In particular for finer time scales, changes in water storage, which has been commonly neglected for long-term temporal scales, may play an important role on hydrological processes. Nevertheless, few studies addressed this question in modifying the Budyko framework, with most of them remaining a steady-state assumption on catchment characteristics. Here, in this work, we derive new analytical formulas of unsteady-state streamflow elasticity in a modified Budyko framework, incorporating both storage change and one specific parameter for catchment properties. We study 78 humid catchments in the USA with simulation data of daily hydrological processes from a probability-distribution-based hydrological model (the modified HyMOD). As indicated by results, the annual storage change ratio is linearly correlated with the annual aridity index, and this relationship can be used to estimate elasticity coefficients with our formulas. The estimated elasticity coefficients perform well in simulating the annual streamflow with the power-law model. For different catchments, variability of the unsteady-state elasticity is higher than that of the steady-state elasticity. Unsteady-state streamflow coefficients show significant linear correlation with catchment properties, such as the average slope, average elevation, and catchment area. This study provides a new analytical approach to investigate the interannual stability of catchments with varying climate and catchment properties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Analytical Model for the Higher Order Moments of Midlatitude Atmospheric Temperature Distributions

Observed distributions of atmospheric temperature are non-Gaussian. Therefore, moments beyond variance are necessary in determining the frequency of extreme temperature events. Here we propose a simple kinematic model for atmospheric mid-latitude temperature variability based on symmetric advection from a non-symmetric background temperature profile. We then use this model to derive analytical expressions for the higher order moments of temperature distributions. Our results show that nonzero skewness and kurtosis arise due to the nonlinearity of the time-mean meridional temperature profile. The analytical model matches an idealized Held-Suarez atmospheric model, indicating nonlinearity of time-mean temperature in latitude is the dominant contribution to nonzero skewness and kurtosis in synoptic temperature variations. Model analysis further shows decrease in higher order moments due to climate change come roughly equally from changes in mixing length and changes in the background temperature profiles.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Analytical and numerical modeling of optical second harmonic generation in anisotropic crystals using #SHAARP package

Electric-dipole optical second harmonic generation (SHG) is a second-order nonlinear process that is widely used as a sensitive probe to detect broken inversion symmetry and local polar order. Analytical modeling of the SHG polarimetry of a nonlinear optical material is essential to extract its point group symmetry and the absolute nonlinear susceptibilities. Current literature on SHG analysis involves numerous approximations and a wide range of (in)accuracies. We have developed an open-source package called the Second Harmonic Analysis of Anisotropic Rotational Polarimetry (#SHAARP.si) which derives analytical and numerical solutions of reflection SHG polarimetry from a single interface (.si) for bulk homogeneous crystals with arbitrary symmetry group, arbitrary crystal orientation, complex and anisotropic linear dielectric tensor with frequency dispersion, a general SHG tensor and arbitrary light polarization. #SHAARP.si enables accurate modeling of polarimetry measurements in reflection geometry from highly absorbing crystals or wedge-shaped transparent crystals. The package is extendable to multiple interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Analytical model for the motion and interaction of two-dimensional active nematic defects

Here, we develop an approximate, analytical model for the velocity of defects in active nematics by combining recent results for the velocity of topological defects in nematic liquid crystals with the flow field generated from individual defects in active nematics. Importantly, our model takes into account the long-range interactions between defects that result from the flows they produce as well as the orientational coupling between defects inherent in nematics. Our work complements previous studies of active nematic defect motion by introducing a linear approximation that allows us to treat defect interactions as two-body interactions and incorporates the hydrodynamic screening length as a tuning parameter. We show that the model can analytically predict bound states between two +1/2 winding number defects, effective attraction between two –1/2 defects, and the scaling of a critical unbinding length between ±1/2 defects with activity. The model also gives predictions for the trajectories of defects, such as the scattering of +1/2 defects by –1/2 defects at a critical impact parameter that depends on activity. In the presence of circular confinement, the model predicts a braiding motion for three +1/2 defects that was recently seen in experiments, as well as stable and ergodic trajectories for four or more defects.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Quantum sensing of paramagnetic analytes by nanodiamonds in levitated microdroplets and aqueous solutions

Nanodiamonds (ND) hosting negatively charged nitrogen-vacancy (NV-) color centers have received attention for applications in magnetic field, electric field, chemical, and bio-sensing. The versatility of these probes is their excellent room-temperature optical and spin properties, along with their small size, functionalized surfaces and resistance to bleaching, making them ideal as nanoscopic sensors in picoliter volumes (e.g. single cells, but also microcompartments and aerosols). For quantitative ND-NV- sensing of paramagnetic analytes in such contexts, however, there remains an incomplete understanding of how factors related to the aqueous phase environment control detection efficiency. To address this, optically detected magnetic resonance (ODMR) is measured in bulk macroscale solutions and single levitated microdroplets as a function of Gd+3 concentration (340 nM to 1.5 mM), nanodiamond size, pH, competitor ions, and ligands. The ODMR response to [Gd+3] is found to be nonlinear, and pH, ND and sample volume dependent; indicating the detection of Gd+3 requires efficient adsorption of the analyte to the diamond surface. Langmuir adsorption isotherms embedded in a quantitative photophysical model links the ODMR response to adsorption thermodynamics of Gd+3. The equilibrium constant for Gd+3 adsorption to a carboxylated ND surface is determined to be (1 ± 0.5) x 105 M-1 corresponding to a free energy of adsorption of (-28 ± 1) kJ mol-1. These results provide general insight into how complex aqueous and microscale environments impact nanodiamond based quantum sensing modalities, and portend their application as quantitative chemical sensors in microenvironments.

Brown, Emily K↗
DOE OSTI2025

Analytical Identification Method of Generalized Short‐Circuit Ratio Using Phasor Measurement Units

This paper introduces a novel analytical approach for the identification of the admittance matrix and the generalized short-circuit ratio (gSCR) in power systems integrated with renewable energy sources. The proposed method leverages voltage and current measurements from phasor measurement units (PMUs) to construct a least squares objective function, which is then solved using matrix calculus and partial derivatives. Unlike conventional optimization algorithms, this approach provides an analytical solution that substantially reduces data requirements, enabling the efficient and accurate identification of the gSCR with smaller datasets. Additionally, its fixed computational complexity allows for real-time updates as new data are collected, ensuring continuous refinement of the system of equations and enabling rapid, precise gSCR calculations. The method also exhibits strong robustness against measurement noise, making it well-suited for practical applications in dynamic power systems. The combination of reduced data requirements, real-time adaptability, noise robustness and fixed computational load establishes this method as a highly efficient and reliable tool for real-time power system stability analysis. Case studies on an EPRI 36-bus system demonstrate the method's effectiveness, highlighting its accuracy in closely matching true gSCR values, even under diverse disturbances and noisy conditions.

Han, Zelei [Hohai University, Nanjing (China)] (OR↗
DOE OSTI2020

Guest Editorial: Advanced Data-Analytics for Power System Operation, Control, and Enhanced Situational Awareness

Along with the smart grid development, modern power systems are entering a ‘data-intensive’ era. A vast volume of data from power grids is being collected through advanced sensing and communication technologies, such as smart metering data, phasor measurement data, as well as meteorological data (e.g., wind speed and solar irradiance) related to renewable power generation. Such data contains comprehensive information about the power system covering equipment's health status, power grid's static and dynamic characteristics, renewable power generation, customers’ electricity usage pattern, etc. Therefore, advanced data-analytics techniques are needed to convert such data to knowledge for practical applications. In line with the trend of widespread data-driven applications in power systems, this Special Issue aims to present state-of-the-art research works on advanced data-analytics for power system's operation, control, and situational awareness. There are in total twenty-six papers accepted for publication in this Special Issue through careful peer reviews and revisions. Under the overarching theme of data-driven applications in power systems, the selected papers are broadly categorised into five topics. The summary of every topic is given below. You are, however, strongly encouraged to read the full paper if interested.

Xu, Yan↗
DOE OSTI2020

Analytic stability boundaries for compressional and global Alfvén eigenmodes driven by fast ions. I. Interaction via ordinary and anomalous cyclotron resonances

Conditions for net fast ion drive are derived for beam-driven, sub-cyclotron compressional (CAE) and global (GAE) Alfvén eigenmodes, such as those routinely observed in spherical tokamaks such as NSTX(-U) and MAST. Both co- and counter-propagating CAEs and GAEs are investigated, driven by the ordinary and anomalous Doppler-shifted cyclotron resonance with fast ions. Whereas prior results were restricted to vanishingly narrow distributions in velocity space, broad parameter regimes are identified in this work which enable an analytic treatment for realistic fast ion distributions generated by neutral beam injection. The simple, approximate conditions derived in these regimes for beam distributions of realistic width compare well to the numerical evaluation of the full analytic expressions for fast ion drive. Additionally, previous results in the very narrow beam case are corrected and generalized to retain all terms in ω/ω ci and |k ∥ /k ⊥ |, which are often assumed to be small parameters but can significantly modify the conditions of drive and damping when they are non-negligible. Lastly, favorable agreement is demonstrated between the approximate stability criterion, simulation results, and a large database of NSTX observations of cntr-GAEs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Comparison of modular analytical wake models to the Lillgrund wind plant

Efficient and accurate wake models are required for wind plant performance modeling and the suite of engineering processes that support wind plant layout, control, and monitoring. Although many analytical and engineering wake models with low computational costs have been proposed, their ability to represent the power production of large wind plants in a wide range of atmospheric conditions is not completely understood. The following validation study reviews the underlying theory for analytical wake models, outlines quality control procedures for observational data, and compares model results with observational data from the Lillgrund Wind Plant. Lillgrund makes a valuable case study for wake modeling because of its regular arrangement and the relatively close spacing of constituent wind turbines, which lead to regular and significant wake interactions within the wind plant and the development of deep array flow conditions. Formulations for the velocity deficit, wake-added turbulence, and wake superposition methods are considered in a modular sense, yielding many possible configurations to represent wind turbine wakes, of which seven are examined in detail. Velocity deficit models that account for flow conditions in the near wake are better able to reproduce power production for wind turbines in the transitional region of the wind plant, where wind turbines experience as many as five wakes from upstream turbines. Finally, in the deep array, where power production reaches asymptotic values and wake statistics become quasi-periodic, wake superposition schemes become the largest driver in error reduction; using the linear or maximum wake superposition methods can reduce the relative root mean square error by as much as 40% in the deep array.

17 WIND ENERGY↗