Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Local time stepping”

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

Disruption event characterization and forecasting in tokamaks

Disruption prediction and avoidance is a critical need for next-step tokamaks, such as ITER. Disruption Event Characterization and Forecasting (DECAF) research fully automates analysis of tokamak data to determine chains of events that lead to disruptions and to forecast their evolution allowing sufficient time for mitigation or complete avoidance of the disruption. Disruption event chains related to local rotating or global magnetohydrodynamic (MHD) modes and vertical instability are examined with warnings issued for many off-normal physics events, including density limits, plasma dynamics, confinement transitions, and profile variations. Along with Greenwald density limit evaluation, a local radiative island power balance theory is evaluated and compared to the observation of island growth. Automated decomposition and analysis of rotating tearing modes produce physical event chains leading to disruptions. A total MHD state warning model comprised of 15 separate criteria produces a disruption forecast about 180 ms before a standard locked mode detector warning. Single DECAF event analyses have begun on KSTAR, MAST, and NSTX/-U databases with thousands of shot seconds of device operation using from 0.5 to 1 × 10 6 tested sample times per device. An initial multi-device database comparison illustrates a highly important result that plasma disruptivity does not need to increase as βN increases. Global MHD instabilities, such as resistive wall modes (RWMs), can give the briefest time period of warning before disruption compared to other physics events. In an NSTX database with unstable RWMs, the mode onset, loss of boundary and current control, and disruption event warnings are found in all cases and vertical displacement events are found in 91% of cases. An initial time-dependent reduced physics model of kinetic RWM stabilization created to forecast the disruption chain predicts instability 84% of the time for experimentally unstable cases with a relatively low false positive rate. Instances of the disruption event chain analysis illustrate dynamics including H–L back transitions for rotating MHD and global RWM triggering events. Furthermore, disruption warnings are issued with sufficient time before the disruption (on transport timescales) to potentially allow active profile control for disruption avoidance, active mode control, or mitigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Harmonic Linearized Navier-Stokes Equations for Transition Prediction in Three-Dimensional Flows

The conventional method to predict the onset of laminar-turbulent transition in convectively unstable boundary-layer flows is based on the logarithmic amplification ratio, the so-called N-factor, of the linear instability waves. To calculate the N-factor, the flow variables are decomposed into a laminar basic state solution and the linear disturbances, which are assumed to be harmonic in time. The most commonly used linear stability analysis approaches include the locally parallel linear stability theory (LST) and the non-local, weakly nonparallel parabolized stability equations (PSE). However, these methods do not account for strong streamwise gradients that are encountered in several configurations of interest, as roughness elements, steps, gaps, or corners. To solve the linear evolution of disturbances along such strongly nonparallel regions, the harmonic linearized Navier-Stokes equations (HLNSE) need to be solved. The discretization of the HLNSE for spanwise/azimuthally inhomogeneous laminar basic states yields a linear system of complex arithmetic with a leading dimension of the order of 107 to 108. A combined multithread and multiprocessor algorithm is implemented for the direct solution of such linear system. Results for a supersonic boundary layer over a three-dimensional roughness patch show good agreement with experimental measurements when the evolution of the instability waves over the roughness patch is included via the HLNSE.

Boundary Layer Stability↗

Interconnection of Distributed Energy Resources in the Indian Context: IEEE 1547-2018 Adaptation for Locally-Appropriate Grid Code Development

The rapid growth of distributed energy resources (DERs) in India and ambitious renewable energy goals set forth by the Indian government have created the need for appropriate distribution interconnection standards to ensure continued system stability, reliability, and resilience. Existing DER interconnection codes in India, at the time of publishing this report, lack many of the critical grid support functionalities for DERs which may protect bulk system stability, such as event ride-through provisions, frequency droop control, or appropriate trip thresholds given the increasing importance of distributed generation. The formal adoption of IEEE 1547-2018 by the Bureau of Indian Standards (BIS) marks a significant step towards achieving appropriate interconnection of DERs, though the implementation of such a standard must be accompanied with careful consideration of the current codes, regulations, and operation of the Indian power system. Studies on the performance of the functionalities within 1547-2018, which properly incorporate local operating conditions, can inform decisions related to the adoption and adaptation of the clauses present in the standard. This report outlines several of these key considerations for stakeholders and provides analyses illustrating the performance of grid support functionalities on typical Indian distribution feeders.

