SEARCH · Engineering Papers

Results for “slicing”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2021

How AI Predicts the Untrained and Unseen

Focus Area: Model predictability improvements (Primary); Data optimization (secondary); Data complexity insights (secondary). The Scientific Challenge: If we believe that a future under extreme conditions will look very differently from today, we can likely agree that ML/AI models trained on past and present datasets will not be adequate to make reliable predictions into the future. This is true for water cycling, as well as biogeochemistry and other Earth system components and behaviors. Additionally, ML/AI models are inherently non-physical. Despite the flourishing success of ML/AI in many applications, such as computer vision, natural language process, and gaming, even the most sophisticated AI models don’t understand the very basic physical laws. Therefore, a natural question is: Can we trust ML/AI based predictions of Earth system behaviors that are fundamentally driven by physical laws? So, are physics models with meticulous process representation a better choice? Not exactly. Physical models, when firmly rooted in first principles, work great at predicting behaviors of systems with a well-defined set of boundary conditions and variables. However, as a complex system, the number of parameters and the degree of complexity and dynamics in processes, coupling, and scale dependent emergent behaviors make the Earth system behaviors very challenging to predict with physical models composed of deterministic laws. In addition, due to the lack of fundamental understandings, physical models often implement empirical correlations derived from observations with biases from locality of data generation. Because correlation is not necessarily causation or comply with first principles, scaling of model predictions beyond locality is often invalid. Beyond the limitation of models, physics or AI, our knowledge of the Earth system is limited by the lack of observational technologies and resources. Insufficient data density, dimensionality and diversity only offer a sliced (or projected) view of the Earth system, e.g. Plato’s Cave analogy, limiting our capability to better understand and represent fundamental processes in models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

High Fidelity Tomographic Density Reconstructions [Slides]

Pressed high explosives have small density variations associated with the pressing process, from a few to sub percent depending on the material. Variations at these levels are known to have an impact on HE performance. Currently, a destructive testing method based on physical sectioning followed by immersion density measurements (slice and dice) is utilized. However, this method has poor spatial resolution limited by the the sectioning size, 1-2 cm. If possible, radiographic methods including Computed Tomography (CT) would be preferred, as they are non-destructive and have higher (less than .1 mm) spatial resolution. However, beam hardening, scattered radiation, and detector effects make density reconstructions at this level extremely challenging.

42 ENGINEERING↗

DOE OSTI2021

Core Characterization of Bakken Shale from the Bedwell 33-52-1-1H Well, Sheridan County, Montana

The computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) at the National Energy Technology Laboratory (NETL) in Morgantown, West Virginia were used to characterize core of the Bedwell 33-52-1-1H well (API 25-091-21920), a wildcat field in Sheridan County, Montana. The primary impetus of this work is to provide a unique dataset to researchers at universities, national laboratories, geological surveys, and other locations for additional analyses. The resultant datasets are presented in this report and can be accessed from NETL's Energy Data eXchange (EDX) online system using the following link: https://edx.netl.doe.gov/dataset/bedwell-33-52-1-1h-well. The equipment and techniques used to characterize the full core were non-destructive, enabling future examinations and analyses to be performed on these cores. However, none of the equipment used was suitable for direct visualization of the pore space in fine-grained structures such as those found in shales; fractures, discontinuities, and millimeter scale features were readily detectable with the methods tested. Imaging with the NETL medical CT scanner was performed on the entire core. Qualitative analysis of the medical CT images, coupled with X-ray fluorescence (XRF), P-wave, and magnetic susceptibility measurements from the MSCL were useful in identifying zones of interest for more detailed analysis. A selection of samples (Table 1) were cored, sliced, and powdered for additional characterization analyses. These analyses did involve destructive subsampling of the full core at discrete sections, and the results are presented here in the context of the larger core description. Selected sections of the core were removed from the full core for detailed analyses using the following destructive techniques: bulk inductively coupled plasma-optical emission spectroscopy (ICP-OES), isolation of carbon and sulfur, Fourier transform infrared (FTIR) spectroscopy, scanning electron microscopy (SEM), and higher resolution CT scan with NETL’s micro-CT scanning systems. The ability to quickly identify key areas for more detailed study with higher resolution will save time and resources in future studies. The combination of methods used provides a multi-scale analysis of the core; the resulting macro and micro descriptions are relevant to many subsurface energy related examinations traditionally performed at NETL.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2022

