Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable order”

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 721 records · Page 40

Binary Optimal Control of Single-Flux-Quantum Pulse Sequences

We introduce a binary, relaxed gradient, trust-region method for optimizing pulse sequences for single flux quanta (SFQ) control of a quantum computer. The pulse sequences are optimized with the goal of realizing unitary gate transformations. Each pulse has a fixed amplitude and duration. Here we model this process as an binary optimal control problem, constrained by Schrödinger’s equation, where the binary variables indicate whether each pulse is on or off. We introduce a first-order trust-region method, which takes advantage of a relaxed gradient to determine an optimal pulse sequence that minimizes the gate infidelity, while also suppressing leakage to higher energy levels. The proposed algorithm has a computational complexity of O(p log(p)), where p is the number of pulses in the sequence. We present numerical results for the H and X gates, where the optimized pulse sequences give gate fidelity’s better than 99.9%, in ≈ 25 trust-region iterations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Cost of Fish Exclusion and Passage Technologies for Hydropower

Hydropower represents a reliable source of renewable energy and accounts for approximately 7% of the total electrical generation in the United States. Future expansion of hydropower is likely to be in the form of either smaller new stream development projects or powering existing non-powered dams. For these new projects to be successful, careful analysis of risks, costs, and uncertainty to offset reduced power production as well as ensuring the protection and safe passage of migratory fish to gain public support, will be required. Exclusion and passage are two common approaches to protect fish from entrainment and impingement at hydropower facilities. The thresholds for entrainment risk and requirements for exclusion and passage often differ depending on the species involved, the characteristics of the facility, and the goals of stakeholders. While the costs associated with environmental mitigations represent a large proportion of the total costs required for the licensing of hydropower facilities, little quantitative information is present within the literature regarding the specific costs of fish exclusion and passage. Working with FOA awardee Natel Energy, scientists at Oak Ridge National Laboratory were tasked with assessing the capital construction costs for downstream fish exclusion and passage infrastructure. This report used keyword searches of an existing environmental mitigation cost data set and manual extraction of additional cost data associated with protection, mitigation, and enhancement (PM&E) measures related to positive barrier screening and passage from regulatory licensing documents available in the Federal Energy Regulatory Commission (FERC) eLibrary. This approach yielded a total of 50 PM&E mitigation measures with estimated capital construction costs pertaining to positive barrier screens, 142 pertaining to passage studies, and 26 pertaining to passage-related studies. PM&E measures associated with positive barrier screens represented <10% of the 171 total FERC project dockets available in the data set. These data were highly skewed toward conventional relicensing projects, as <7% were associated with new stream development (NSD) projects. Results from these data indicate highly variable costs associated with fish screening, with flow-normalized costs one to two orders of magnitude higher for screening with the highest exclusion capability (≤0.09 in. spacing) compared with coarser screening (1 to 2 in.). Furthermore, estimated capital costs of passage infrastructure were positively related to the scale of the project based on installed capacity for some, but not all, types of passage. These data provide an initial baseline for estimating exclusion and passage costs for hydropower development and may help developers consider options for more fish-friendly generation technologies, though gaps remain relating to a lack of data, particularly for NSD projects. More data may still be available within the FERC eLibrary, but significant effort will be required to manually identify and extract the data for future analyses.

13 HYDRO ENERGY↗

PSTN-051: Survey Strategy and Cadence Choices for the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST)

A summary of survey strategy and cadence choices, simulated and evaluated by the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) Scheduler Team, prepared for the Survey Cadence and Optimization Committee (SCOC). A large telescope survey, covering the entire visible sky repeatedly every few days in multiple bandpasses over the course of ten years, is the core idea of the LSST. An area of about 20,000 square degrees observed under a wide range of conditions to deep coadded limiting magnitudes in bandpasses ugrizy enables cosmological studies and studies of the Milky Way structure with unprecedented precision; the same survey, when cadenced well, can serve to open new windows into our understanding of transient and variable stars, and extend our knowledge of small bodies throughout the Solar System by orders of magnitude. The outlines of, and some basic necessary requirements for these goals are outlined in the LSST Science Requirements Document (SRD). Finding options for the observing strategy to meet more detailed needs of an even wider range of science goals, as well as building the LSST Scheduler and Metrics Analysis Framework, has been the work of the LSST Scheduler Team with support and input from the astronomical community, including the Community Observing Strategy Evaluation Paper (COSEP), the Call for White Papers, numerous performance metrics, and guidance from the LSST Science Advisory Committee in their Recommendations for Operations Simulator Experiments. This report will help enable the SCOC to make its recommendation for the preferred survey strategy.

