Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “monotonicity”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

A Comparison of ARTEMIS Observations and Particle-in-cell Modeling of the Lunar Photoelectron Sheath in the Terrestrial Magnetotail

As an airless body in space with no global magnetic field, the Moon is exposed to both solar ultraviolet radiation and ambient plasmas. Photoemission from solar UV radiation and collection of ambient plasma are typically opposing charging currents and simple charging current balance predicts that the lunar dayside surface should charge positively; however, the two ARTEMIS probes have observed energydependent loss cones and high-energy, surface-originating electron beams above the dayside lunar surface for extended periods in the magnetosphere, which are indicative of negative surface potentials. In this paper, we compare observations by the ARTEMIS P1 spacecraft with a one dimensional particle-in-cell simulation and show that the energy-dependent loss cones and electron beams are due to the presence of stable, non-monotonic, negative potentials above the lunar surface. The simulations also show that while the magnitude of the non-monotonic potential is mainly driven by the incoming electron temperature, the incoming ion temperature can alter this magnitude, especially for periods in the plasma sheet when the ion temperature is more than twenty times the electron temperature. Finally, we note several other plasma phenomena associated with these non-monotonic potentials, such as broadband electrostatic noise and electron cyclotron harmonic emissions, and offer possible generation mechanisms for these phenomena.

Lunar↗
NASA NTRS2014

Influence of Mixed Mode I-Mode II Loading on Fatigue Delamination Growth Characteristics of a Graphite Epoxy Tape Laminate

Mixed mode I-mode II interlaminar tests were conducted on IM7/8552 tape laminates using the mixed-mode bending test. Three mixed mode ratios, G(sub II)/G(sub T) = 0.2, 0.5, and 0.8, were considered. Tests were performed at all three mixed-mode ratios under quasi-static and cyclic loading conditions, where the former static tests were used to determine initial loading levels for the latter fatigue tests. Fatigue tests at each mixed-mode ratio were performed at four loading levels, Gmax, equal to 0.5G(sub c), 0.4G(sub c), 0.3G(sub c), and 0.2G(sub c), where G(sub c) is the interlaminar fracture toughness of the corresponding mixed-mode ratio at which a test was performed. All fatigue tests were performed using constant-amplitude load control and delamination growth was automatically documented using compliance solutions obtained from the corresponding quasi-static tests. Static fracture toughness data yielded a mixed-mode delamination criterion that exhibited monotonic increase in Gc with mixed-mode ratio, G(sub II)/G(sub T). Fatigue delamination onset parameters varied monotonically with G(sub II)/G(sub T), which was expected based on the fracture toughness data. Analysis of non-normalized data yielded a monotonic change in Paris law exponent with mode ratio. This was not the case when normalized data were analyzed. Fatigue data normalized by the static R-curve were most affected in specimens tested at G(sub II)/G(sub T)=0.2 (this process has little influence on the other data). In this case, the normalized data yielded a higher delamination growth rate compared to the raw data for a given loading level. Overall, fiber bridging appeared to be the dominant mechanism, affecting delamination growth rates in specimens tested at different load levels and differing mixed-mode ratios.

Ratcliffe, James G.↗
NASA NTRS2018

Effect of Ex Vivo Ionizing Radiation on Static and Fatigue Properties of Mouse Vertebral Bodies

