Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Relativity”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Numerical Integration of Viscoelastic Models

Accurate modeling of viscoelasticity remains an important consideration for a variety of materials (e.g. polymers and inorganic glasses). As such, over the previous decades a substantial body of work has been dedicated to developing appropriate constitutive models for viscoelasticity ranging from initial considerations of linear thermoviscoelasticity to more complex non-linear formulations incorporating fictive temperatures or potential energy clocks including the use of both internal state variable(ISV) and hereditary integral representations. Nonetheless, relatively limited (in comparison to plasticity) attention has been paid to the numerical integration of such schemes. In terms of integral based formulations, Taylor et al. first considered the problem of the integration of a linear viscoelasticity model. That work focused on the integration of the hereditary integrals and demonstrated improved performance of the new scheme with a custom finite element code over an existing finite difference reference. Chambers and Becker, using a free volume based shift factor, also considered the integration of the hereditary integrals and the impact on the problem of a pressurized thick-walled cylinder and developed an adaptive scheme to bound the error. Chambers later developed three-point Gauss and composite integration schemes for the hereditary integrals and noted improved accuracy. With respect to ISV-based schemes, formulations for the non-linear Schapery model have been proposed. However, in those efforts greater attention was paid to convergence of the non-linear solution scheme than impact of numerical integration. Various authors (e.g. Holzapfel and Simo and Hughes) have also studied the use of convolution integrals with differential forms of ISVs for temperature-independent formulations. Regardless, while the "potential energy clock" (PEC) and "simplified potential energy clock"(SPEC) models have been used to study a variety of non-linear responses (e.g.), limited attention has been paid to the numerical performance. As will be discussed later, the "clock" at the center of the formulations includes temperature and complex history dependence making the numerical integration of such a model even more challenging. Thus, in the current work an initial effort towards characterizing the numerical integration of the constitutive model through simplified problems is performed. To that end, in Section 2 the theory of the model is briefly presented while the numerical integration is discussed in Section 3. Results of various studies characterizing the numerical behavior and performance are then given in Section 4. Finally, some concluding remarks and thoughts for follow on works are provided in Section 5.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Wave-powered water pump for upwelling in aquaculture: Numerical model and ocean test

Wave-powered upwelling can increase the productivity and survivability of several aquaculture species. This enhancement is due to transporting cold, nutrient-rich ocean water, typically found lower in the water column, to the surface. Macroalgaes, like kelp, exhibit increased growth from these altered conditions. The University of New Hampshire’s (UNH) wave-powered water pump (wave pump) is a point absorber wave energy converter (WEC) that uses ocean waves to create relative motion between a spar buoy and a concentric float which drives an internal pump. A numerical model of the wave pump was developed using WEC-Sim to predict device performance in the ocean. Wave pump performance was evaluated during a five day ocean test near Appledore Island in Maine in March 2023, where volumetric flow rate, relative distance between spar and float, and wave conditions were measured. These data were then used for numerical model validation. The ocean deployment recorded the device’s performance in a variety of sea states, with average significant wave heights up to 0.7 m. The ocean test data were compared to the WEC-Sim numerical model of the device with favorable results. Average values of device stroke period, stroke height, and flow rate agreed between the ocean test and model data to within approximately 16 to 22%. Furthermore, the validated numerical model provides a valuable tool for improving the design and developing a commercial-scale, wave-powered water pump for use in aquaculture.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Unexpected Scaling of Peak Acceleration for a Yielding Mass-Spring System Subjected to a Triangular Base Acceleration Pulse

The use of bounding scenarios is a common practice that greatly simplifies the design and qualification of structures. However, this approach implicitly assumes that the quantities of interest increase monotonically with the input to the structure, which is not necessarily true for nonlinear structures. This paper surveys the literature for observations of nonmonotonic behavior of nonlinear systems and finds such observations in both the earthquake engineering and applied mechanics literature. Numerical simulations of a single degree-of-freedom mass-spring system with an elastic–plastic spring subjected to a triangular base acceleration pulse are then presented, and it is shown that the relative acceleration of this system scales nonmonotonically with the input magnitude in some cases. The equation of motion for this system is solved symbolically and an approximate expression for the relative acceleration is developed, which qualitatively agrees with the nonmonotonic behavior seen in the numerical results. The nonmonotonicity is investigated and found to be a result of dynamics excited by the discontinuous derivative of the base acceleration pulse, the magnitude of which scales nonmonotonically with the input magnitude due to the fact that the first yield of the spring occurs earlier as the input magnitude is increased. The relevance of this finding within the context of defining bounding scenarios is discussed, and it is recommended that modeling be used to perform a survey of the full range of possible inputs prior to defining bounding scenarios.