79 ASTRONOMY AND ASTROPHYSICS↗

A simple variable coupler for the cryogenic test of SRF cavities

The cryogenic rf tests of SRF cavities in vertical cryostats is typically carried out using fixed-length antennae to couple rf power into the cavity and to probe the energy stored into the cavity. Although variable couplers have been designed, built and used in the past, they are often a complex, costly, not very reliable auxiliary component to the cavity test. In this contribution we present the design and implementation of a simple variable rf antenna which has about 50 mm travel, allowing to obtain about four orders of magnitude variation in Qext -value. The motion of the antenna is driven by a motorized linear feedthrough outside of the cryostat. The antenna can easily be mounted on the most common type of cavity flanges.

Ciovati, G.↗

A Stratigraphic Approach to Characterize the Deposition and Storage of Organic Matter in Reservoir Sediments.

The relationship between carbon burial and sedimentation in reservoirs is unknown, exposing gaps in our fundamental understanding of the transport, processing, and deposition of sediment and organic matter in fluvial and lacustrine systems and contributing to uncertainty in our understanding of the net impact of dams to the global carbon budget. The 2011-2014 removal of two large dams on the Elwha River, Washington State, the largest dam removal yet completed globally, created extensive cutbank exposures of reservoir sediments, allowing the first characterization of the facies architecture of sediments through direct observation in reservoirs worldwide and providing an unparalleled opportunity to 1) assess the relationship between environmental influences, such as and changes in sediment supply, and their expression in the stratigraphic record, 2) assess the relationship between sedimentation processes and detrital organic carbon deposition and storage, and the importance of coarse-grained organic matter and woody debris to the total carbon budget of a reservoir, and 3) apply the insight gained from these reservoirs to evaluate current global estimates of carbon storage in reservoirs and develop a conceptual model of carbon burial in reservoirs to guide further research, as defined by characteristic stratigraphic “types”. Former Lake Mills, the younger, upstream reservoir, was characterized by a tripartite, subaerial Gilbert-style delta which prograded >1 km into the main reservoir from 1927 to 2011. Sediments were composed of coarse-grained topset beds, steeply dipping foreset beds, and a fine-grained, gently dipping prodelta. While individual event horizons were discernible in fine-grained sediments of former Lake Mills, their number and spacing did not correspond to known drawdown or flood events. Former Lake Aldwell, impounded from 1913 to 2011, was initially defined by the rapid progradation of a Gilbert-style, subaerial delta prior to the upstream completion of Glines Canyon Dam. However, the 1927 closure of Glines Canyon Dam upstream caused the delta to evolve to a fine-grained, mouth-bar type delta indicative of low, finer-grained sediment. This evolution, combined with a previously-unrecognized landslide deposit into the upper delta plain, suggests that understanding the exogenic influences on reservoir sedimentation is critical to interpretation and prediction of the sedimentation within individual systems. Former Lake Mills accumulated ~330 Gg of, with depositional-zone average accumulation rates from 229 to 9262 gCm -2 yr -1 , while Former Lake Aldwell accumulated ~ 91 Gg (263 to 2414 gCm -2 yr -1 ). Carbon storage in both reservoirs was dominated by heterogeneous, coarse organic matter and woody debris in the coarse-grained delta slope and relatively coarse-grained prodelta regions of the reservoirs, with little storage in the gravel-dominated, subaerial delta plains. Carbon accumulation in fine-grained lacustrine and prodelta sediments was relatively homogeneous, but turbidity flows from the Gilbert-style delta slope in former Lake Mills delivered significantly more carbon to the prodelta than the mouth-bar style delta of former Lake Aldwell. C:N ratios support interpretation of most organic matter in both reservoirs as allochthonous. Sampling schemes based only on lacustrine and/or prodelta would underestimate of total carbon accumulation by up to 30% in former Lake Aldwell, but the overestimate by up to 47% in former Lake Mills. Global estimates of carbon sequestration rates in reservoir sediments vary by three orders of magnitude, while individual-reservoir estimates vary by four orders of magnitude and over only 37 reservoirs and a literature review of predictive variables suggests weak or contradictory relationships. A conceptual stratigraphic framework of four unique reservoir types suggests that organic matter deposition is intrinsically tied to sedimentation processes and that patterns of carbon storage vary systematically with the stratigraphy of reservoir sediments. Deltaically-dominated reservoirs (whether Gilbert style or shoalwater) appear to store most carbon in their deltaic and prodelta regions, while thalweg-style reservoirs exhibit a bimodal distribution, with allochthonous carbon preferentially routed along the former river thalweg and autochthonous deposited on the former floodplain. Lacustrine-style reservoirs are dominated by suspended sediment deposition and thus relatively homogeneous, but literature suggests these reservoirs are more variable than typically measured. Current methods of reservoir sampling fail to account for this systematic variation and tend to be biased toward fine sediment, suggesting that global reservoir carbon storage is underestimated.

