SEARCH · Engineering Papers
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.
The project supported the development of new methodologies for performing matrix computations that form key building blocks in modern scientific computing, such as low rank approximation of matrices, and efficient representations of global operators that arise in simulations of physical phenomena.
This project, Activity Planning with Resources for the Exploration of Space (APRES), uses a mixed-integer linear program (MILP) to solve planning problems. This work enables APRES to interpret a model file and output a solution with improved human readability. A plan model is optimized using a MILP solver and the best solution is taken. Once a plan is generated, it is parsed allowing it to retain only desired information and modified for swift human readability.
LCLS-II is an upgrade project for the linear coherent light source (LCLS) at SLAC. The LCLS-II linac consists of thirty-five 1.3 GHz and two 3.9 GHz superconducting RF (SRF) continuous wave (CW) cryomodules with high quality factor cavities. Cryomodules were produced at Fermilab and at Jefferson Lab in collaboration with SLAC. Fermilab successfully completed the assembly, testing and delivery of seventeen 1.3 GHz and three 3.9 GHz cryomodules. LCLS-II-HE is a planned upgrade project to LCLS-II. The LCLS-II-HE linac will consist of twenty-three 1.3 GHz cryomodules with high gradient and high quality factor cavities. This paper presents LCLS-II-HE cryomodule production plans, emphasizing the improvements done based on the challenges, mitigations, and lessons learned from LCLS-II.
A brute-force, long-time gyrokinetic simulation of plasma profile evolution in magnetic fusion devices is not desirable due to large computational resource requirements and a possible accumulation of numerical error. The equation-free projective integration method of Keverekidis et al. [Commun. Math. Sci. 1(4), 715–762 (2003)] is one of the outstanding candidates in projecting micro-scale simulations to a longer timescale. However, its application to tokamak plasma has not been fruitful due to the appearance of spurious transient oscillations in the lifting process, which are present when the kinetic simulations are initialized with a simplified model distribution function and which make the kinetic simulations to deviate from the desired paths. In this work, a kinetically informed lifting algorithm is added to the equation-free projective integration method, which is then verified in the electrostatic gyrokinetic particle-in-cell code XGCa [R. Hager and C. S. Chang, Phys. Plasmas 23, 042503 (2016)] for a neoclassical ion heat transport problem with adiabatic electrons. This new lifting operator is demonstrated to control spurious transients, enabling an over four-times reduction in the overall computing time in the time-evolution of the ion temperature profile in an axisymmetric toroidal plasma. Further reduction in the computing time is found to be limited due to the stability properties of the linear least squares projective integrator.
Extrusion-based AM processes, including material extrusion and directed energy deposition, construct objects by continuously depositing a relatively small amount of molten feedstock to a specific location. Nonextrusion AM processes, on the other hand, tend to form each layer by linearly or areally projecting either energy onto a vat of photocurable liquid or a binding agent into a bed of loose powder. Such systems vary significantly from extrusion-based AM processes but do share similarities. The general approach to path planning is the same, but nonextrusion processes typically favor the use of bitmap representations, instead of polygons, for pathing computation and do not use a g-code file to command the printer. In this chapter, a high-level discussion of these differences will be covered.
Gaussian processes (GPs) are ubiquitously used in sciences and engineering as metamodels. Standard GPs, however, can only handle numerical or quantitative variables. Here we introduce latent map Gaussian processes (LMGPs) that inherit the attractive properties of GPs and are also applicable to mixed data which have both quantitative and qualitative inputs. The core idea behind LMGPs is to learn a continuous, low-dimensional latent space or manifold which encodes all qualitative inputs. To learn this manifold, we first assign a unique prior vector representation to each combination of qualitative inputs. We then use a low-rank linear map to project these priors on a manifold that characterizes the posterior representations. As the posteriors are quantitative, they can be directly used in any standard correlation function such as the Gaussian or Matern. Hence, the optimal map and the corresponding manifold, along with other hyperparameters of the correlation function, can be systematically learned via maximum likelihood estimation. Through a wide range of analytic and real-world examples, we demonstrate the advantages of LMGPs over state-of-the-art methods in terms of accuracy and versatility. In particular, we show that LMGPs can handle variable-length inputs, have an explainable neural network interpretation, and provide insights into how qualitative inputs affect the response or interact with each other. We also employ LMGPs in Bayesian optimization and illustrate that they can discover optimal compound compositions more efficiently than conventional methods that convert compositions to qualitative variables via manual featurization.
Although used in the design and costing of large projects such as linear colliders and fusion tokamaks, the theory of vacuum arcs and gradient limits is not well understood. Almost 120 years after the isolation of vacuum arcs, the exact mechanisms of the arcs and the damage they produce are still being debated. We describe our simple and general model of the vacuum arc that can incorporate all active mechanisms and aims to explain all relevant data. Our four stage model, is based on experiments done at 805 MHz with a variety of cavity geometries, magnetic fields, and experimental techniques as well as data from Atom Probe Tomography and failure analysis of microelectronics. The model considers the trigger, plasma formation, plasma evolution and surface damage phases of the RF arc. This paper also examines how known mechanisms can explain the observed sharp field dependence, fast breakdown times and observed surface damage. We update the model and discuss new features while also pointing out where new data would be useful in extending the model to a wider range of frequencies.
Fermilab’s PIP-II project’s superconducting linear accelerator will drive the next generation of particle accelerators through a revolution in beam intensity. Key to beam intensification are the high efficiencies of niobium-tin superconducting radiofrequency (SRF) cavities operating at cryogenic temperatures near 5 Kelvin. A multifaceted approach is employed to achieve and maintain the extreme temperature. Vacuum provides the first barrier to thermal convection. Physical thermal intercept zones further isolate the exterior shell vacuum vessel at 300 Kelvin from the supercooled beamline. The first thermal intercept lies just inside the exterior vacuum vessel forming the 40 Kelvin zone. This 40K thermal shield is the focus of this investigation. In preparation for operation, the thermal shield is cooled from 300K to 40K in a slow process over 2-1/2 days. The factors limiting the cooling rate are a high thermal gradient which produces thermal strain and high mechanical stress. D ecreasing the time required to reach operational temperatures is the desired research goal. Physical testing data and finite element method (FEM) computer models from the Single Spoke Resonator – 1 (SSR1) and Linac Coherent Light Source (LCLS-II) cryomodules will be evaluated to determine needed design changes. The thermal shield design for the High Beta 650 MHz (HB650) will be modified and analyzed using FEM computer models to determine maximum cooling rate with a target of 50 thermal cycles of life. Maximum allowable stress will be determined through low cycle fatigue calculations.
Integrating biomass-derived intermediates into traditional petroleum refineries presents unique challenges and opportunities, requiring innovative analysis approaches that account for biofuel producers and refiners. Consequently, teams within the National Renewable Energy Laboratory (NREL) have developed a comprehensive refinery integration analysis framework that combines experimental data, detailed techno-economic analyses of bio-conversion pathways, and economic projections within refinery linear programming (LP) optimization models. These models enable the identification of promising bio-integration strategies tailored to specific refinery configurations, economic conditions, and production goals. They also capture upstream and downstream impacts, highlighting critical bottlenecks and research opportunities for increasing the share of biogenic feedstocks in traditional refining operations. Refinery models are also packaged into a broader bio-economy optimization framework which enables biofuel supply chain optimization with standalone biofuel production and refinery co-processing/repurposing options. This presentation discusses promising refinery bio-integration strategies along with key challenges and opportunities identified using NREL's refinery and bio-economy optimization frameworks.
Standard means of representing finite rotation in rigid-body kinematics, including orientation angles, Euler parameters, and Rodrigues parameters, are reviewed and compared. General kinematical relations for a beam theory that treats arbitrarily large rotation are then presented. The standard methods of representing finite rotations are applied to these kinematical expressions, and comparison is made among the standard methods and additional methods found in the literature, such as quasi-coordinates and linear combinations of projection angles. The method of Rodrigues parameters is shown to stand out for both its simplicity and generality when applied to beam kinematics, a result that is really missing from the literature.
This project developed a higher-fidelity model of a recently envisioned small spacecraft propulsion system for precision pointing and proximity control. Plasmonic force propulsion harnesses solar light focused onto plasmon reactive subwavelength nanostructures to accelerate and expel nanoparticle propellant via strong optical forces. The goal of the project was to show that plasmonic space propulsion can provide the level of proximity and attitude control envisioned for future NASA nano/picosatellite missions, a level that is better than state-of-the-art approaches. We achieved this goal by showing that plasmonic force thrusters are feasible for a range of advanced mission concepts requiring swarm formations in a deep space environment. We performed three case studies that evaluated the performance of the plasmonic force propulsion thruster in a deep space, microsatellite swarm formation. These case studies assumed the propulsion system could generate thrust at the level predicted from our Phase 1 study (1.6 μN). Through these cases we were able to analyze the concept within a mission specific context through detailed orbital dynamics calculations. Results indicate that, with the Phase 1 estimated thrust level, the approach is promising for providing attitude control to swarm formation spacecraft. Further, we achieved goals related to technology development. Specifically, we experimentally demonstrated nanoparticle acceleration due to plasmonic forces with asymmetric nanostructures excited by focused laser light. Additionally, we investigated the thrust sensitivity and nanoparticle propellant injection dependencies upon thermal effects. As a result of our study, plasmonic force propulsion is at an early TRL 3. Active research and design has been conducted analytically and in the laboratory. Furthermore, practical applications such as the three case studies have been identified for the scientific basic principles that were observed. Future efforts related to fundamental understanding of these techniques should focus on 1) developing a standalone array of asymmetric nanostructures that can effectively interact with a stream or reservoir of particles or 2) experimentally evaluate a dielectrophoretic injector for nanoparticle propellant. The main limitation discovered about plasmonic propulsion regards performance estimates significantly below the Phase 1 estimations. Specifically, original assumptions in the Phase 1 project (notably, a linear array of asymmetric nanostructures) is not a viable approach to achieving significant acceleration, high exhaust velocity, of nanoparticles. More specifically, we assumed in Phase 1 that nanoparticles would be accelerated in series by a long linear array of asymmetric nanostructures. That is, the acceleration of the nanoparticle would build and increase with the kick received by each subsequent nanostructure. This is fundamentally flawed. The potential profile of a single nanostructure is such that it prohibits this phenomenon. The potential energy associated with the plasmon-generated dielectrophoretic force is a potential well, which is good for trapping nanoparticles, but cannot provide significant acceleration of particles to expel them out and away from the nanostructure. Further, a nanoparticle expelled from the first nanostructure would need to overcome the potential barrier for entry into the next nanostructure accelerating stage. Fundamentally, this effect means that a linear array of nanostructures is not a viable accelerating structure. Correspondingly then, acceleration can, or should, only be provided by one nanostructure, and the net acceleration and thrust force of a single nanostructure is small (~cm/s exhaust velocities, sub-nN level thrust vs. the 100’s m/s, μN originally envisioned). While our experiments demonstrated acceleration and manipulation of a nanoparticle using laser light in aqueous environment, the achievable energy and momentum addition to the nanoparticle from a single nanostructure stage is too low for useful propulsion. In terms of thrust prediction, the estimated thrust of 1.6 μN in Phase 1 is reduced to a few nN of thrust with this new insight and understanding of the concept. This thrust level is too small to achieve attitude control of swarms as originally envisioned.
A new technique for projecting high-dimensional data to low-dimensional spaces, called locally linear embedding (LLE) has recently been introduced.
Propose a method for preparing path entangled states with a definite photon number larger than two that relies on projective measurements.
The Neutron Induced Fission Fragment Tracking Experiment (NIFFTE) collaboration has performed measurements with a fission time projection chamber (fissionTPC) to study the fission process by reconstructing full three-dimensional tracks of fission fragments and other ionizing radiation. The amount of linear momentum imparted to the fissioning nucleus by the incident neutron can be inferred by measuring the opening angle between the fission fragments. Using this measured linear momentum, fission fragment angular distributions can be converted to the center-of-mass frame for anisotropy measurements. Angular anisotropy is an important experimental observable for understanding the quantum mechanical state of the fissioning nucleus and vital to determining detection efficiency for cross section measurements. Neutron linear momentum transfer to fissioning 235 U, 238 U, and 239Pu and fission fragment angular anisotropy of 235 U and 238 U as a function of neutron energies in the range 130 keV–250 MeV are presented.
The author has identified the following significant results. The pattern of linear systems within the project area has been extended into the western foothill belt of the Sierra Nevada. The chief pattern of linear features in the western Sierran foothill belt trends about N. 10 - 15 deg W., but in the vicinity of the Feather River the trend of the features abruptly changes to about N. 50-60 deg W and appears to be contiguous across the Sacramento Valley with a similar system of linear features in the Coast Ranges. The linear features in the Modoc Plateau and Klamath Mt. areas appear unrelated to the systems detected in the Coast Ranges of Sierran foothill belt. Although the change in trend of the Sierran structural features has been previously suggested and the interrelationship of the Klamath Mt. region with the northern Sierra Nevadas has been postulated, the data obtained from the ERTS-1 imagery strengthens these notions and provides for the first time evidence of a direct connection of the structural trends within the alluviated part of the Sacramento Valley. In addition rocks of Pleistocene and Holocene age are offset by some of the linear features seen on ERTS-1 imagery and hence may record the latest episode of geologic deformation in north-central California.
Tomography has made a revolutionary impact on the physical, biological and medical sciences. The mathematical foundation of tomography is to reconstruct a three-dimensional (3D) object from a set of two-dimensional (2D) projections. As the number of projections that can be measured from a sample is usually limited by the tolerable radiation dose and/or the geometric constraint on the tilt range, a main challenge in tomography is to achieve the best possible 3D reconstruction from a limited number of projections with noise. Over the years, a number of tomographic reconstruction methods have been developed including direct inversion, real-space, and Fourier-based iterative algorithms. Here, we report the development of a real-space iterative reconstruction (RESIRE) algorithm for accurate tomographic reconstruction. RESIRE iterates between the update of a reconstructed 3D object and the measured projections using a forward and back projection step. The forward projection step is implemented by the Fourier slice theorem or the Radon transform, and the back projection step by a linear transformation. Our numerical and experimental results demonstrate that RESIRE performs more accurate 3D reconstructions than other existing tomographic algorithms, when there are a limited number of projections with noise. Furthermore, RESIRE can be used to reconstruct the 3D structure of extended objects as demonstrated by the determination of the 3D atomic structure of an amorphous Ta thin film. We expect that RESIRE can be widely employed in the tomography applications in different fields. Finally, to make the method accessible to the general user community, the MATLAB source code of RESIRE and all the simulated and experimental data are available at https://zenodo.org/record/7273314.
Not Available
Abstract. The sea level contribution of the Antarctic ice sheet constitutes a large uncertainty in future sea level projections. Here we apply a linear response theory approach to 16 state-of-the-art ice sheet models to estimate the Antarctic ice sheet contribution from basal ice shelf melting within the 21st century. The purpose of this computation is to estimate the uncertainty of Antarctica's future contribution to global sea level rise that arises from large uncertainty in the oceanic forcing and the associated ice shelf melting. Ice shelf melting is considered to be a major if not the largest perturbation of the ice sheet's flow into the ocean. However, by computing only the sea level contribution in response to ice shelf melting, our study is neglecting a number of processes such as surface-mass-balance-related contributions. In assuming linear response theory, we are able to capture complex temporal responses of the ice sheets, but we neglect any self-dampening or self-amplifying processes. This is particularly relevant in situations in which an instability is dominating the ice loss. The results obtained here are thus relevant, in particular wherever the ice loss is dominated by the forcing as opposed to an internal instability, for example in strong ocean warming scenarios. In order to allow for comparison the methodology was chosen to be exactly the same as in an earlier study (Levermann et al., 2014) but with 16 instead of 5 ice sheet models. We include uncertainty in the atmospheric warming response to carbon emissions (full range of CMIP5 climate model sensitivities), uncertainty in the oceanic transport to the Southern Ocean (obtained from the time-delayed and scaled oceanic subsurface warming in CMIP5 models in relation to the global mean surface warming), and the observed range of responses of basal ice shelf melting to oceanic warming outside the ice shelf cavity. This uncertainty in basal ice shelf melting is then convoluted with the linear response functions of each of the 16 ice sheet models to obtain the ice flow response to the individual global warming path. The model median for the observational period from 1992 to 2017 of the ice loss due to basal ice shelf melting is 10.2 mm, with a likely range between 5.2 and 21.3 mm. For the same period the Antarctic ice sheet lost mass equivalent to 7.4 mm of global sea level rise, with a standard deviation of 3.7 mm (Shepherd et al., 2018) including all processes, especially surface-mass-balance changes. For the unabated warming path, Representative Concentration Pathway 8.5 (RCP8.5), we obtain a median contribution of the Antarctic ice sheet to global mean sea level rise from basal ice shelf melting within the 21st century of 17 cm, with a likely range (66th percentile around the mean) between 9 and 36 cm and a very likely range (90th percentile around the mean) between 6 and 58 cm. For the RCP2.6 warming path, which will keep the global mean temperature below 2 ∘C of global warming and is thus consistent with the Paris Climate Agreement, the procedure yields a median of 13 cm of global mean sea level contribution. The likely range for the RCP2.6 scenario is between 7 and 24 cm, and the very likely range is between 4 and 37 cm. The structural uncertainties in the method do not allow for an interpretation of any higher uncertainty percentiles. We provide projections for the five Antarctic regions and for each model and each scenario separately. The rate of sea level contribution is highest under the RCP8.5 scenario. The maximum within the 21st century of the median value is 4 cm per decade, with a likely range between 2 and 9 cm per decade and a very likely range between 1 and 14 cm per decade.