42 ENGINEERING↗
DOE OSTI2022

Soils and sediments host Thermoplasmata archaea encoding novel copper membrane monooxygenases (CuMMOs)

Copper membrane monooxygenases (CuMMOs) play critical roles in the global carbon and nitrogen cycles. Organisms harboring these enzymes perform the first, and rate limiting, step in aerobic oxidation of ammonia, methane, or other simple hydrocarbons. Within archaea, only organisms in the order Nitrososphaerales (Thaumarchaeota) encode CuMMOs, which function exclusively as ammonia monooxygenases. From grassland and hillslope soils and aquifer sediments, we identified 20 genomes from distinct archaeal species encoding divergent CuMMO sequences. These archaea are phylogenetically clustered in a previously unnamed Thermoplasmatota order, herein named the Ca. Angelarchaeales. The CuMMO proteins in Ca. Angelarchaeales are more similar in structure to those in Nitrososphaerales than those of bacteria, and contain all functional residues required for general monooxygenase activity. Ca. Angelarchaeales genomes are significantly enriched in blue copper proteins (BCPs) relative to sibling lineages, including plastocyanin-like electron carriers and divergent nitrite reductase-like (nirk) 2-domain cupredoxin proteins co-located with electron transport machinery. Ca. Angelarchaeales also encode significant capacity for peptide/amino acid uptake and degradation and share numerous electron transport mechanisms with the Nitrososphaerales. Ca. Angelarchaeales are detected at high relative abundance in some of the environments where their genomes originated from. While the exact substrate specificities of the novel CuMMOs identified here have yet to be determined, activity on ammonia is possible given their metabolic and ecological context. The identification of an archaeal CuMMO outside of the Nitrososphaerales significantly expands the known diversity of CuMMO enzymes in archaea and suggests previously unaccounted organisms contribute to critical global nitrogen and/or carbon cycling functions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Bayesian merging of numerical modeling and remote sensing for saltwater intrusion quantification in the Vietnamese Mekong Delta

Saltwater intrusion has become one of the most concerning issues in the Vietnamese Mekong Delta (VMD) due to its increasing impacts on agriculture and food security of Vietnam. Reliable estimation of salinity plays a crucial role to mitigate the impacts of saltwater intrusion. Here, this study developed a hybrid technique that merges satellite imagery with numerical simulations to improve the estimation of salinity in the VMD. The salinity derived from Landsat images and by numerical simulations was fused using the Bayesian inference technique. The results indicate that our technique significantly reduces the uncertainties and improves the accuracy of salinity estimates. The Nash–Sutcliffe coefficient is 0.74, which is much higher than that of numerical simulation (0.63) and Landsat estimation (0.6). The correlation coefficient between the ensemble and measured salinity is relatively high (0.88). The variance of the ensemble salinity errors (5.0 ppt 2 ) is lower than that of Landsat estimation (10.4 ppt 2 ) and numerical simulations (9.6 ppt 2 ). The proposed approach shows a great potential to combine multiple data sources of a variable of interest to improve its accuracy and reliability wherever these data are available.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

An asynchronous parallel high-throughput model calibration framework for crystal plasticity finite element constitutive models

Crystal plasticity finite element model (CPFEM) is a powerful numerical simulation in the integrated computational materials engineering toolboxes that relates microstructures to homogenized materials properties and establishes the structure–property linkages in computational materials science. However, to establish the predictive capability, one needs to calibrate the underlying constitutive model, verify the solution and validate the model prediction against experimental data. Bayesian optimization (BO) has stood out as a gradient-free efficient global optimization algorithm that is capable of calibrating constitutive models for CPFEM. Here in this paper, we apply a recently developed asynchronous parallel constrained BO algorithm to calibrate phenomenological constitutive models for stainless steel 304 L, Tantalum, and Cantor high-entropy alloy.

304L stainless steel↗
DOE OSTI2021

Improved accuracy and transferability of molecular-orbital-based machine learning: Organics, transition-metal complexes, non-covalent interactions, and transition states