Stratton, Laurel↗

A Methodology for Loading the Advanced Test Reactor Driver Core for Experiment Analysis

The Advanced Test Reactor (ATR) at the Idaho National Laboratory (INL) is one of the premier irradiation facilities in the world. The ATR hosts experiments from a variety of laboratories, companies, and disciplines. Because of its flexibility in providing customized irradiation conditions to experiments, such as coolant chemistry, temperature, and reactor power, analysis for experiments, both from a programmatic and safety perspective, can be a challenging endeavor. Because of these experiment variables, there is always an effort to improve upon the existing analysis methodology in order to provide the sponsor of any given experiment with the most accurate predictions possible. The paper discusses the effort to improve upon the methodology for loading the ATR driver core for experiment analysis. Experiment analysis at INL is generally broken up into three separate disciplines: neutronics, thermal-hydraulics, and structural. The focus of this paper is neutronics analysis. After the experiment has a preliminary design, the neutronics analyst will create a model in the Monte Carlo N-Particle (MCNP) transport code to calculate fuel and material heating rates from MCNP tallies and provide those to the thermal-hydraulics analyst, and calculate reactivity worths, perturbation on the axial fission profile of the ATR driver core by the experiment, void worths, and experiment backup worths. Also, there are usually several programmatic parameters that the experiment sponsor requests that the neutronics analyst calculates. These can include heat generation rates, burnup and displacements per atom (DPA). In order to accurately analyze these parameters, an important input into the physics analyst's MCNP model is the loading of the ATR driver core. The ATR driver core consists of 40, 19-plate highly-enriched uranium (HEU) fuel elements, arranged in a serpentine shape, as pictured in Figure 1. The ATR driver core is responsible for the vast majority of the thermal power of the reactor and the nuclear interactions with the experiments being irradiated in ATR. A unique feature of ATR is the 'power tilt' across the reactor. A power tilt means that the Northwest, Northeast, Southwest, and Southeast Lobes of ATR frequently operate at different powers, providing different irradiation conditions for the experiments in and around the lobes of the reactor. In order to achieve lobe power splits, different fuel elements for each of the 40 positions are chosen by ATR Reactor Engineering so that lobe power splits can be maintained during an operating cycle. A typical fuel loading for an ATR cycle consists of 14 to 18 new fuel elements and 2 to 26 depleted fuel elements. These loadings are typically chosen from three different types of fuel elements. All three ATR fuel element types share a common design, spelled out in Chapter 4 of the ATR Safety Analysis Report (SAR). The fuel element meat is HEU UAlx dispersion fuel sandwiched between aluminum cladding with a length of 121.92 cm (4 ft), with 19 plates per fuel element. The three variations of this common design are designated by the first two letters in their serial numbers: XA, YA, and NB. XA elements are most commonly used. These elements have nominally 1075 g of uranium-235 and 0.66 g of boron-10. The NB elements have the same uranium-235 loading as the XA elements, but no boron in any plates. These elements are used when excess reactivity is needed. The YA elements have the same loading as the XA elements, except that the outermost plate (Plate 19) has no boron or uranium. These elements are used to mitigate aged reflector damage. In order to accurately model the ATR driver core, all three types of fuel elements must be accounted for in the MCNP model. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling the Spectral Energy Distribution of the Active Galactic Nucleus inside NGC 4395