SCALE Analysis of a Fluoride Salt-Cooled High-Temperature Reactor in Support of Severe Accident Analysis

As part of a US Nuclear Regulatory Commission–sponsored project to assess the modeling and simulation capabilities for accident progression, source term, and consequence analysis for advanced reactor technologies with SCALE and MELCOR, SCALE was used for the modeling and simulation of a fluoride salt-cooled high-temperature reactor (FHR). Based on the preconceptual design for a small modular 236 MWth FHR developed by the University of California, Berkeley (PB-FHR-Mk1), a SCALE model of the PB-FHR-Mk1 reactor core was developed. The reactor was modeled at equilibrium state with different fuel compositions in different regions of the reactor. An iterative approach was used to interpolate and mix the burnup-dependent fuel compositions obtained through the depletion calculation of a core slice model. After demonstrating the applicability of SCALE’s multigroup (MG) approach for the simulation of the PB-FHR-Mk1, the resulting equilibrium core was studied in terms of the power profile, the flux profile, temperature reactivity coefficients, and the xenon reactivity. Furthermore, the tritium production rate in the salt coolant was determined, and the dependence of the one-group cross sections on the burnup and location in the reactor core was studied. The results obtained with SCALE were post-processed to provide the MELCOR team with the core inventory and decay heat of the equilibrium core, a zone-wise power profile, temperature feedback coefficients, the tritium production rate, and the xenon worth.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

MPACT Theory Manual (V.4.3)

This manual presents the theory underlying the three-dimensional (3D) whole-core, pin-resolved neutron transport calculation methodologies employed in the MPACT code. MPACT’s primary goal is to provide accurate sub-pin power distributions in a computationally efficient manner. To accomplish this, MPACT offers several different transport method options. The 2D/1D method, in which 3D problems are decomposed into an axial stack of radial “slices,” is currently the most commonly used. In this option, two-dimensional (2D) planar solutions are provided by the method of characteristics (MOC), and axial solutions are provided via one-dimensional (1D) approximate diffusion, or P3 solutions. The radial and axial solutions are coupled by (i) axial and radial transverse leakages, and (ii) a global 3D coarse mesh finite difference (CMFD) solve, which provides both acceleration and stability to the solution iteration scheme. The subsequent chapters of this manual present a range of topics, including MOC, CMFD, axial nodal transport solvers, 2D/1D, self-shielding, depletion, thermal-hydraulics, and transient methods. Development of the underlying theory of the 2D/1D method is presented, diagrams are included to highlight important algorithmic flow, and important concepts are discussed as appropriate. This manual is intended to be self-sufficient, but references to published articles and other materials are included for further reading. MPACT is a relatively new code, with new capabilities and many computational methods that did not exist until recently. This manual is an even newer document, with chapters written by several different code contributors working under time constraints. For these reasons, the manual is not yet complete and finalized. Future versions of this manual will address current deficiencies.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Search for Rare Processes: Mu2e and KOTO

The Standard Model (SM) of particle physics has been very successful in describing the fundamental forces, the properties of elementary particles, and predicting particles prior to their discovery, for example the top quark and the Higgs boson. Yet the SM does not explain all phenomena, and discoveries of physics processes beyond the SM will allow progress in our understanding of the fundamental nature of matter and forces. The search for rare processes is one of the paths to look for new physics. Two of these searches are for muon to electron conversion and the decay $\text{K}_\text{L}\rightarrow \pi^0 \nu \bar{\nu}$. Our main effort was in the Mu2e experiment under construction at Fermilab, which has the goal of searching for the neutrinoless conversion of a muon into an electron at a single event sensitivity of 2.87 x 10-17 and, in the process, provide evidence of physics processes beyond the Standard Model. We developed the firmware to be used in the Mu2e tracker readout, assisted in the construction and operation of the Mu2e Vertical Slice Test, and provided input to studies of cosmic event reconstructions in Mu2e using simulated data. Under this reward we also concluded our effort in two KOTO experiment in J-PARC, Japan and set a limit on the rate of the rare decay, $\text{K}_\text{L}\rightarrow \pi^0 \nu \bar{\nu}$ using the data collected in the 2016-2018 period. We also maintained and operated the DAQ system used for data collection thru 2021.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Second Target Station Target Segment PIE Dose Estimates