1547↗

Extraordinarily large permittivity modulation in zinc oxide for dynamic nanophotonics

The dielectric permittivity of a material encapsulates the essential physics of light-matter interaction into the material's local response to optical excitation. Photo-induced modulation of the permittivity can enable an unprecedented level of control over the phase, amplitude, and polarization of light. Therefore, the detailed dynamic characterization of technology-relevant materials with substantially tunable optical properties and fast response times is a crucial step to realize tunable optical devices. This work reports on the extraordinarily large permittivity changes in zinc oxide thin films (up to similar to 3.6 relative change in the real part of the dielectric permittivity at 1600 nm wavelength) induced by optically generated free carriers. We demonstrate broadband reflectance modulation up to 70% in metal-backed oxide mirrors at the telecommunication wavelengths, with picosecond-scale relaxation times. The epsilon near zero points of the films can be dynamically shifted from 8.5 mm to 1.6 mm by controlling the pump fluence. Additionally, the modulation can be selectively enhanced at specific wavelengths employing metal-backed zinc oxide disks while maintaining picosecond-scale switching times. This work provides insights into the free-carrier assisted permittivity modulation in zinc oxide and could enable the realization of novel dynamic devices for beam-steering, polarizers, and spatial light modulators.

36 MATERIALS SCIENCE↗

Iterative solution of transonic flows over airfoils and wings, including flows at Mach 1

A new method of calculating transonic flows based on a 'rotated' difference scheme is described. It is suitable for the calculation of both two- and three-dimensional flows without restriction on the speed at infinity and is well adapted to computer use. The Murman procedure is modified to eliminate any assumptions about the direction of flow when constructing the difference scheme. The proper directional property is obtained by rotating the difference scheme to conform with the local stream direction. In the hyperbolic region retarded difference formulas are used for all contributions to the streamwise second derivative, producing a correctly oriented positive artificial viscosity. In the absence of a simple implicit scheme in the hyperbolic and elliptic regions, the concept of iterations as steps in artificial time is introduced. Computer testing of this procedure provides numerical confirmation of the existence and uniqueness of weak solutions of the potential equation when a suitable entropy inequality is enforced.

Jameson, A.↗

An Indicator-based Approach to Sustainable Management of Natural Resources (Chapter 12)

Assessing the sustainability of natural resource management choices for agricultural and forest lands requires quantification of potential changes to a set of environmental and socioeconomic indicators selected to characterize reference scenarios relative to projected future scenarios. Correctly framing the questions with local stakeholders is a critical first step in the sustainability assessment, and the questions that can be addressed are often limited by data availability. Selecting and prioritizing indicators with stakeholders to address their needs and concerns improves the likelihood of investment in monitoring and evaluation of those indicators over time. Computational techniques for analyzing interactions between the selected indicators are inherently affected by the scales and formats of the assembled indicator datasets. Data analytics have the potential to improve understanding of the potential synergies and tradeoffs involved with meeting multiple environmental and socioeconomic goals simultaneously, but timely and appropriate indicator datasets are not always available—even in this new era of “big data.” Continued improvements in data science and data analytics are needed to broaden understanding and acceptance of problems and to provide valuable information for natural resource management. Advances in these areas will enable society to design future landscapes that meet multiple objectives, including the provisioning of agricultural and forest resources along with a variety of ecosystem services (e.g., clean water and healthy soils).

Parish, Esther↗

Martian seismicity