We study the broadband spectral energy distribution (SED) of the prototypical low-mass active galactic nucleus (AGN) in NGC 4395. We jointly model the optical through mid-IR SED with a combination of galaxy and AGN light, and find that on arcsecond scales the AGN dominates at most wavelengths. However, there is still some ambiguity about emission from the galaxy, owing partially to the strong short-term variability of the black hole. We investigate the use of smooth- and clumpy-torus models in order to disentangle the nuclear IR emission, as well as exploring the use of poloidal wind emission to account for the blue spectral slope observed in the near-IR. Even when simultaneously fitting the full optical–IR spectral range, we find that degeneracies still remain in the best-fit models. We conclude that high-spatial-resolution and wider-wavelength coverage with the JWST is needed to understand the mid-IR emission in this complex, highly variable object, which is the best nearby example to provide a blueprint to finding other low-mass AGNs via their mid-IR emission in the future.

79 ASTRONOMY AND ASTROPHYSICS↗

Predicting laser weld reliability with stochastic reduced-order models. Predicting laser weld reliability

Laser welds are prevalent in complex engineering systems and they frequently govern failure. The weld process often results in partial penetration of the base metals, leaving sharp crack-like features with a high degree of variability in the geometry and material properties of the welded structure. Furthermore, accurate finite element predictions of the structural reliability of components containing laser welds requires the analysis of a large number of finite element meshes with very fine spatial resolution, where each mesh has different geometry and/or material properties in the welded region to address variability. We found that traditional modeling approaches could not be efficiently employed. Consequently, a method is presented for constructing a surrogate model, based on stochastic reduced-order models, and is proposed to represent the laser welds within the component. Here, the uncertainty in weld microstructure and geometry is captured by calibrating plasticity parameters to experimental observations of necking as, because of the ductility of the welds, necking – and thus peak load – plays the pivotal role in structural failure. The proposed method is exercised for a simplified verification problem and compared with the traditional Monte Carlo simulation with rather remarkable results.

36 MATERIALS SCIENCE↗

An arbitrary-order Runge–Kutta discontinuous Galerkin approach to reinitialization for banded conservative level sets

Here, we present a Runge-Kutta discontinuous Galerkin method for solving conservative reinitialization in the context of the conservative level set method. This represents an extension of the method recently proposed by Owkes and Desjardins [21], by solving the level set equations on the refined level set grid and projecting all spatially-dependent variables into the full basis used by the discontinuous Galerkin discretization. By doing so, we achieve the full k+1 order convergence rate in the L1 norm of the level set field predicted for RKDG methods given kth degree basis functions when the level set profile thickness is held constant with grid refinement. Shape and volume errors for the 0.5-contour of the level set, on the other hand, are found to converge between first and second order. We show a variety of test results, including the method of manufactured solutions, reinitialization of a circle and sphere, Zalesak's disk, and deforming columns and spheres, all showing substantial improvements over the high-order finite difference traditional level set method studied for example by Herrmann. We also demonstrate the need for kth order accurate normal vectors, as lower order normals are found to degrade the convergence rate of the method.

97 MATHEMATICS AND COMPUTING↗

Modeling the SHG activities of diverse protein crystals