The purpose of the analysis described in this document is to support post irradiation examination (PIE) future decisions as related to the Second Target Station (STS) Project. The STS Project is discussing plans for future PIE studies, and the radionuclide inventories and dose rates included in this study will aid decision makers as they consider the need for various PIE-related machinery. Three configurations of the lasagna-type target segment are evaluated for this analysis: a single bare tungsten plate, a slice of the segment containing a single tungsten plate, and the upstream half of a segment. The lasagna type target segment used in this analysis has a precipitation hardenable nickel alloy (UNS N07718 or inconel 718) shroud, copper (C10100) cladding, and pure tungsten plates. The target segment is evaluated for 10 years of operation followed by two decay periods: 6 months and 2 years. For each decay period, decay gamma dose rates are tallied 30 cm from the upstream, downstream, top, bottom, and sides of the configurations. Radionuclide inventories are also provided in this report separated by material. The highest dose rates calculated in this study are 2872 rem/hr 30 cm away from the top and bottom of the half-target configuration after 6 months of decay.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Literature Review of Recycling Polypropylene and Polyamide 12 Powders for Selective Laser Sintering

Additive manufacturing (AM) is the creation of three-dimensional parts by adding material layer-by-layer based on two dimensional “slices” of a CAD file, without molds or tooling. AM has broken the relationship between part complexity and manufacturing cost. As a general rule, for conventional manufacturing, the more complex the part produced, the more costly its manufacturing. For example, in the production of a bracket for a satellite, moving from conventional manufacturing to AM allowed the part to be consolidated from 4 parts to just one, as well as reducing the weight of the part by 35%. A further benefit, AM can be more cost effective for small lot sizes. To illustrate, consider a case study of white board marker caps presented by Klahn et al. The cost to produce 1,000 units of a new design using SLS was about an order of magnitude lower than conventional manufacturing, as shown in Figure 1.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Photoinitiated Reactions of Molecules and Radicals in Molecular Beams

The UV photochemistry of organic molecules is a fundamental process that governs reactions in the atmosphere, synthetic chemistry, processing of organic aerosols, and biological damage in living tissues. In many organic molecules, the ensuing evolution involves pathways that are in competition, giving rise to different products, isomerization, coupling to other electronic states, and secondary reactive collisions. Because photodissociation is usually fast (picoseconds to microseconds), these processes are far from equilibrium and are controlled by kinetic competition and dynamical forces. The work accomplished was focused primarily on the photochemistry of alpha-keto carboxylic acids, a group of acids produced from natural sources, which are implicated in aerosol formation and biological processes. In the atmosphere, they are destroyed mainly by solar radiation. Studying their photochemistry has been surprisingly difficult because of the complexity of their excited electronic states, and the effect of collisions and secondary reactions. The second project investigated the lifetime of excited electronic states of pyrazine and picoline and had both experimental and theoretical aspects. The work focused on a model molecule, pyruvic acid (PA), because it was hypothesized that it should be a good source of the unstable carbene, methylhydroxycarbene (MHC). PA has an internal hydrogen bond that controls the evolution of its decomposition. The decomposition was studied following excitation to two of its lowest excited states reached by laser irradiation at 351 and 193 nm. The photodissociation dynamics in the absence and presence of collisions was monitored and compared. Understanding the UV photochemistry requires the use of complementary experimental approaches. Two methods were enlisted, which together generated a comprehensive and detailed set of results: (i) The time-sliced velocity map imaging (SVMI) instrument at USC was exploited to determine kinetic energy release of fragments, fast dissociation timescales, and internal state distributions of fragments for which Resonance Enhanced Multiphoton Ionization (REMPI) schemes exist; and (ii) The multiplexed photoionization mass spectrometer (MPIMS) setup developed at the Sandia Combustion Research Facility was used for product discovery, achieved by exploiting tunable narrowband VUV radiation at the Advanced Light Source (ALS), and to follow in real time their subsequent unimolecular and bimolecular reactions. The experimental conditions ranged from collisionless molecular beams to study nascent products to flow reactors of variable pressures to study subsequent bimolecular reactions of the products. The goals stated above were successfully accomplished. By exciting PA to its lowest excited state, MHC was identified as the only primary product and its isomerization to vinylalcohol and acetaldehyde was directly observed in real time. Moreover, we reported the first bimolecular reaction of MHC, which was with acetaldehyde, and identified its reaction product. This paves the way for the study of other reactions of this important carbene intermediate. When excitation was carried out at 193 nm, which imparted much higher energy to PA, many more products were directly observed in real time, some deriving from three-body dissociation. Quantitative branching ratios and secondary reactions of radical products were also determined and analyzed. The studies of pyrazine and picoline focused on the second excited state of these molecules and showed (via ionization studies) that their excited states were very short lived (<100 fs) because of efficient couplings to lower electronic states via vibronic coupling. Theoretical work on modeling the absorption spectrum and decay mechanisms of the excited states are in progress in collaboration with theoreticians.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Reconstruction of the 4D beam matrix