We report molecular-orbital-based machine learning (MOB-ML) provides a general framework for the prediction of accurate correlation energies at the cost of obtaining molecular orbitals. The application of Nesbet’s theorem makes it possible to recast a typical extrapolation task, training on correlation energies for small molecules and predicting correlation energies for large molecules, into an interpolation task based on the properties of orbital pairs. We demonstrate the importance of preserving physical constraints, including invariance conditions and size consistency, when generating the input for the machine learning model. Numerical improvements are demonstrated for different datasets covering total and relative energies for thermally accessible organic and transition-metal containing molecules, non-covalent interactions, and transition-state energies. MOB-ML requires training data from only 1% of the QM7b-T dataset (i.e., only 70 organic molecules with seven and fewer heavy atoms) to predict the total energy of the remaining 99% of this dataset with sub-kcal/mol accuracy. This MOB-ML model is significantly more accurate than other methods when transferred to a dataset comprising of 13 heavy atom molecules, exhibiting no loss of accuracy on a size intensive (i.e., per-electron) basis. It is shown that MOB-ML also works well for extrapolating to transition-state structures, predicting the barrier region for malonaldehyde intramolecular proton-transfer to within 0.35 kcal/mol when only trained on reactant/product-like structures. Finally, the use of the Gaussian process variance enables an active learning strategy for extending the MOB-ML model to new regions of chemical space with minimal effort. We demonstrate this active learning strategy by extending a QM7b-T model to describe non-covalent interactions in the protein backbone–backbone interaction dataset to an accuracy of 0.28 kcal/mol.

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

Diagnostics for multiple frequency heating and investigation of underlying processes

The development of new facilities routinely challenges ion source designers to build and operate sources that can achieve ever higher beam intensities and energies. As shown in this work, electron cyclotron resonance ion sources have proven to be extremely capable in meeting these challenges through the production of intense beams of medium and high-charge state ions. As performance boundaries are pushed, source stability becomes an issue as does the technology required to meet the challenge. Multiple frequency heating, the simultaneous use of two or more plasma heating frequencies, is a powerful tool in meeting the simultaneous need of intensity and stability. Relatively straightforward to utilize, the technique has been employed at numerous facilities to increase beam current and achievable charge state while also stabilizing the plasma. Its application has expanded the operational boundaries of existing and next generation sources, demonstrating that these devices have not yet achieved their full operational potential. To better understand the underlying physics, the diagnostics used to probe the source operational boundaries and the plasma properties have become increasingly sophisticated. In concert with detailed modeling, they are beginning to provide insight into the heating mechanism and, with that, the prospect of future advances.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Connecting primordial gravitational waves and dark energy

Cosmic acceleration manifested in the early universe as inflation, generating primordial gravitational waves detectable in the cosmic microwave background (CMB) radiation. Cosmic acceleration is occurring again at present as dark energy, detectable in cosmic distance and structure surveys. We explore the intriguing idea of connecting the two occurrences through quintessential inflation by an α-attractor potential without a cosmological constant. For this model we demonstrate robustness of the connection 1 + w 0 ≈ 4/(3N 2 r) between the present day dark energy equation of state parameter w 0 and the primordial tensor to scalar ratio r for a wide range of initial conditions. Analytic and numerical solutions produce current thawing behavior, resulting in a tight relation w a ≈ -1.53(1 + w 0 )≈ -0.2 (4 × 10 -3 /r). Upcoming CMB and galaxy redshift surveys can test this consistency condition. In conclusion, within this model, lack of detection of a dark energy deviation from Λ predicts a higher r, and lack of detection of r predicts greater dark energy dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Newton versus the machine: solving the chaotic three-body problem using deep neural networks

ABSTRACT Since its formulation by Sir Isaac Newton, the problem of solving the equations of motion for three bodies under their own gravitational force has remained practically unsolved. Currently, the solution for a given initialization can only be found by performing laborious iterative calculations that have unpredictable and potentially infinite computational cost, due to the system’s chaotic nature. We show that an ensemble of converged solutions for the planar chaotic three-body problem obtained using an arbitrarily precise numerical integrator can be used to train a deep artificial neural network (ANN) that, over a bounded time interval, provides accurate solutions at a fixed computational cost and up to 100 million times faster than the numerical integrator. In addition, we demonstrate the importance of training an ANN using converged solutions from an arbitrary precise integrator, relative to solutions computed by a conventional fixed precision integrator, which can introduce errors in the training data, due to numerical round-off and time discretization, that are learned by the ANN. Our results provide evidence that, for computationally challenging regions of phase space, a trained ANN can replace existing numerical solvers, enabling fast and scalable simulations of many-body systems to shed light on outstanding phenomena such as the formation of black hole binary systems or the origin of the core collapse in dense star clusters.

Breen, Philip G.↗
DOE OSTI2024