During the Viking mission to Mars, the seismometer on Lander II collected approximately 0.24 earth years of observational data, excluding periods of time dominated by wind-induced Lander vibration. The 'quiet-time' data set contains no confirmed seismic events. A proper assessment of the significance of this fact requires quantitative estimates of the expected detection rate of the Viking seismometer. The first step is to calculate the minimum magnitude event detectable at a given distance, including the effects of geometric spreading, anelastic attenuation, seismic signal duration, seismometer frequency response, and possible poor ground coupling. Assuming various numerical quantities and a Martian seismic activity comparable to that of intraplate earthquakes, the appropriate integral gives an expected annual detection rate of 10 events, nearly all of which are local. Thus only two to three events would be expected in the observational period presently on hand and the lack of observed events is not in gross contradiction to reasonable expectations. Given the same assumptions, a seismometer 20 times more sensitive than the present instrument would be expected to detect about 120 events annually.

Goins, N. R.↗

Modeled sub-hourly nitrate concentrations in subsurface water across a salt marsh system in Elkhorn Slough, California

Salt marshes are dynamic hydrologic systems where terrestrial groundwater, terrestrial surface water, and seawater mix due to bi-directional flows and pressure gradients. Salt marshes are known to remove excess nitrate from water thereby improving water quality in coastal systems. However, our knowledge of the short-term processes (e.g., intra-tidal) in salt marshes is limited. Further, we lack knowledge about how these short-term dynamics may change as a function of seasonality. To address this knowledge gap, we develope a sampling methodology that allowed us to continuously measure nitrate in subsurface water at three different depths across three marsh elevations at sub-hourly time steps. Here we present a data table with the results of nitrate in porewater across the three salt marsh elevations (information in the data set). We analyzed nitrate with a s::can spectrophotometer. We used the UV-Vis absorbance of the samples to calculate nitrate concentrations using a local calibration library. More information about the methodology can be found in Grande et al (2023). Seasonal and Tidal Variations in Hydrologic Inputs Drive Salt Marsh Porewater Nitrate Dynamics, Hydrological Processes.

54 ENVIRONMENTAL SCIENCES↗

Elevated temperature fatigue testing of metals

Material characterization and evaluation conducted for the purpose of calculating fatigue crack initiation lives of components operating at elevated temperatures are discussed. The major technology areas needed to perform a life prediction of an aircraft turbine engine hot section component and the steps required for life prediction are outlined. These include: the determination of the operating environment, the calculation of the thermal and mechanical loading of the component, the cyclic stress strain and creep behavior of the material required for structural analysis, the structural analysis to determine the local stress strain temperature time response of the material at the critical location in the component, and from a knowledge of the fatigue, creep, and failure resistance of the material, a prediction of the life of the component.

Hirschberg, M. H.↗

High-Resolution Mapping of Lunar Crustal Magnetic Fields: Correlations with Albedo Markings of the Reiner Gamma Class