For a variety of medical and scientific reasons, human bones can be exposed to a wide range of ionizing radiation levels. In vivo radiation therapy (0.05 kGy) is used in cancer treatment, and ex vivo irradiation (25-35 kGy) is used to sterilize bone allografts. Ionizing radiation in these applications has been shown to increase risk of fracture, decrease bone quality and degrade collagen integrity. Past studies have investigated the deleterious effects of radiation on cortical or trabecular bone specimens individually, but to date no studies have examined whole bones containing both cortical and trabecular tissue. Furthermore, a clear relationship between the dose and the mechanical and biochemical response of bone's extracellular matrix has yet to be established for doses ranging from cancer therapy to allograft sterilization (0.05-35 kGy). To gain insight into these issues, we conducted an ex vivo radiation study to investigate non-cellular (i.e. matrix) effects of ionizing radiation dose on vertebral whole bone mechanical properties, over a range of radiation doses (0.05-35 kGy), with a focus on any radiation-induced changes in collagen. With underlying mechanisms of action in mind, we hypothesized that any induced reductions in mechanical properties would be associated with changes in collagen integrity. METHODS: 20-week old female mice were euthanized and the lumbar spine was dissected using IACUC approved protocols. The lumbar vertebrae (L1- S1) were extracted from the spine via cuts through adjacent intervertebral discs, and the endplates, posterior processes, surrounding musculature, and soft tissues were removed (approx. 1.5mm diameter, approx. 2mm height). Specimens were randomly assigned to one of five groups for ex vivo radiation exposure: x-ray irradiation at 0.05, 1, 17, or 35 kGy, or a 0 kGy control. Following irradiation, the vertebrae were imaged using microcomputed tomography (micro-CT) and then subjected to either monotonic compressive loading to failure or uniform cyclic compressive loading. During cyclic testing, samples were loaded in force control to a force level that corresponded to a strain of 0.46%, as determined in advance by a linearly elastic micro-CT-based finite element analysis for each specimen. Tests were stopped at imminent fracture, defined as a rapid increase in strain. The main outcome for the monotonic test was the strength (maximum force); for cyclic testing it was the fatigue life (log of the number of cycles of loading at imminent failure). A fluorometric assay was used on the S1 vertebrae to measure the number of non-enzymatic collagen crosslinks[4]. A one-way ANOVA was performed on mechanical properties and collagen crosslinks; means were compared with controls using Dunnett's method, with a Tukey-Kramer post-hoc analysis when significance was found (p < or = 0.05). RESULTS: Compared to the unirradiated control group, the concentration of non-enzymatic collagen crosslinks was significantly increased for all irradiated groups (p < 0.0001), and being higher by at least 50% (Figure 1a). By contrast, the radiation effects on the collagen were only evident at the higher doses. For irradiation exposures of 17 kGy or more, strength decreased substantially as the radiation level was increased, but no effect was evident below 17 kGy (Figure 1b). There was no significant change in the stiffness or maximum displacement for any radiation dose (p>0.05). The finite element analysis prescribed force level for cyclic loading exceeded the measured (monotonic) strength of the 17 and 35 kGy irradiated groups (mean +/- SD, 20.6 +/- 5.6 N; 13.2 +/- 3.7 N, respectively) and therefore these groups were eliminated from the fatigue study. The fatigue life for the 0.05 and 1 kGy groups were similar to each other and were not statistically significantly different from the control group (Figure 1c).

radiation↗
NASA NTRSDate not supplied

Southern Hemisphere Winter Storms Respond Differently to Low and High CO2 Forcings

Mid-latitude storms have large climate impacts affecting the variability of precipitation, temperature, and winds in the extra-tropics. In the Southern Hemisphere, climate models project an intensification of winter storms and a poleward shift of summer storms by the end of the 21st century. However, previous studies, using idealized models, have shown a non-monotonic behavior of the storms with increasing temperatures. Thus, with the unabated emissions of greenhouse gases into the atmosphere, one might expect a similar non-monotonic behavior beyond the 21st century. To examine this, we here investigate the response of mid-latitude winter storms in the Southern Hemisphere to increasing CO 2 concentrations (up to eight times the preindustrial CO 2 levels) using two state-of-the-art climate models. We find that at high CO 2 levels, winter storms no longer exhibit an intensification across the entire extratropics, as projected by the end of the 21st century. Instead, winter storms shift poleward, similar to the projected response of summer storms by 2100, including a weakening of the storms at low-mid latitudes and an intensification at mid-high latitudes. Out results stresses the importance of investigating the mid latitude flow beyond 2100, to account for any non-monotonic response to increasing GHGs.

CO2 forcings↗
NASA NTRS2024

Decreased Northern Hemisphere Precipitation from Consecutive CO 2 Doublings Is Associated with Significant AMOC Weakening