A Phenomenological Thermal Spray Wall Interaction Modeling Framework Applied to a High Temperature Ignition Assistant Device

Airborne compression ignition engines must operate with reliable ignition systems to achieve proper ignition at every cycle, particularly at high altitudes. Glow-plug-based ignition-assistant (IA) devices can provide the necessary energy to preheat the fuel and ensure ignitability of the fuel-air mixture. Ignitability of liquid sprays can be facilitated via direct impingement onto the hot IA surface, however this comes with adverse effects on the IA durability. Therefore, optimizing an IA’s design requires detailed understanding of the physics of fuel spray impingement of superheated surfaces. While spray impingement on relatively low wall temperatures has been extensively studied and appropriate numerical models have been proposed through the years, fundamental understanding of high-speed liquid spray impingement on superheated walls is still elusive. This work aims to formulate a phenomenological thermal spray-wall interaction framework for modeling the film-boiling-induced heat transfer, atomization, and dispersion of fuel spray droplets impinging onto a superheated IA device. A qualitative comparison of the new phenomenological model is performed against optical experiments from the literature of an F-24 fuel spray injected onto an IA device located 12 mm away from the injector tip. The temperature of the IA was set at 1400 K. The fuel injection pressure was 400 bar, while the ambient gas pressure and temperature were 30 bar and 800 K, respectively. The performance of the phenomenological model is evaluated in comparison with two other state-of-art models from the literature. A qualitative analysis of the different spray and fuel-air mixture characteristics is performed to outline the differences in the predictions offered by the new phenomenological model and the two state-of-art spray-wall interaction models.

droplet dispersion↗
DOE OSTI2025

Microstructure Scale Lithium-Ion Battery Modeling: Part I. On Through-Plane Heterogeneity, Impact of Mesh Representation, and Differences between Macro- and Microscale Models

Li-ion battery performance and degradation are strongly correlated with the electrode microstructures and can be modeled at different scales, each with their own limitations. Herein, we compare predictions achieved with a macro- and a micro-scale model, that is, respectively, neglecting or considering the microstructural heterogeneity of the composite electrodes, on virtual numerically generated and real microstructures. While both models are in relative agreement at the low charge rates, differences arise for fast charging scenarios and especially for the real, highly heterogenous, microstructures. The microscale model predicts that electrolyte concentration saturation and depletion, respectively, at the back of the cathode and of the anode are exacerbated, and that lithium plating occurs earlier for real microstructures. The present work also indicates that the mesh representation significantly impacts the microscale model predictions, and consequently that microscale models should add surface area as a parameter to consider explicitly surface roughness. This article is the first of a series, with subsequent entries further investigating in-plane heterogeneities, lithium plating, and the impact of microstructure representativity on model predictions.

25 ENERGY STORAGE↗
DOE OSTI2019

Nadir Pressure Exposure Estimates for the Existing Turbines at McNary Dam

This report documents a numerical modeling study conducted to assess the pressure environment and relative biological performance of the existing turbine runner design at the U.S. Army Corps of Engineers (USACE) McNary Dam to establish a baseline for comparison to planned replacement runner designs. Computational fluid dynamics (CFD) models of the turbine unit constructed for this study were analyzed using a streamtrace-based technique called the Biological Performance Assessment (BioPA), which was developed at the Pacific Northwest National Laboratory (PNNL).

42 ENGINEERING↗
DOE OSTI2025

Ab initio equation of state for vanadium

Accurate materials’ equations of state (EOS) are essential for understanding materials properties as well as for use in multiphysics simulations. In particular, hydrodynamics simulations are based on three fundamental conservation laws (mass, momentum and energy) that form an under-determined system of equations. The equation of state serves as an additional closure relation between thermodynamic variables for a given material that enables numerical hydrodynamics simulation. In this report, we focus on the development of an ab initio EOS for the body centered cubic (BCC) phase of Vanadium (V) for eventual integration into a multiphase EOS in the OpenSesame EOS database.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Tidal Origin for a Three-body Resonance in Kepler-221