The origins of the diversity in the SHG signal from protein crystals are investigated and potential protein-crystal coverage by SHG microscopy is assessed. A symmetry-additive ab initio model for second-harmonic generation (SHG) activity of protein crystals was applied to assess the likely protein-crystal coverage of SHG microscopy. Calculations were performed for 250 proteins in nine point-group symmetries: a total of 2250 crystals. The model suggests that the crystal symmetry and the limit of detection of the instrument are expected to be the strongest predictors of coverage of the factors considered, which also included secondary-structural content and protein size. Much of the diversity in SHG activity is expected to arise primarily from the variability in the intrinsic protein response as well as the orientation within the crystal lattice. Two or more orders-of-magnitude variation in intensity are expected even within protein crystals of the same symmetry. SHG measurements of tetragonal lysozyme crystals confirmed detection, from which a protein coverage of ∼84% was estimated based on the proportion of proteins calculated to produce SHG responses greater than that of tetragonal lysozyme. Good agreement was observed between the measured and calculated ratios of the SHG intensity from lysozyme in tetragonal and monoclinic lattices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Measurements and interpretations of W ± Z production cross-sections in pp collisions at $\sqrt{s}=13$ TeV with the ATLAS detector

Measurements of integrated and differential cross-sections for W ± Z production in proton-proton collisions are presented. The data collected by the ATLAS detector at the Large Hadron Collider from 2015 to 2018 at a centre-of-mass energy of $\sqrt{s}=13$ TeV are used, corresponding to an integrated luminosity of 140 fb −1 . The W ± Z candidate events are reconstructed using leptonic decay modes of the gauge bosons into electrons or muons. The integrated cross-section per lepton flavour for the production of W ± Z is measured in the detector fiducial region with a relative precision of 4%. The measured value is compared with the Standard Model prediction at a precision of up to next-to-next-to-leading-order in QCD and next-to-leading-order in electroweak. Cross-sections for W + Z and W − Z production and their ratio are presented. The W ± Z production is also measured differentially as functions of various kinematic variables, including new observables sensitive to CP-violation effects. All measurements are compared with state-of-the-art Standard Model predictions from fixed-order calculations or Monte Carlo generators based on next-to-leading-order matrix elements interfaced with parton showers. An effective field theory interpretation of the measurements is performed, considering both CP-conserving and CP-violating dimension-6 operators modifying the W ± Z production. In the absence of observed deviations from the Standard Model, limits on CP-conserving Wilson coefficients are extracted using the transverse mass of the W ± Z system. For CP-violating coefficients a machine learning approach is used to construct an observable with enhanced sensitivity to CP-violation effects.

hadron-hadron scattering↗

A deep learning enabler for nonintrusive reduced order modeling of fluid flows

Here in this paper, we introduce a modular deep neural network (DNN) framework for data-driven reduced order modeling of dynamical systems relevant to fluid flows. We propose various DNN architectures which numerically predict evolution of dynamical systems by learning from either using discrete state or slope information of the system. Our approach has been demonstrated using both residual formula and backward difference scheme formulas. However, it can be easily generalized into many different numerical schemes as well. We give a demonstration of our framework for three examples: (i) Kraichnan-Orszag system, an illustrative coupled nonlinear ordinary differential equation, (ii) Lorenz system exhibiting chaotic behavior, and (iii) a nonintrusive model order reduction framework for the two-dimensional Boussinesq equations with a differentially heated cavity flow setup at various Rayleigh numbers. Using only snapshots of state variables at discrete time instances, our data-driven approach can be considered truly nonintrusive since any prior information about the underlying governing equations is not required for generating the reduced order model. Our a posteriori analysis shows that the proposed data-driven approach is remarkably accurate and can be used as a robust predictive tool for nonintrusive model order reduction of complex fluid flows.

42 ENGINEERING↗

Prognostication of LED Remaining Useful Life and Color Stability in the Presence of Contamination