Previous studies found many climate properties such as northern hemisphere (NH) surface temperature and precipitation respond non-monotonically when CO2 is increased from 1× to 8×CO2 relative to pre-industrial levels. Here, we explore the robustness of the non-monotonicity in the NH precipitation response in 11 coupled climate models. Eight models show a decrease in NH precipitation under repeated CO2 doubling, indicating that the non-monotonic response is a common but not universal result. Although common, the critical CO2 level where the NH precipitation decrease first occurs differs widely across models, ranging from 2×CO2 to 8×CO2. These models also show a prominent weakening in the Atlantic meridional overturning circulation (AMOC) at the same critical CO2 level, with the AMOC weakening leading the precipitation decrease. The sensitivities of NH precipitation and the AMOC to CO2 doublings are positively correlated, especially when the AMOC weakens beyond 10 Sv. This suggests that the differences in models' AMOC response can explain their contrasting NH precipitation responses, where models with a large AMOC weakening have decreased NH precipitation. Regionally, this decrease in NH precipitation is the most prominent over the North Atlantic, Europe and the tropical Pacific. Our results suggest that special care must be taken with the use of pattern scaling to inform regional climate decision-making.

Atlantic meridional overturning circulation↗
DOE OSTI2023

Stochastic projective splitting

Here, we present a new, stochastic variant of the projective splitting (PS) family of algorithms for inclusion problems involving the sum of any finite number of maximal monotone operators. This new variant uses a stochastic oracle to evaluate one of the operators, which is assumed to be Lipschitz continuous, and (deterministic) resolvents to process the remaining operators. Our proposal is the first version of PS with such stochastic capabilities. We envision the primary application being machine learning (ML) problems, with the method’s stochastic features facilitating “mini-batch” sampling of datasets. Since it uses a monotone operator formulation, the method can handle not only Lipschitz-smooth loss minimization, but also min–max and noncooperative game formulations, with better convergence properties than the gradient descent-ascent methods commonly applied in such settings. The proposed method can handle any number of constraints and nonsmooth regularizers via projection and proximal operators. We prove almost-sure convergence of the iterates to a solution and a convergence rate result for the expected residual, and close with numerical experiments on a distributionally robust sparse logistic regression problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Landmark-Warped Emulators for Models with Misaligned Functional Response

Many computer models output functional data, and in some cases, these functional data have similar, but misaligned, shape characteristics. In this paper, we introduce a general approach for building emulators for computer models that output misaligned functional data when key values in the functional response (landmarks) can be easily identified. This approach has two main parts: modeling the aligned (using the landmarks) functional data, and modeling the functions that map the misaligned data to the aligned space (warping functions). As the warping functions are required to be monotonic, we give special attention to modeling monotonic functional response data. We discuss how our approach can be easily applied for a variety of typical emulators, such as Gaussian processes, Bayesian multivariate adaptive regression splines, and Bayesian additive regression trees, and how sensitivity analysis can be performed. We demonstrate our approach by building emulators for two applications: (1) a high-energy-density physics computer model used to simulate inertial confinement fusion ignition experiments, where model outputs are highly misaligned, and (2) a multiphysics continuum hydrocode used to simulate high-velocity impact experiments, where model outputs are only slightly misaligned. In case (1) traditional methods cannot be applied, while in (2) they can be applied, but the proposed method performs significantly better.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Oxidation and pyrolysis of methyl propyl ether