The widely used transverse parameters characterizing particle beams are the Twiss parameters. These parameters can be measured experimentally but they do not fully characterize the beam since they do not account for possible correlations in particle distribution between two transverse coordinates. These correlations may occur due to uncompensated magnetic field at the cathode or misalignment of focusing quadrupoles in the transport beamline. We test a novel diagnostic for diagnosing full 4D beam matrix which may be used to identify such imperfections. The diagnostic is based on transporting the beam through the beamline which includes a quadrupole and a skew quadrupole magnets and measuring the resulting 2D beam distribution at the screen downstream. Such a measurement can be viewed as measuring a 2D projection of the 4D distribution. Different settings of the quads provide measurements of different slices of the phase space. The reconstruction of the original beam matrix from a number of measurements is done using machine learning algorithm, which provides a fast and reliable way of reconstruction for an arbitrary configuration of the scanning beamline. In August 2024, we set up the diagnostic beamline to perform a quadrupole scan of the beam. The setup includes a skew quadrupole, a regular quadrupole, and a screen. The images on the screen were post-processed to remove experimental artifacts and enhance contrast by eliminating background noise outside the core of the distribution=. The rms parameters of the distribution were then calculated and used as inputs for the reconstruction algorithm. This algorithm attempts to determine the initial beam matrix that produces expected images on the screen closely matching the observed images across all quadrupole settings. The algorithm found a solution in which the expected rms parameters closely align with the observations. Validation of the results is planned for FY25.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Resonance Trajectories and Island Tune Spectra in the Tune Domain

Island tune spectra are introduced, beginning with curves of resonance trajectories plotted in the (Q x , Q I ) tune domain plane (where Q x is horizontal tune and Q I is island tune) for many resonance orders N. A slice in the (Q x , Q I ) plane delivers a Q I -spectrum at a fixed horizontal tune, with one spectral line for every appropriate value of N.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Alpha-Quartz Plastic Strength Investigation Via Diffraction Experiments on Novaculite Using A D-Dia and Elastic Plastic Self- Consistent Interpretation [Thesis}