The reliability of LED products may be affected by both luminous flux drop and color shift. Previous research on the topic focuses on either luminous maintenance or color shift. However, luminous flux degradation usually takes very long time to observe in LEDs under normal operating conditions. In this paper, the impact of a VOC (volatile organic compound) contaminated luminous flux and color stability are examined. As a result, both luminous degradation and color shift had been recorded in a short time. Test samples are white, phosphorconverted, high-power LED packages. Absolute radiant flux is measured with integrating sphere system to calculate the luminous flux. Luminous flux degradation and color shift distance were plotted versus aging time to show the degradation pattern. A prognostic health management (PHM) method based on the state variables and state estimator have been proposed in this paper. In this PHM framework, unscented kalman filter (UKF) was deployed as the carrier of all states. During the estimation process, third order dynamic transfer function was used to implement the PHM framework. Both of the luminous flux and color shift distance have been used as the state variable with the same PHM framework to exam the robustness of the method. Predicted remaining useful life is calculated at every measurement point to compare with the tested remaining useful life. The result shows that state estimator can be used as the method for the PHM of LED degradation with respect to both luminous flux and color shift distance. The prediction of remaining useful life of LED package, made by the states estimator and data driven approach, falls in the acceptable errorbounds (20%) after a short training of the estimator.

Lall, Pradeep↗

Prognostics Health Management Model for LED Package Failure Under Contaminated Environment

The reliability consideration of LED products includes both luminous flux drop and color shift. Previous research either talks about luminous maintenance or color shift, because luminous flux degradation usually takes very long time to observe. In this paper, the impact of a VOC (volatile organic compound) contaminated luminous flux and color stability are examined. As a result, both luminous degradation and color shift had been recorded in a short time. Test samples are white, phosphor-converted, high-power LED packages. Absolute radiant flux is measured with integrating sphere system to calculate the luminous flux. Luminous flux degradation and color shift distance were plotted versus aging time to show the degradation pattern. A prognostic health management (PHM) method based on the state variables and state estimator have been proposed in this paper. In this PHM framework, unscented kalman filter (UKF) was deployed as the carrier of all states. During the estimation process, third order dynamic transfer function was used to implement the PHM framework. Both of the luminous flux and color shift distance have been used as the state variable with the same PHM framework to exam the robustness of the method. Predicted remaining useful life is calculated at every measurement point to compare with the tested remaining useful life. The result shows that state estimator can be used as the method for the PHM of LED degradation with respect to both luminous flux and color shift distance. The prediction of remaining useful life of LED package, made by the states estimator and data driven approach, falls in the acceptable error-bounds (20%) after a short training of the estimator.

Lall, Pradeep↗

Synchrotron X-Ray Microdiffraction Investigation of Scaling Effects on Reliability for Through-Silicon Vias for 3-D Integration

Synchrotron x-ray microdiffraction has been applied to TSV characterization in various studies for nondestructive inspection with submicron resolution due to its high beam intensity and penetration depth. Here, the application of this technique to TSV investigations is examined and the correlation of the plastic deformation to the microstructure and extrusion behavior along with the effect of TSV dimensional scaling is examined. It is shown that the variability of the copper microstructure and resulting TSV behavior requires a larger number of samples in order to report statistically significant observations. The role of the microstructure in creating statistical scatter is demonstrated through microdiffraction measurements of grain orientation correlated with the observed peak widening, which shows that degraded TSV reliability is largely due to the high elastic anisotropy of copper. After taking the statistical variations into account, the scaling effect was clearly observed, with larger plastic deformation in 2μm diameter TSVs than in 5μm diameter TSVs consistent with microstructure variations. This is confirmed by TSV extrusion measurements, which show that the magnitude and statistical spread of the via extrusion for the 2μm diameter TSVs is higher than that of the 5μm diameter TSVs. These findings, confirmed by thermomechanical simulation, demonstrate first that large sample sizes are required in copper TSV investigations due to high variability, which is not improved with scaling.

36 MATERIALS SCIENCE↗

Variability, timescales, and nonlinearity in climate responses to black carbon emissions