During the last eight months of the Lunar Prospector mission (December 1999-July 1999), the spacecraft was placed in a relatively low-altitude (15-30-km perapsis), near-polar orbit that allowed high-resolution mapping of crustal magnetic fields. We report here initial studies of the correlation of locally strong magnetic anomalies with unusual, swirl-like albedo markings of the Reiner Gamma class. Based on this correlation, which is known from earlier studies of Apollo subsatellite magnetometer data, it has been proposed that the swirls represent regions whose higher albedos have been preserved via deflection of the solar-wind ion bombardment by strong crustal fields. This model in turn depends on the hypothesis that solar-wind implanted H is at least one component of the process that optically matures exposed silicate surfaces in the inner solar system . Specifically, it is hypothesized that implanted H acts as an effective reducing agent to enhance the rate of production of nanophase metallic Fe particles from preexisting silicates during micrometeoroid impacts. According to the model, the curvilinear shapes of these albedo markings are caused, at least in part, by the geometry of ion deflections in a magnetic field. The improved resolution and coverage of the Prospector data allow more detailed mapping of the fields, especially on the lunar farside. This permits a more quantitative test of whether all albedo markings of this class are associated with strong local magnetic fields.Only if the latter condition is met can the solar-wind deflection hypothesis he valid. The basic procedure for mapping crustal magnetic fields using Lunar Prospector magnetometer data follows that developed for analysis of Apollo subsatellite magnetometer data. The specific mapping steps are (1) selection of mission time intervals suitable for mapping crustal fields; these are limited essentially either to times when the Moon is in a lobe of the geomagnetic tail or to times when the Moon is in the solar wind but the spacecraft is in the lunar wake; the data are transformed to a radial, east, and north coordinate system with measurements given as a function of spacecraft latitude, longitude, and altitude; (2) visual editing of individual orbit segments selected for minimal external field disturbances; (3) minimization of remaining low-frequency external fields for individual orbit data segments by quadratic detrending; and (4) two-dimensional filtering of individual orbit segments to produce a vector field map along the slightly curved surface defined by the spacecraft altitude; maps of the three field components (radial, east, and north), the field magnitude, and the spacecraft altitude are constructed. For data obtained at low to middle latitudes, the horizontal resolution of the field maps is limited by the orbit-track separation (about 30 km at the equator). Maps of the field magnitude have been constructed within limited selenographic regions based mainly on data acquired in March, April, and May of 1999. This was a time period when the orbit plane was nearly aligned with the Sun-Moon line so that field mapping was possible at times when the Moon was in the solar wind as well as when the Moon was in the geomagnetic tail. Most of the coverage is across the lunar farside. However, a shows an example of a field map produced from solar-wind wake data for a region including Reiner Gamm on western Oceanus Procellarum (location: 58.5W, 7.5N). The contour interval is 3 nT and the mean spacecraft altitude is 18 km to within the accuracy allowed by the resolution of the map (30 km or about 1 deg.); strong magnetic anomalies correlate closely with swirl locations. Individual orbit profiles (whose resolution along the orbit track is comparable to the spacecraft altitude of 18 km) also demonstrate a good correlation of field magnitude with surface albedo. In order to investigate the correlation of magnetic fields with the location of swirl features, we have reexamined available lunar imagery (Lunar Orbiter, Apollo, and Clementine) to identify and map swirl locations within regions where swirls have previously been mapped. In these images, swirls were distinguished from other high-albedo features such as crater rays by their curvilinear shapes and increased visibility in forward-scattered light. Digital maps of swirls identified by all available imagery were then superposed on maps of the field magnitude at the spacecraft altitude. Based upon analysis of these composite magnetic/geologic maps, we draw the preliminary conclusion that swirl features are associated with magnetic anomalies revealed by Lunar Prospector. Detailed maps of these swirl features are currently being constructed for the magnetically strong regions antipodal to the Imbrium, Serenitatis, and Crisium Basins. Additional information contained in the original,

Hood, L. L.↗

Quantum Sensing of Electron Transfer Pathways in Natural Photosynthesis Using Time-Resolved High-Field Electron Paramagnetic Resonance/Electron–Nuclear Double Resonance Spectroscopy

Photosynthetic integral membrane reaction center (RC) proteins capture and convert sunlight into chemical energy via efficient charge separation achieved through a series of rapid, photoinitiated electron transfer steps. These fast electron transfers create an entangled spin qubit (radical) pair that contains detailed information about the weak magnetic interactions, structure, and dynamics of localized protein environments involved in charge separation events. In this work, we investigate how these entangled electron spin qubits interact with nuclear spins of the protein environment using the high spectral resolution of 130 GHz electron paramagnetic resonance (EPR) and electron-nuclear double resonance (ENDOR). Spectroscopic interrogation enabled the observation and probing of protons located in the electron transfer pathway between the membrane-spanning electron pair P + Q A - (where P + is the primary donor, a special pair of bacteriochlorophylls, and Q A is the primary quinone acceptor) in the bacterial RC. Spectroscopic analysis reveals hydrogen-bonding interactions involved in regulating the route that light-induced electrons travel through the RC protein during charge separation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermal characterization of the build chamber in electron beam melting