X-Ray Tomography. Porosity is commonly measured using mercury injection (MI) or water immersion porosimeter (WIP). Both MI and WIP utilize pressure to fill open voids with mercury or water respectively. A measure of the volume change of the sample is then used to estimate the percentage of open voids. However, for the purpose of rheological study it is important to get an accurate measure of open and closed voids. X-ray tomography has been selected as it offers a three-dimensional view into the sample that can quantify all pores limited only by the voxel size which is in the micron range for this study. Radiographs for tomography were collected at the Material Science and Technology Division at Los Alamos National Lab using the Carl Zeiss Xradio 520 instrument and Scout-and-Scan version 16.1 operating software. Two samples were imaged, the starting material and the deformed sample SiO2_65. 3001 radiographs were taken of the starting material with a 6 second exposure time using a 4x objective lens. The x-ray beam was set to 60 kilovoltage peak (kVp) and 5 watts. 1901 radiographs were taken of SiO2_65 with a 25 second exposure time using a 10x objective lens. The x-ray beam was set to 80 kVp and 7 watts. Radiograph files were analyzed by Brian Patterson using Avizo. Void and inclusion volumes were output by voxel sized (1.03 μm) slices, used to calculate a total percentage volume for the starting material.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Novel Concepts for High Gradient Acceleration (Final Technical Report)

We have conducted an intensive, pioneering program to demonstrate novel concepts for achieving high gradient acceleration at frequencies from the conventional microwave bands up to the millimeter wave /THz bands. High gradient accelerators hold the promise of smaller and less costly accelerators for applications ranging from the largest scale accelerators used for discovery science down to the smallest accelerators used for industrial, homeland security and medical applications. The research consisted of two major research thrusts: 1.) Structure-based wakefield accelerator (SWFA) research in collaboration with the Argonne Wakefield Accelerator research group and 2.) Millimeter Wave / THz high gradient acceleration in collaboration with SLAC. The specific goals of the research program were: Design novel metallic metamaterial structures that increase the beam-wave coupling for the accelerator mode and reduce the effect of high order modes; test novel metamaterial structures to achieve higher output power, > 1 GW at X-Band (11.7 GHz), in test at the Argonne Wakefield Accelerator (AWA); determine the experimental breakdown threshold for nanosecond-scale pulses at X-Band in testing at the AWA; test a 110 GHz accelerator structure with a field emission electron gun, built at SLAC, using pulses from a 1 MW, 110 GHz gyrotron; design, build and test a 110 GHz quasi-optical, resonant-ring pulse compressor to compress microsecond pulses from the 1 MW gyrotron into > 20 MW, 5 ns output pulses for accelerator structure testing. The proposed research program built on our successes in our research program including: Generation of 510 MW, 2.1 ns (FWHM) pulses at 11.7 GHz from a metallic metamaterial structure in test at the Argonne Wakefield Accelerator using a train of eight 65 MeV electron bunches spaced at 1.3 GHz with a total charge of 280 nC. The metamaterial structure consisted of 100 copper unit cells each consisting of a “wagon-wheel” plate and a spacer plate with a total structure length of 0.2 m. The 510 MW pulse generated an on-axis wakefield of 130 MV/m that could be used to accelerate a trailing witness bunch. Demonstration of coupling of an unprecedented rf power level of 575kW into a 110 GHz accelerator structure using a quasi-optical setup. The standing structure consisted of a central copper cavity located between two matching cavities fed by a TM01 mode. The 6 ns input pulses were sliced from 3 microsecond pulses from the gyrotron using a laser-driven silicon switch. We obtained an unprecedented high gradient up to 230MV/m corresponding to a peak surface electric field of more than 520 MV/m.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

Developing and Managing Data Acquisition Software Using Spack

The Data Acquisition systems of particle physics experiments regularly push the boundaries of high-throughput computing, demanding low-latency collection of data from thousands of devices, collating data into time-sliced events, processing these events and making trigger decisions, and writing the selected data streams to disk. To accomplish these tasks, the DAQ Engineering and Operations department at Fermilab leverages multiple software libraries and builds reusable DAQ frameworks on top. These libraries must be delivered in well-defined bundles and are thoroughly tested for compatibility and functionality before being deployed to live detectors. We have several techniques used to ensure that a consistent set of dependencies can be delivered and re-created at need. We must also support active development of DAQ software components, ideally in an environment as close as possible to that of the detectors. This development often occurs across multiple packages which have to be built in concert and features tested in a consistent and reproducible manner. I will present our scheme for accomplishing these goals using Spack environments, bundle packages, and Github Actions-based CI.