Black carbon (BC) particles exert a potentially large warming influence onthe Earth system. Reductions in BC emissions have attracted attention as apossible means to moderate near-term temperature changes. For the first time,we evaluate regional climate responses, nonlinearity, and short-termtransient responses to BC emission perturbations in the Arctic, midlatitudes,and globally based on a comprehensive set of emission-driven experimentsusing the Community Earth System Model (CESM). Surface temperature responsesto BC emissions are complex, with surface warming over land from midlatitudeBC perturbations partially offset by ocean cooling. Climate responses do notscale linearly with emissions. While stronger BC emission perturbations havea higher burden efficiency, their temperature sensitivity is lower. BCimpacts temperature much faster than greenhouse gas forcing, with transienttemperature responses in the Arctic and midlatitudes approaching aquasi-equilibrium state with a timescale of 2–3 years. We find largevariability in BC-induced climate changes due to background model noise. As aresult, removing present-day BC emissions results in discernible surfacetemperature changes for only limited regions of the globe. In order to betterunderstand the climatic impacts of BC emissions, both the drivers ofnonlinear responses and response variability need to be assessed acrossclimate models.

54 ENVIRONMENTAL SCIENCES↗

Equivalence of ADM Hamiltonian and Effective Field Theory approaches at next-to-next-to-leading order spin1-spin2 coupling of binary inspirals

The next-to-next-to-leading order spin1-spin2 potential for an inspiralling binary, that is essential for accuracy to fourth post-Newtonian order, if both components in the binary are spinning rapidly, has been recently derived independently via the ADM Hamiltonian and the Effective Field Theory approaches, using different gauges and variables. Here we show the complete physical equivalence of the two results, thereby we first prove the equivalence of the ADM Hamiltonian and the Effective Field Theory approaches at next-to-next-to-leading order with the inclusion of spins. The main difficulty in the spinning sectors, which also prescribes the manner in which the comparison of the two results is tackled here, is the existence of redundant unphysical spin degrees of freedom, associated with the spin gauge choice of a point within the extended spinning object for its representative worldline. After gauge fixing and eliminating the unphysical degrees of freedom of the spin and its conjugate at the level of the action, we arrive at curved spacetime generalizations of the Newton-Wigner variables in closed form, which can also be used to obtain further Hamiltonians, based on an Effective Field Theory formulation and computation. Finally, we make use of our validated result to provide gauge invariant relations among the binding energy, angular momentum, and orbital frequency of an inspiralling binary with generic compact spinning components to fourth post-Newtonian order, including all known sectors up to date.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On the Preconditioning of a Newton-Krylov Solver for a High-Order reconstructed Discontinuous Galerkin Discretization of All-Speed Compressible Flow with Phase Change for Application in Laser-Based Additive Manufacturing

This dissertation focuses on the development of a fully-implicit, high-order compressible ow solver with phase change. The work is motivated by laser-induced phase change applications, particularly by the need to develop large-scale multi-physics simulations of the selective laser melting (SLM) process in metal additive manufacturing (3D printing). Simulations of the SLM process require precise tracking of multi-material solid-liquid-gas interfaces, due to laser-induced melting/ solidi cation and evaporation/condensation of metal powder in an ambient gas. These rapid density variations and phase change processes tightly couple the governing equations, requiring a fully compressible framework to robustly capture the rapid density variations of the ambient gas and the melting/evaporation of the metal powder. For non-isothermal phase change, the velocity is gradually suppressed through the mushy region by a variable viscosity and Darcy source term model. The governing equations are discretized up to 4th-order accuracy with our reconstructed Discontinuous Galerkin spatial discretization scheme and up to 5th-order accuracy with L-stable fully implicit time discretization schemes (BDF2 and ESDIRK3-5). The resulting set of non-linear equations is solved using a robust Newton-Krylov method, with the Jacobian-free version of the GMRES solver for linear iterations. Due to the sti nes associated with the acoustic waves and thermal and viscous/material strength e ects, preconditioning the GMRES solver is essential. A robust and scalable approximate block factorization preconditioner was developed, which utilizes the velocity-pressure (vP) and velocity-temperature (vT) Schur complement systems. This multigrid block reduction preconditioning technique converges for high CFL/Fourier numbers and exhibits excellent parallel and algorithmic scalability on classic benchmark problems in uid dynamics (lid-driven cavity ow and natural convection heat transfer) as well as for laser-induced phase change problems in 2D and 3D.

42 ENGINEERING↗