Over the course of the last two decades, traditional models of planet formation have been repeatedly challenged by the emerging census of extrasolar planets. Key among them is the orbital architecture problem: while standard models of orbital migration predict resonant orbits for short-period objects, most planets do not appear to lie in orbital resonances. Here, we show that the four-planet system Kepler-221, not previously recognized to have active orbital resonances, has a three-body commensurability relation unique within the Kepler sample. Using a suite of numerical experiments as well as a perturbative analysis, we demonstrate that this system likely began as a resonant chain and proceeded to undergo large-scale divergence away from resonance, under the action of tidal dissipation. Our results further indicate that obliquity tides, driven by a secular spin–orbit resonance and mutual inclination, are an excellent candidate for driving this orbital divergence, and that the high tidal luminosity may also explain the anomalous size of planet b, which lies within the Fulton radius gap.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Numerical study of heat transfer enhancement in laminar and turbulent flow in enhanced tubes with uniform and non-uniform thermal boundary conditions

Various types of tube enhancements have been proposed in order to increase heat transfer while minimizing any associated increases in pressure drop. This paper presents a numerical study of two different types of enhanced tubes: the relatively common twisted tape insert tube, and a novel tube with airfoil shaped pins on the inner surface. Each tube is studied in both the laminar and turbulent flow regime. In addition to the typical constant temperature and constant flux boundary conditions, a non-constant heat flux boundary condition is applied on one half of the tube, a situation which may be encountered when absorbing solar radiation. In general, each type of enhanced tube offers a thermal performance factor greater than unity at certain Re. In addition, the non-uniform heat flux boundary condition appears to function similarly as the uniform heat flux boundary condition.

Sandlin, Matthew J.↗
DOE OSTI2020

Expedited Walkdowns and Preparation of Cost Estimates for Decontamination and Demolition of Excess Facilities - 20335

Lawrence Livermore National Security, LLC (LLNS) is the prime contractor responsible for managing and operating Lawrence Livermore National Laboratory (LLNL). The mission of the LLNL is to strengthen security of the United States through development and application of world-class science and technology to enhance the nation's defense, reduce the global threat from terrorism and weapons of mass destruction, and respond with vision, quality, integrity, and technical excellence to scientific issues of national importance. To accomplish its mission, LLNL must plan and manage its campus space and facilities to optimize use of its relatively small one-square-mile footprint. LLNL currently has numerous excess facilities that require decontamination and demolition (D and D). The LLNL Legacy Facility Program is responsible for stewardship and risk reduction programs to effectively manage these excess assets at LLNL. The Legacy Facility Program identified 10 excess facilities at LLNL that were candidates for demolition based on risks from further degradation of the facilities, potential exposure to hazardous materials, or need for removal to make way for new program facilities. The buildings range from a structure built in 1943 to serve as a US Navy drill hall during World War II, to a facility that provided neutrons to study material properties for the fusion energy program. The National Nuclear Security Administration (NNSA) appropriated FY19 funding for the development of Class 3 estimates for these 10 excess facilities that would be used to develop a Program Management Plan for funding consideration. Class 3 estimates are not conceptual, nor overly detailed, but adequate to be used for budget and appropriation purposes. Using existing strategic sourcing agreements, LLNL was able to organize a team of senior subject matter experts (SMEs) with diverse technical backgrounds across the Department of Energy (DOE) Complex to complete extensive process knowledge (PK) reviews and walk-downs of each excess facility and to deliver accurate and complete estimates under an expedited schedule (i.e., within eight weeks). To facilitate the expedited schedule, LLNL was fully prepared for the arrival of the JGMS team, which consisted of the following companies: J.G. Management Systems, Inc. (JGMS); Strata-G, LLC; and Michael Baker International. LLNL provided background and PK information on the excess facilities, pictures, supporting historical documentation, and radiological survey history as well as escorted access to the facilities. The JGMS team has previous experience supporting Program Management Plans for excess facilities at the Y-12 National Security Complex. The team mobilized to LLNL and performed the facility walk-downs from January 22-25, 2019. Using a methodical evaluation process, expertise, and experience, the team was able to deliver walk-down draft reports and estimates to LLNL before their February 8, 2019 deadline. LLNL then in turn reviewed, completed the estimates by adding waste management costs and other LLNL costs, and delivered the estimates to NNSA by February 15, 2019. The completed final reports, including draft D and D estimates, were delivered to LLNL before their March 6, 2019 deadline. The work was successfully completed ahead of the expedited schedule required by LLNL, resulting in a commendation from the client. This paper describes how the partnership and collaboration as a fully integrated team comprised of the prime contractor, LLNS, and the JGMS team successfully met the needs of the US government, resulting in the performance of expedited site walk-downs, the development of summary evaluation reports, and the preparation of draft Class 3 estimates and preliminary schedules for the D and D of 10 excess facilities at LLNL. The efforts of the team will enable the D and D of the first excess facility in the 2020-2021 time frame. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