Electron beam powder bed fusion, commonly termed electron beam melting (EBM), offers great versatility in multiple-part processes and can produce high quality as-built components due to low residual stresses. The EBM process is complex and requires careful thermal management, including uniform and consistent preheating of the powder bed, in order to ensure quality and consistency of the product. However, most of the simulations in the literature focus on the selective melting stage of the process. As of today, optimal conditions for initial pre-heating temperatures are only available for specific materials. Thus, in order to extend the EBM technology to other desired build materials, a much better understanding of the pre-heating stages is required. In this work, numerical and experimental approaches are combined in order to investigate the effects and sensitivities of heat removal from the build plate during EBM pre-heating stages. For this purpose, a carefully reconstructed numerical model of the build chamber of an ARCAM Q20 + machine is developed. It includes all main parts of the chamber and all relevant heat transfer mechanisms, whereas special attention is paid to radiation heat exchange between various bodies. In order to validate the model, dedicated experiments are performed, in which a system of thermocouples is installed in the build chamber, allowing direct measurement of the local temperatures of the start-plate and heat shields. A good agreement between the simulation and experimental findings is achieved, leading to a better understanding of the thermal processes characteristic to the pre-heating stages. This basic analysis is followed by a representative pre-heating case, where a powder bed is also considered. The energy required to obtain the desired pre-heating temperatures is evaluated, and the role of the powder bed in heat transfer within the chamber is assessed. The pre-heating stage, simulated in the present work, is supposed to create proper conditions for sintering and consequent melting of the powder, which are highly dependent on the local temperatures and heat transfer features. Thus, the findings of the reported work present a step towards a better understanding of the thermal processes that characterize EBM. The reported model can be further used to provide realistic boundary condition inputs for other meso- or macro-scale models as a function of time and geometry. The model can serve also for verification of machine settings, i.e., jump-safe and melt-safe ones which actually provide desired preheat, and for development of settings for new powders.

36 MATERIALS SCIENCE↗

An efficient numerical model for predicting residual stress and strain in parts manufactured by laser powder bed fusion

Abstract Computational modeling of additively manufactured structures plays an increasingly important role in product design and optimization. For laser powder bed fusion processes, the accurate modeling of stress and distortion requires large amount of computational cost due to very localized heat input and evolving complex geometries. The current study takes advantage of a graphics processing unit accelerated explicit finite element analysis code and approximated heat conduction analysis to predict the macroscopic thermo-mechanical behavior in laser selective melting. Adjacent layers and tracks were lumped to reduce the number of time steps and elements in the finite element model. The effects of track and layer grouping on prediction accuracy and solution efficiency are investigated to provide a guidance for a cost-effective simulation. Thin-wall builds from Inconel alloy 625 (IN625) powders were simulated by applying the developed modeling approach to get the detailed residual stress and distortion at a computational speed 50 times higher than conventional approach. Under repeated heating and cooling cycles, a high tensile stress was produced near surfaces of a build due to a larger shrinkage on surface than that in central area. It is also shown that horizontal stresses concentrate near the root and top layers of the IN625 build. The predicted residual elastic strain distribution was validated by the experimental measurement using x-ray synchrotron diffraction.

36 MATERIALS SCIENCE↗

Predicting Accumulations of Ice on Aerodynamic Surfaces

LEWICE is a computer program that predicts the accumulation of ice on two-dimensional aerodynamic surfaces under conditions representative of the flight of an aircraft through an icing cloud. The software first calculates the airflow surrounding the body of interest, then uses the airflow to compute the trajectories of water droplets that impinge on the surface of the body. The droplet trajectories are also used to compute impingement limits and local collection efficiencies, which are used in subsequent ice-growth calculations and are also useful for designing systems to protect against icing. Next, the software predicts the shape of accumulating ice by modeling transfers of mass and energy in small control volumes. The foregoing computations are repeated over several computational time steps until the total icing exposure time is reached. Results of computations by LEWICE have been compared with an extensive database of measured ice shapes obtained from experiments, and have been shown to closely approximate those shapes under most conditions of interest to the aviation community.

Bidwell, Colin↗