The ignition, oxidation, and pyrolysis chemistry of methyl propyl ether (MPE) was probed experimentally at several different conditions, and a comprehensive chemical kinetic model was constructed to help understand the observations, with many of the key parameters computed using quantum chemistry and transition state theory. Experiments were carried out in a shock tube measuring time variation of CO concentrations, in a flow tube measuring product concentrations, and in a rapid compression machine (RCM) measuring ignition delay times. The detailed reaction mechanism was constructed using the Reaction Mechanism Generator software. Sensitivity and flux analyses were used to identify key rate and thermochemical parameters, which were then computed using quantum chemistry to improve the mechanism. Validation of the final model against the 1-20 bar 600-1500 K experimental data is presented with a discussion of the kinetics. The model is in excellent agreement with most of the shock tube and RCM data. Strong non-monotonic variation in conversion and product distribution is observed in the flow-tube experiments as the temperature is increased, and unusually strong pressure dependence and significant heat release during the compression stroke is observed in the RCM experiments. These observations are largely explained by a close competition between radical decomposition and addition to O2 at different sites in MPE; this causes small shifts in conditions to lead to big shifts in the dominant reaction pathways. The validated mechanism was used to study the chemistry occurring during ignition in a diesel engine, simulated using Ignition Quality Test (IQT) conditions. At the IQT conditions, where the MPE concentration is higher, bimolecular reactions of peroxy radicals are much more important than in the RCM.

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

Modeling the effect of wind speed and direction shear on utility‐scale wind turbine power production

Abstract Wind speed and direction variations across the rotor affect power production. As utility‐scale turbines extend higher into the atmospheric boundary layer (ABL) with larger rotor diameters and hub heights, they increasingly encounter more complex wind speed and direction variations. We assess three models for power production that account for wind speed and direction shear. Two are based on actuator disc representations, and the third is a blade element representation. We also evaluate the predictions from a standard power curve model that has no knowledge of wind shear. The predictions from each model, driven by wind profile measurements from a profiling LiDAR, are compared to concurrent power measurements from an adjacent utility‐scale wind turbine. In the field measurements of the utility‐scale turbine, discrete combinations of speed and direction shear induce changes in power production of −19% to +34% relative to the turbine power curve for a given hub height wind speed. Positive speed shear generally corresponds to over‐performance and increasing magnitudes of direction shear to greater under‐performance, relative to the power curve. Overall, the blade element model produces both higher correlation and lower error relative to the other models, but its quantitative accuracy depends on induction and controller sub‐models. To further assess the influence of complex, non‐monotonic wind profiles, we also drive the models with best‐fit power law wind speed profiles and linear wind direction profiles. These idealized inputs produce qualitative and quantitative differences in power predictions from each model, demonstrating that time‐varying, non‐monotonic wind shear affects wind power production.

Energy & Fuels↗
DOE OSTI2022

$c$-functions in flows across dimensions

We explore the notion of c-functions in renormalization group flows between theories in different spacetime dimensions. We discuss functions connecting central charges of the UV and IR fixed point theories on the one hand, and functions which are monotonic along the flow on the other. First, using the geometric properties of the holographic dual RG flows across dimensions and the constraints from the null energy condition, we construct a monotonic holographic c-function and thereby establish a holographic c-theorem across dimensions. Second, we use entanglement entropies for two different types of entangling regions in a field theory along the RG flow across dimensions to construct candidate c-functions which satisfy one of the two criteria but not both. In due process we also discuss an interesting connection between corner contributions to the entanglement entropy and the topology of the compact internal space. As concrete examples for both approaches, we holographically study twisted compactifications of 4d N =4 SYM and compactifications of 6d N = (2, 0) theories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Entanglement Cost for Infinite-Dimensional Physical Systems

We prove that the entanglement cost equals the regularized entanglement of formation for any infinite-dimensional quantum state ρ ΑΒ with finite quantum entropy on at least one of the subsystems A or B. This generalizes a foundational result in quantum information theory that was previously formulated only for operations and states on finite-dimensional systems. The extension to infinite-dimensional systems is nontrivial because the conventional tools for establishing both the direct and converse bounds, i.e., strong typicality, monotonicity, and asymptotic continuity, are no longer directly applicable. To address this problem, we construct a new entanglement dilution protocol for infinite-dimensional states implementable by local operations and a finite amount of one-way classical communication (one-way LOCC), using weak and strong typicality multiple times. We also prove the optimality of this protocol among all protocols, even under infinite-dimensional separable operations, by developing an argument based on alternative forms of monotonicity and asymptotic continuity of the entanglement of formation for infinite-dimensional states. Along the way, we derive a new integral representation for the quantum entropy of infinite-dimensional states, which we believe to be of independent interest. Our results allow us to fully characterize an important operational entanglement measure—the entanglement cost—for all infinite-dimensional physical systems.