Flumerfelt, Eric [Fermilab]↗

DOE OSTI2025

Design of a Microbunched Electron Cooler Energy Recovery Linac

Microbunched electron Cooling (MBEC), a type of Coherent electron Cooling (CeC), is a possible way to cool high energy protons; such an electron cooler can be driven by an energy recovery linac (ERL). The beam parameters of this design are based on cooling 275 and 100 GeV protons at the Electron-Ion Collider (EIC), requiring 150 and 55\,MeV electrons, respectively. If implemented, a high energy cooler would serve to increase the average luminosity of the collider by mitigating the emittance growth caused by various processes. This ERL is designed to deliver a bunch charge of 1\,nC, an average current of 100\,mA, and strict requirements on the transverse emittance, slice energy spread, and longitudinal distribution profile. This paper covers the current state of the design.

Deitrick, K. [Thomas Jefferson National Accelerato↗

DOE OSTI2026

Tribal Engagement in Transmission Planning

The purpose of this white paper is to provide Tribes with information that (a) helps assess whether transmission access is important for their energy goals, and (b) provides guidance on how to engage in transmission planning, if and when doing so would help the Tribe achieve its energy goals. DOE's Office of Indian Energy notes that "Indian Country contains vast untapped energy resources" (OIE 2023). Individual tribes express a wide array of energy needs and goals, ranging from simple access to basic electricity service and energy efficiency tools, to the development of utility-scale generating plants. Access to transmission is a key part of this multifaceted tribal energy picture. About 2.3 percent of the nation's transmission miles are on tribal lands (OIE 2023). Even within that slice, however, circumstances vary widely. Much of that transmission is for local network delivery (138kV or less). There are some large-capacity 500kV lines, but in many cases the Tribe has limited access to them even when they run across the reservation. Some of the Tribes with the largest land area have little transmission, and some such as Navajo Nation and Hopi Tribe have many homes with no access to electricity. This white paper begins with a high-level summary of transmission issues that can affect a Tribe's energy goals. The discussion includes findings from major transmission-related studies, including the 2023 National Transmission Needs Study (GDO 2023). It then discusses strategies that could help a Tribe safeguard its interests when those interests intersect transmission infrastructure planning.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2025

Surrogate Modelling of 3rd Integer Resonant Extraction at Fermilab Delivery Ring

We present an ongoing work in which a surrogate model is being developed to reproduce the response dynamics of the third-integer resonant extraction process in the Delivery Ring (DR) at Fermilab. This effort is in pursuit of smoothly extracting circulating beam to the Mu2e Experiment s production target, wherein the goal is to extract a uniform slice of the circulating $1e12$ protons in the DR over 25,000 turns (43~ms). The DR contains 3 harmonic sextupoles which excite a third-integer resonance as well as three fast, tune-ramping quadrupole magnets which drive the horizontal tune towards the $29/3$ resonance. In our initial work the surrogate model trains on a semi-analytical simulation provided in the same format as live data. Using Reinforcement Learning (and other potential ML methods), the trained surrogate acts as the environment in which a simple ML control agent could learn to dynamically adjust the quadrupole ramp at 430 break points within the 43 microsecond spill window. The control agent will be hosted on a dedicated Arria 10 FPGA, introducing its own requirements on control agent architecture. In this work we report the accuracy and fidelity of surrogate models in comparison to the response dynamics of the physics simulator.

Narayanan, Aakaash [Fermilab] (ORCID:0000000157944↗

DOE OSTI2024

Beam Dynamics of the Strong Hadron Cooler ERL at the Electron-Ion Collider

The Strong Hadron Cooler Energy Recovery Linac (SHC-ERL) for the Electron-Ion Collider (EIC) is required to produce an electron beam with a bunch charge of 1 nC, an average current of 100 mA, and a beam energy of both 150 and 55 MeV, with strict requirements for the transverse emittance, slice energy spread, and longitudinal distribution. This paper covers the design in detail, including preliminary considerations of higher order effects and diagnostics.

Deitrick, K.↗