Thermal Runaway of Nonflammable Localized High‐Concentration Electrolytes for Practical LiNi 0.8 Mn 0.1 Co 0.1 O 2 |Graphite‐SiO Pouch Cells

Abstract With continuous improvement of batteries in energy density, enhancing their safety is becoming increasingly urgent. Herein, practical high energy density LiNi 0.8 Mn 0.1 Co 0.1 O 2 |graphite‐SiO pouch cell with nonflammable localized high concentration electrolyte (LHCE) is proposed that presents unique self‐discharge characteristic before thermal runaway (TR), thus effectively reducing safety hazards. Compared with the reference electrolyte, pouch cell with nonflammable LHCE can increase self‐generated heat temperature by 4.4 °C, increase TR triggering temperature by 47.3 °C, decrease the TR highest temperature by 71.8 °C, and extend the time from self‐generated heat to triggering TR by ≈8 h. In addition, the cell with nonflammable LHCE presents superior high voltage cycle stability, attributed to the formation of robust inorganic‐rich electrode–electrolyte interphase. The strategy represents a pivotal step forward for practical high energy and high safety batteries.

36 MATERIALS SCIENCE↗

Optimization of 3D controlled ELM-free state with recovered global confinement for KSTAR with n = 1 resonant magnetic field perturbation

Mitigation of deleterious heat flux from edge-localized modes (ELMs) on fusion reactors is often attempted with 3D perturbations of the confining magnetic fields. However, the established technique of resonant magnetic perturbations (RMPs) also degrades plasma performance, complicating implementation on future fusion reactors. Here, we introduce an adaptive real-time control scheme on the KSTAR tokamak as a viable approach to achieve an ELM-free state and simultaneously recover high-confinement (βN ~ 1.91, β p ~ 1.53, and H 98 ~ 0.9), demonstrating successful handling of a volatile complex system through adaptive measures. We show that, by exploiting a salient hysteresis process to adaptively minimize the RMP strength, stable ELM suppression can be achieved while actively encouraging confinement recovery. This is made possible by a self-organized transport response in the plasma edge which reinforces the confinement improvement through a widening of the ion temperature pedestal and promotes control stability, in contrast to the deteriorating effect on performance observed in standard RMP experiments. These results establish the real-time approach as an up-and-coming solution toward an optimized ELM-free state, which is an important step for the operation of ITER and reactor-grade tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Direct Subsurface Measurements through Precise Micro Drilling

Wellbore integrity is a significant problem in the U.S. and worldwide, which has serious adverse environmental and energy security consequences. Wells are constructed with a cement barrier designed to last about 50 years. Indirect measurements and models are commonly used to identify wellbore damage and leakage, often producing subjective and even erroneous results. The research presented herein focuses on new technologies to improve monitoring and detection of wellbore failures (leaks) by developing a multi-step machine learning approach to localize two types of thermal defects within a wellbore model, a prototype mechatronic system for automatically drilling small diameter holes of arbitrary depth to monitor the integrity of oil and gas wells in situ, and benchtop testing and analyses to support the development of an autonomous real-time diagnostic tool to enable sensor emplacement for monitoring wellbore integrity. Each technology was supported by experimental results. This research has provided tools to aid in the detection of wellbore leaks and significantly enhanced our understanding of the interaction between small-hole drilling and wellbore materials.

42 ENGINEERING↗

Elevated temperature fatigue testing of metals

The major technology areas needed to perform a life prediction of an aircraft turbine engine hot section component are discussed and the steps required for life prediction are outlined. These include the determination of the operating environment, the calculation of the thermal and mechanical loading of the component, the cyclic stress-strain and creep behavior of the material required for structural analysis, and the structural analysis to determine the local stress-strain-temperature-time response of the material at the critical location in the components. From a knowledge of the fatigue, creep, and failure resistance of the material, a prediction of the life of the component is made. Material characterization and evaluation conducted for the purpose of calculating fatigue crack initiation lives of components operating at elevated temperatures are emphasized.

Hirschberg, M. H.↗