Complexity↗
DOE OSTI2022

Micromechanics of Fracture Propagation During Multistage Stress Relaxation and Creep in Brittle Rocks

Time-dependent rock deformation caused by the initiation and growth of fractures leads to the weakening of the rock mass. Understanding the fracturing mechanisms involved in the time-dependent behavior in brittle rocks is very important and to achieve this goal, a systematic series of three types of experiments was performed on double-flawed prismatic Barre granite specimens under unconfined compression. The first series aimed to identify the failure mechanism in the short-term failure mode under monotonic loading, the, second series involved multistage relaxation (constant strain) experiments to analyze the damage at different strain levels, and the third series explored the fracture propagation under multistage creep (constant load) experiments. The spatial and temporal evolution of cracking mechanisms were evaluated using the acoustic emission (AE) and two-dimensional digital image correlation (2D-DIC) techniques to observe the whole crack growth process as well as the accumulated inelastic strain at the specified region of interest. Results suggest that in the case of multistage creep experiments, the time to failure was less compared to the multistage relaxation, when loaded above the crack damage threshold (CD) estimated from the monotonic testing. Furthermore, the frequency magnitude distribution of the AE events generated in the three loading conditions followed the Gutenberg Richter model. A relatively lower b-value was obtained for the creep experiments, indicative of high energy AE events and faster crack growth. In addition, the AE and DIC results also revealed high evolution of tensile cracks at-different stages of creep and relaxation compared to shear and mixed-mode cracks.

58 GEOSCIENCES↗
DOE OSTI2021

An experimental and kinetic modeling study of NO x sensitization on methane autoignition and oxidation

An experimental and kinetic modeling study of the influence of NO x (i.e. NO 2 , NO and N 2 O) addition on the ignition behavior of methane/‘air’ mixtures is performed. Ignition delay time measurements are taken in a rapid compression machine (RCM) and in a shock tube (ST) at temperatures and pressures ranging from 900–1500 K and 1.5–3.0 MPa, respectively for equivalence ratios of 0.5–2.0 in ‘air’. The conditions chosen are relevant to spark ignition and homogeneous charge compression ignition engine operating conditions where exhaust gas recirculation can potentially add NO x to the premixed charge. The RCM measurements show that the addition of 200 ppm NO 2 to the stoichiometric CH 4 /oxidizer mixture results in a factor of three increase in reactivity compared to the baseline case without NOx for temperatures in the range 600–1000 K. However, adding up to 1000 ppm N 2 O does not show any appreciable effect on the measurements. The promoting effect of NO 2 was found to increase with temperature in the range 950–1150 K, while the sensitization effect decreases at higher pressures. The experimental results measured are simulated using NUIGMech1.2 comprising an updated NO x sub-chemistry in this work. A kinetic analysis indicates that the competition between the reactions ĊH 3 + NO 2 ↔ CH 3 Ȯ + NO and ĊH 3 + NO 2 (+M) ↔ CH 3 NO 2 (+M), the former being a propagation reaction and the latter being a termination reaction governs NO x sensitization on CH 4 ignition. Recent calculations by Matsugi and Shiina (A. Matsugi, H. Shiina, J. Phys. Chem. A. 121 (2017) 4218–4224) for the nitromethane formation reaction CH 3 + NO 2 (+M) ↔ CH 3 NO 2 (+M), together with the recently calculated rate constants for HONO/HNO 2 reactions significantly improve ignition delay time predictions in the temperature range 600–1000 K. Furthermore, the experiments with NO addition reveal a non-monotonous sensitization impact on CH 4 ignition at lower temperatures with NO initially acting as an inhibitor at low NO concentrations and then as a promoter as NO concentrations increase in the mixture. This non-monotonous trend is attributed to the role of the chain-termination reaction ĊH 3 + NO 2 (+M) ↔ CH 3 NO 2 (+M) and the impact of NO on the transition to the chain-branching steps CH 2 O + HȮ 2 ↔ HĊO + H 2 O 2 , H 2 O 2 (+M) ↔ ȮH + ȮH (+M), HĊO ↔ CO + Ḣ followed by CO + O 2 ↔ CO 2 + Ö and Ḣ + O 2 ↔ Ö + ȮH. NUIGMech1.2 is systematically validated against the new ignition delay measurements taken here together with species measurements and high temperature ignition delay time data available in the literature for CH 4 /oxidizer mixtures diluted with NO 2 /N 2 O/NO and is observed to accurately capture the sensitization trends.

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

Woven ceramic matrix composite surrogate model based on physics-informed recurrent neural network

A recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear constitutive behavior of woven ceramic matrix composites (CMCs) driven by matrix damage at multiple length scales. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the multiscale generalized method of cells (MSGMC) approach coupled with a matrix damage model. This coupling permits simulating the nonlinear behavior of woven CMCs based on constituent response at the micro-, meso-, and macroscales. The multiscale repeating unit cell is loaded under non-monotonic conditions including multiple load / unload cycles and tension / compression. The fiber volume fraction as well as the intra- and intertow void volume fractions are also varied in the generation of training data. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input strain sequence and fiber and void volume fractions, the resulting stress versus strain response while satisfying physical constraints such as positive semi-definiteness of the tangent stiffness matrix and linear elastic unloading. Further, the trained surrogate model effectively matches the stress versus strain response and successfully predicts the tangent modulus throughout the loading regime. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex architectures, nonlinear multiaxial material response, and under non-monotonic loading conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Testing and modeling of CLT-to-CLT joints made using hardwood dowels and screws with varying spacing

With the increasing popularity of mass timber, utilizing large-diameter wood screws and hardwood dowels can offer an easy and cost-efficient way to join structural elements. However, the existing research on hardwood dowels is limited to single fastener joints under monotonic loading. No research is available on the group effects of joints with multiple hardwood dowels. Here, this paper presents the results of 96 monotonic and cyclic single-shear plane CLT-to-CLT joint tests with hardwood dowels (Red Oak and Yellow Birch, 25.4 mm diameter) and self-tapping wood screws (8 mm diameter) in Grand-Fir CLT, where screws were inserted both perpendicular and at a 45-degree angle to the grain for comparison to hardwood dowels. The study evaluated the mechanical properties of the joints by changing fastener spacing. Mechanical properties such as yielding and peak strength, displacement, elastic and yielding stiffness, and ductility were evaluated. The results showed equal strength and stiffness properties for hardwood dowels comparable with large-diameter screws. Additionally, the hardwood dowel joints demonstrated moderate ductility properties. The strength and stiffness of the tested joints were compared with analytical equations found in the literature by considering the group action factor of the fasteners. Lastly, a nonlinear force-displacement model was presented and compared to the experimental results and analytical calculations, which gives the possibility to predict and optimize the mechanical behavior of hardwood dowel and screw joints with varying numbers and spacing between fasteners.

36 MATERIALS SCIENCE↗
DOE OSTI2021

On Anomaly Detection for Transactive Energy Systems with Competitive Market

Two anomaly-detection criteria are proposed for transactive energy systems with competitive markets. Participants of transactive energy systems seek an optimal power allocation through hybrid economic-control methods to facilitate the integration of various types of distributed energy resources to power distribution systems. In transactive energy systems, every participant is assumed to be a rational entity, and consumers have diminishing marginal utility and suppliers have increasing marginal cost. With the first proposed anomaly-detection criterion, the monotonicity of marginal cost and marginal utility are examined. The impact of line flow constraints is also taken into consideration. Then, the second anomaly-detection criterion is proposed for TESs with marginal cost and marginal utility which change faster than a certain rate. The second criterion is more accurate than the first one for TES with marginal cost and marginal utility which change faster than a certain rate, but it requires the knowledge of that rate. Neither criteria requires more data than those necessary to find the optimal power allocation and the market-clear price in a transactive energy system. Therefore, the proposed criteria do not disclose any more data than necessary. As the monotonicity of marginal cost and marginal utility in a TES with competitive markets results in convex objective functions in an optimization problem and strongly convex ones when marginal cost and marginal utility changes faster than a certain rate, the two detection criteria are also applicable to anomaly detection of general convex optimization problems. Simulations are carried out to show the efficacy of the proposed criteria to detect anomalies caused by cyberattacks.

Wang, Peng↗
DOE OSTI2024

A block-spectral adaptive H-/$p$-refinement strategy for shock-dominated problems

An adaptive H-/p-refinement strategy using a novel sensor is devised and tested in a block-spectral compressible Euler code equipped with adaptive-mesh refinement (AMR) and high-order flux-reconstruction numerics. At each Gauss quadrature point (or solution point) within each spectral block (or mesh element) the discrete velocity jump ΔU = ∂U/∂y 1 Δy 1 + ∂V/∂y 2 Δy 2 + ∂W/∂y 3 Δy 3 is calculated and normalized by the local speed of sound, a. Here, the grid spacing, Δx i , is calculated in each direction as the distance between auxiliary Gauss-Lobatto points, staggered relative to the solution points. The polynomial order is increased from p = 0 to p = p max in regions of weak compression, (ΔU/a) crit < ΔU/a < 0 and kept at p = p max in regions of flow expansion ΔU/a ≥ 0, while staying at the H = 0 base mesh level. Regions experiencing strong compressions, i.e. ΔU/a < (ΔU/a) crit , are H-refined up to H = H max where H max is applied at the location of maximum compression, ΔU/a = min(ΔU/a) in the domain, while keeping p = 0 to guarantee robustness and monotonicity of the solution in the H refined region. The critical value of (ΔU/a) crit = -0.06 is found to effectively separate smooth and non-smooth solution regions, supported by a 1D detonation initiation test case in ideal gas and a shock-to-detonation transition in high explosives. Using this value, the Sod shock tube, Shu-Osher problem, double Mach reflection and a 2D detonation in a high-explosive are simulated with the proposed adaptive H-/p-refinement. In the Sod shock tube case, p-refinement resolves the (weak) contact discontinuity while H-refinement enhances the grid resolution in the shock exploiting the monotonicity of the p = 0 reconstruction. For the Shu-Osher problem, p-refinement captures the small-scale oscillations trailing the shock that would be otherwise attenuated, while H-refinement triggered by the ΔU-sensor appropriately tracks the shock. In the double Mach reflection problem, H-refinement confines the numerical diffusion around the reflected shock while p-refinement recaptures many physical features trailing the shock. Finally, in the 2D high-explosive detonation case, H-refinement follows the leading shock and resolves the curvature of the detonation wave, while p-refinement adds resolution to the trailing reaction zone. Finally, the proposed methodology is tested in a detonation-wave propagation test case in high-explosives with numerical predictions comparing favorably against experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Spin canting of Ni/CoO/Fe films grown on curved MgO(0 0 1) substrate

Using element-resolved x-ray magnetic circular dichroism (XMCD) and x-ray magnetic linear dichroism (XMLD) measurements, we determined the spin orientations of Ni, CoO and Fe films in Ni/CoO/Fe films grown on a curved MgO(0 0 1) substrate. We find that the vicinal surface of MgO(0 0 1) substrate results in a spin canting towards out-of-plane direction in the Ni and CoO films as a result of the interfacial coupling. The Ni spin canting angle increases monotonically with the vicinal angle in the studied range of 0–17° and the CoO spin canting angle increases more rapidly towards saturation only at a few degrees of the vicinal angle. The uniaxial magnetic anisotropy induced in the Ni layer by the Ni/CoO interfacial coupling is quantitatively determined and is shown to increase monotonically with the vicinal angle. Our result provides a new pathway for tailoring the spin orientation by modifying the substrate surface symmetry in combining with the ferromagnetic/antiferromagnetic interfacial interaction in thin-film based spintronic devices.

36 MATERIALS SCIENCE↗