Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuous”

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 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

U-splines: Splines over unstructured meshes

U-splines are a novel approach to the construction of a spline basis for representing smooth objects in Computer-Aided Design (CAD) and Computer-Aided Engineering (CAE). A spline is a piecewise-defined function that satisfies continuity constraints between adjacent cells in a mesh. U-splines differ from existing spline constructions, such as Non-Uniform Rational B-splines (NURBS), subdivision surfaces, T-splines, and hierarchical B-splines, in that they can accommodate local variation in cell size, polynomial degree, and smoothness simultaneously over more varied mesh configurations. Mixed cell types (e.g., triangle and quadrilateral cells in the same mesh) and T-junctions are also supported, although the continuity of interfaces with triangle and tetrahedral cells is limited in the present work. The U-spline algorithm introduces a new technique for using local null space solutions to construct basis functions for the global spline null space problem. The U-spline construction is presented for curves, surfaces, and volumes with higher dimensional generalizations possible. Lastly, a set of requirements are given to ensure that the U-spline basis is positive, forms a partition of unity, is complete, and is locally linearly independent.

42 ENGINEERING↗
DOE OSTI2023

Modular machine learning-based elastoplasticity: Generalization in the context of limited data

The development of highly accurate constitutive models for materials that undergo path-dependent processes continues to be a complex challenge in computational solid mechanics. Challenges arise both in considering the appropriate model assumptions and from the viewpoint of data availability, verification, and validation. Recently, data-driven modeling approaches have been proposed that aim to establish stress-evolution laws that avoid user-chosen functional forms by relying on machine learning representations and algorithms. However, these approaches not only require a significant amount of data but also need data that probes the full stress space with a variety of complex loading paths. Furthermore, they rarely enforce all necessary thermodynamic principles as hard constraints. Hence, they are in particular not suitable for low-data or limited-data regimes, where the first arises from the cost of obtaining the data and the latter from the experimental limitations of obtaining labeled data, which is commonly the case in engineering applications. In this work, we discuss a hybrid framework that can work on a variable amount of data by relying on the modularity of the elastoplasticity formulation where each component of the model can be chosen to be either a classical phenomenological or a data-driven model depending on the amount of available information and the complexity of the response. The method is tested on synthetic uniaxial data coming from simulations as well as cyclic experimental data for structural materials. The discovered material models are found to not only interpolate well but also allow for accurate extrapolation in a thermodynamically consistent manner far outside the domain of the training data. This ability to extrapolate from limited data was the main reason for the early and continued success of phenomenological models and the main shortcoming in machine learning-enabled constitutive modeling approaches. Training aspects and details of the implementation of these models into Finite Element simulations are discussed and analyzed.

42 ENGINEERING↗
DOE OSTI2020

Recent developments in catalyst-related PEM fuel cell durability

Cost and durability remain the two major barriers to the widespread commercialization of polymer electrolyte membrane fuel cell (PEMFC)-based power systems, especially for the most impactful but challenging fuel cell electric vehicle (FCEV) application. Commercial FCEVs are now on the road; however, their PEMFC systems do not meet the cost targets established by the U.S. Department of Energy, primarily due to the high platinum loading needed on the cathode to achieve the requisite performance and lifetime. Furthermore, while the activities of a number of commercial Pt-based alloy cathode catalysts exceed the beginning-of-life (BOL) targets, these activities, and the overall cathode performance, degrade via a variety of mechanisms described herein. Degradation is mitigated in current FCEVs by utilizing a cathode catalyst with a lower BOL activity (e.g., much lower transition metal alloy content and larger BOL nanoparticle size), necessitating higher catalyst loadings, and through the utilization of system controls that avoid conditions known to exacerbate degradation processes, such as limiting the fuel cell stack voltage range. The design and development of active and robust materials and eliminating the need for vehicle mitigation strategies would greatly simplify the operating system, allowing for greater transient operation, avoiding large hybridization, and curtailing of fuel cell power. Although system mitigation strategies have provided the near-term pathway for FCEV commercialization, material-specific solutions are required to further reduce costs and improve operability and efficiency. Future material developments should focus on stabilization of the electrode structure and minimization of the catalyst particle susceptibility to dissolution caused by oxide formation and reduction over PEMFC cathode-relevant operating potentials plus minimization of support corrosion. Ex situ accelerated stress tests have provided insight into the processes responsible for material and performance degradation and will continue to provide useful information on the relative stability of materials and benchmarks for robust and stable materials-based solutions not requiring system mitigation strategies to achieve adequate lifetime.

25 ENERGY STORAGE↗
DOE OSTI2026

Flame Chemistry Workshop: a perspective on challenges and strategic actions in combustion experiments and chemical kinetics modeling

Continued progress in the development of predictive models for combustion chemistry—including ignition and flame behavior, species evolution, and combustion system performance—relies on overcoming enduring and emerging challenges in experimental measurements, theoretical formulations, and chemical kinetics mechanism construction. As combustion science continues to coincide with advances in sustainable fuels development, plasma technologies, and automated modeling capabilities, the need for coordinated, community-driven strategies is essential. The Flame Chemistry Workshop (FCWS), held biennially before the International Symposium on Combustion, serves as a dedicated platform to identify, consolidate, and address these challenges in a structured and collaborative manner. This perspective arises from discussions at the 7th FCWS in Milan, Italy (2024), and presents a collective view of the critical barriers currently limiting progress. Across the five technical domains discussed during the 7th FCWS – sustainable fuels combustion, advanced diagnostics for combustion measurements, experiments and modeling in plasma combustion, artificial intelligence and automated methods for theory and mechanisms generation, and chemical kinetic models—a series of persistent and emerging scientific challenges were identified, highlighting the need for deeper integration between three areas: theory, experiments, and modeling. In conclusion, the present article concisely describes present challenges that were identified in each of the technical domains in an effort to streamline and coordinate solutions to accelerate progress in combustion science.

Chemical kinetics↗
DOE OSTI2020

A direct transcription-based multiple shooting formulation for dynamic optimization

The growing need for fast and efficient solution techniques for solving dynamic optimization problems is driven by a broad spectrum of applications in scheduling and control. We suggest a novel framework for dynamic optimization that utilizes a multiple shooting “backbone” with discrete rather than continuous subproblems, thereby eliminating need for repeated time-integration. A Lagrangian relaxation (LR)-based decomposition scheme is proposed, which dualizes the state continuity requirements between subproblems and enables parallel solution of the problem. We demonstrate the applicability of our method on two case studies: the Van der Pol oscillator and a batch reactor.

42 ENGINEERING↗
DOE OSTI2024

Investigation of the notch sensitivity of tailorable long fiber discontinuous prepreg composite laminates

Tailorable discontinuous fiber composite laminates provide relative formability beyond that of continuous fiber laminates, while achieving improved mechanical performance over comparable stochastic systems. Here, in this work, the notch sensitivity of engineered prepreg platelet molded composite (PPMC) laminates is investigated using the open-hole tension (OHT) test and compared to available data for stochastic PPMCs and continuous fiber laminates made with the same material. The press-formed thermoplastic composites (AS4/PEKK) were molded with a quasi-isotropic stacking sequence. The discontinuous PPMC laminate was found to be notch insensitive with OHT strengths ranging from 145.4 MPa (CV $=$ 7%) for d/w $=$ 0.5 to 229.3 MPa (CV $=$ 9%) for d/w $=$ 0.25. The highly ordered meso-structure of the engineered PPMC laminate yields comparatively excellent mechanical properties for relatively thin laminates in contrast to stochastic systems. Both net- and gross-section failures were observed for d/w $=$ 0.25, which suggests that the engineered PPMC laminates studied here maintain a degree of inherent, internal stress concentrations that compete with those caused by geometric features such as a circular hole. Computational simulations of the OHT tests with explicitly represented platelets were found to be in good agreement with experimental measurements. The progressive failure analysis was used to conduct a numerical investigation of the stacking sequence and platelet meso-morphology.

36 MATERIALS SCIENCE↗
DOE OSTI2021

An innovative digital image correlation technique for in-situ process monitoring of composite structures in large scale additive manufacturing

As additive manufacturing (AM) continues to develop and become a standardized manufacturing method, there will be a continued need to provide in-situ monitoring during the manufacturing of polymer composite printed components. Thermal residual stress is a primary cause of failures such as interlayer disbonds or delamination, micro cracking, and dimensional instability, which can occur during or after the build. Here, we report a novel digital image correlation (DIC) adaptation to monitor thermal residual stresses during the entire print process for large-scale AM. In this work, DIC has been investigated (a) by the natural speckle produced by the polymer surface for correlation, (b) to monitor AM build, and (c) to evaluate the effect of thermal residual stress on warpage of the printed component. The natural speckle pattern of the AM material resulted in a respectable 3.57% error compared to the traditional painted speckle pattern of 3.05% error. DIC measured a 190% increase in vertical displacement at the edge of the wall compared to the center, indicating warpage during AM. This work is a step towards a non-intrusive residual stress measuring technique using DIC for large-scale AM.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quasi-Helmholtz decomposition, Gauss' laws and charge conservation for finite element particle-in-cell

Development of particle-in-cell (PIC) methods using finite element based methods (FEMs) have been a topic of renewed interest; this has largely been driven by (a) the ability of finite element methods to better model geometry, (b) better understanding of function spaces that are necessary to represent all Maxwell quantities, and (c) more recently, the fundamental rubrics that should be obeyed in space and time so as to satisfy Gauss' laws and the equation of continuity. In that vein, methods have been developed recently that satisfy these equations and are agnostic to time stepping methods. While this development is indeed a significant advance, it should be noted that implicit FEM transient solvers support an underlying null space that corresponds to a gradient of a scalar potential ∇Φ(r) (or t∇Φ(r) in the case of wave equation solvers). While explicit schemes do not suffer from this drawback, they are only conditionally stable with time step sizes that are mesh dependent and very small. Furthermore, the null space produces spurious charge that can corrupt the desired physics of a PIC simulation. The way to overcome this bottleneck, and indeed, satisfy all four Maxwell's equation is to use a quasi-Helmholtz formulation on a tessellation. In the re-formulation presented, we strictly satisfy the equation of continuity and Gauss' laws for both the electric and magnetic flux densities. Results illustrating the efficacy of this scheme will be demonstrated by analyzing non-neutral beams with space-charge and the adiabatic expansion of a neutral plasma with realistic parameters (Debye length and real mass ratios).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Rethinking the ill-posedness of the spectral function reconstruction — Why is it fundamentally hard and how Artificial Neural Networks can help

Reconstructing hadron spectral functions through Euclidean correlation functions are of the important missions in lattice QCD calculations. However, in a Källen–Lehmann (KL) spectral representation, the reconstruction is observed to be ill-posed in practice. It is usually ascribed to the fewer observation points compared to the number of points in the spectral function. In this paper, by solving the eigenvalue problem of continuous KL convolution, we show analytically that the ill-posedness of the inversion is fundamental and it exists even for continuous correlation functions. We discussed how to introduce regulators to alleviate the predicament, in which include the Artificial Neural Networks (ANNs) representations recently proposed by the Authors in another study. The uniqueness of solutions using ANNs representations is manifested analytically and validated numerically. Reconstructed spectral functions using different regularization schemes are also demonstrated, together with their eigen-mode decomposition. We observe that components with large eigenvalues can be reliably reconstructed by all methods, whereas those with low eigenvalues need to be constrained by regulators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A provably stable numerical method for the anisotropic diffusion equation in confined magnetic fields

We present a novel numerical method for solving the anisotropic diffusion equation in magnetic fields confined to a periodic box which is accurate and provably stable. We derive energy estimates of the solution of the continuous initial boundary value problem. A discrete formulation is presented using operator splitting in time with the summation by parts finite difference approximation of spatial derivatives for the perpendicular diffusion operator. Weak penalty procedures are derived for implementing both boundary conditions and parallel diffusion operator obtained by field line tracing. We prove that the fully-discrete approximation is unconditionally stable. Discrete energy estimates are shown to match the continuous energy estimate given the correct choice of penalty parameters. A nonlinear penalty parameter is shown to provide an effective method for tuning the parallel diffusion penalty and significantly minimises rounding errors. Several numerical experiments, using manufactured solutions, the “NIMROD benchmark” problem and a single island problem, are presented to verify numerical accuracy, convergence, and asymptotic preserving properties of the method. Finally, we present a magnetic field with chaotic regions and islands and show the contours of the anisotropic diffusion equation reproduce key features in the field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

TRANSP integrated modeling code for interpretive and predictive analysis of tokamak plasmas

This paper provides a comprehensive review of the TRANSP code, a sophisticated tool for interpretive and predictive analysis of tokamak plasmas, detailing its major capabilities and features. It describes the equations for particle, power, and momentum balance analysis, as well as the poloidal field diffusion equations. The paper outlines the spatial and time grids used in TRANSP and details the equilibrium assumptions and solvers. Various models for heating and current drive and radiation, including updates to the NUBEAM model, are discussed. The handling of large-scale events such as sawtooth crashes and pellet injections is examined, along with the predictive capabilities for advancing plasma profiles. The integration of TRANSP with the ITER Integrated Modeling and Analysis Suite (IMAS) is highlighted, demonstrating enhanced data access and analysis capabilities. Additionally, the paper discusses best practices and continuous integration techniques to enhance TRANSP's robustness. The suite of TRANSP tools, designed for efficient data analysis and simulation, further supports the optimization of tokamak operations and coupling with other tokamak codes. Continuous development and support ensure that TRANSP remains a major code for the analysis of experimental data for controlled thermonuclear fusion, maintaining its critical role in supporting the optimization of tokamak operations and advancing fusion research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

An editorial to the Special Issue on “Severe climate Risks”

The history of this Special Issue (https://www.sciencedirect.com/special-issue/10JD7LNJNQ0) indirectly dates back to the early 1990s, when the signature of the United Nations Framework Convention on Climate Change kicked-off an international political process based on one overarching and foundational principle: to avoid “dangerous anthropogenic interference with the climate system” at the global level. More than three decades later, such a principle remains central, though complementary aims made their way through the climate negotiation process, such as the importance of ensuring equity and justice, to give just one example here. Scientific knowledge also considerably progressed and we know more about the range of risks that climate change imposes and will continue to impose to the biosphere and humankind, worldwide and at all territorial levels. It is also clear that societal responses to these risks —“climate adaptation” as we know it— are increasingly happening, but definitely not at the pace of climate risk trends (Berrang-Ford et al., 2021, Erisken et al., 2021, Olazabal and Ruiz De Gopegui, 2021, Magnan et al., 2023a, Reckien et al., 2023, UNEP, 2023). As a result, concerns have emerged over the recent years in both the scientific and policy arenas around the idea that societies may not be able to address all climate risks, and that limits to adaptation and induced residual risks need to be considered more seriously. Such concerns further highlight the continuing importance of the imperative to minimise dangerous anthropogenic interference with the climate system, at any scale. But what does “dangerous interference” mean? How can we decide that we are entering the “dangerous” space, compared to a broader range of climate risks that would qualify as problematic but not necessarily “dangerous”? Who should make such a decision? Which conditions drive risk severity over time, including in the future? And what would be the environmental, economic, social and cultural implications of prioritising some climate risks over others? The Intergovernmental Panel on Climate Change (IPCC) was a pioneer in addressing such questions through the development of the “Key Risks” framing that describes those climate risks having the potential to become dangerous or “severe” over the course of this century (Schellnhuber et al., 2006, Schneider et al., 2007, Oppenheimer et al., 2014, O’Neill et al., 2022). The Fifth and Sixth assessment cycles (AR5 and AR6) went a step further by identifying about 120 Key Risks across regions and sectors, and clustering them into 8 “Representative Key Risks” covering a range of geographical systems (low-lying coasts, and to terrestrial and ocean ecosystems), sectors (critical infrastructure, living standards, human health, food security, and water security) and human dimensions (peace and mobility) (Oppenheimer et al., 2014, O’Neill et al., 2022). This Special Issue was born of the efforts of a range of authors, during the development of the IPCC AR6 main Assessment Report between 2019 and 2022, to characterise Key Risks and Representative Key Risks, and advance knowledge on what shapes “severe climate risks” conceptually as well as in the real-world. The series of papers forming this Special Issue is not intended to cover the topic exhaustively, but rather to give readers an overview through the following narrative: defining “severe climate risks” is highly challenging (Magnan et al., 2023b), but knowledge is expanding on the driving climate hazards (Tebaldi et al., 2023) and their implications on geographical systems, sectors and human dimensions, using here food security (Mirzabaev et al., 2023), human mobility (Gilmore et al., 2024) and peace (Buhaug et al., 2023) as illustrative examples. The overall intention is to support especially decision-makers, whatever the scale or sector considered, in asking severity-driven questions to identify sector- and context-specific “priority” risks from climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Optimization-based technoeconomic comparison of multi-stage membrane distillation configurations for hypersaline produced water desalination

Unconventional oil and gas production raises concerns regarding sustainable management of high salinity wastewaters generated in this process. Membrane distillation (MD) is a thermal desalination process capable of treating hypersaline brines such as produced water. The low single pass recovery in MD systems operated in a single stage requires a large recycle stream to achieve the desired recovery, resulting in high energy consumption and operating cost. Multistage configurations in continuous recirculation operation mode offer the potential to reduce the energy intensity of MD systems. However, rigorous analysis is needed to assess the performance of MD configurations when operating in multi-stage mode. Here, we present an optimization-based comparison of economic and energetic performance for five configurations of MD (including DCMD, AGMD, PGMD, CGMD, and VMD) operating in multi-stage continuous recirculation mode. Our findings demonstrate that multistage operation reduces the treatment cost and energy intensity of all MD configurations compared to single stage MD, with the greatest benefit for VMD and PGMD. AGMD with five stages outperforms other configurations with 3.58 $/m 3 feed treatment cost, followed by VMD, CGMD, DCMD, and PGMD with 3.8, 4.2, 5.4, and 9.06$/m 3 feed treatment costs corresponding to twelve, eight, eight, and sixteen stages, respectively.

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

A review of black carbon in snow and ice and its impact on the cryosphere

Black carbon (BC) has emerged as an important short-lived climate forcer. Due to its light absorption properties, BC can darken the snow/ice surface, affect the energy balance, and further lead to acceleration of the cryospheric melting. By reviewing the recent published literatures, we present an overview of the historical changes, spatial distribution of BC in snow/ice, and how these changes are related to the melting of the cryosphere (e.g., glaciers, snow cover, and sea ice). Ice core records show an abrupt increase of BC concentrations that began in the 1850s and continued throughout the 20th century, which is consistent with an increase in BC emissions owing to industrialization. A decrease in the amount of BC since 1970s in Arctic and European ice cores has been partially attributed to the Clean Air Act. However, in the Himalayas, BC records show a continuous increase during this period. BC concentrations in snow and ice in the mid-latitude regions are one to two orders of magnitude higher than those in the polar regions. In particular, BC concentrations in aged snow and granular ice in ablation areas of mountain glaciers are one to three orders of magnitude higher than those in fresh snow or snowpits in the glacier accumulation areas due to BC accumulation during melting season. BC in the surface snow/ice is responsible for about 20% of the albedo reduction in the Tibetan Plateau during glacier melt season. Globally, observations and modeling results indicate that radiative forcing (RF) induced by BC in snow and ice is highest in the mid-latitudes, ranging from several W m–2 in fresh snow to hundreds of W m–2 in aged snow and granular ice in glacier ablation areas. The large BC-in-snow RF and associated snow albedo feedback lead to an acceleration in the total glacier melt (approximately 20%) and/or a reduction in the duration of the snow cover by several days, resulting in an increase of glacier discharge. Given our limited understanding of the role of BC in cryospheric melting, it is important to synthesize the existing research on the multi-scale processes related to BC in snow and ice to identify the gaps in our understanding of these processes and to propose a path forward to improve the quality of our observations of the aforementioned phenomena to fill these gaps.

Kang, Shichang↗
DOE OSTI2026

A U.S. scientific community review of carbon cycle science gaps and opportunities to better support earth system science and carbon management

Greenhouse gas (GHG) emissions continue to grow, while natural carbon reservoirs are becoming increasingly vulnerable to anthropogenic pressures, climate extremes, and disturbance. These changes are impacting humans, ecosystems, and natural resources worldwide. Tracking and mitigating GHG emissions require a pivot to operational monitoring of regional carbon flux and stock changes. The current GHG observing system is addressing needs at two distinct scales: 1) Local scale (< 1 km), related to anthropogenic point source emissions, and 2) global scales (> 1000 km), related to land and ocean carbon sinks. More focus on intermediate (10–1000 km) scales is needed to more effectively monitor progress in reducing carbon emissions, enhancing removals, and maintaining sinks. Representatives from carbon cycle biomass and flux communities across United States government agencies and academic institutions met in September 2024 to discuss the rationale and scientific context for more effectively implementing an operational system for GHG monitoring in support of urban and national carbon management needs. To guide development of this system, we propose a multi-tiered global spaceborne observing framework for carbon flux and stock, prioritizing: 1) frequent GHG partial columns for carbon emissions and removals; 2) continuous time series and data fusion of biomass from Lidar and Synthetic Aperture Radar (SAR) for carbon stocks, and 3) expanded coverage of tropical, high latitude, and oceanic regions to monitor carbon cycle tipping points and feedbacks. This system should be complemented by expanded surface and airborne networks for oceanic and terrestrial/aquatic ecosystems for calibration, ground truthing, and study of under-sampled regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Approximating a linear multiplicative objective in watershed management optimization

Implementing management practices in a cost-efficient manner is critical for regional efforts to reduce the amount of pollutants entering the Chesapeake Bay. We study the problem of selecting a subset of practices that minimizes pollutant load—subject to budgetary and environmental constraints—as simulated in a widely used regulatory watershed model. Mimicking the computation of pollutant load in the regulatory model, we formulate this problem as a continuous optimization model with a linear multiplicative objective function and linear constraints. To lay the groundwork for incorporating additional stakeholder requirements in the future, especially those that would require integer variables, we present and study a continuous linear optimization model that approximates the nonlinear model. The linear model, which requires an exponential number of variables, arises naturally as an alternative model for the same underlying physical process. We examine the theoretical behavior of these optimization models and investigate restrictions of the linear model to handle its large number of variables. Through extensive computational tests on real and randomly generated instances, we demonstrate that the linear model and its restrictions provide optimal solutions close to those of the nonlinear model in practice, despite poor approximation properties in the worst case. We conclude that the linear model—together with our approach to handling its large number of variables—provides a viable framework from which to extend the optimization model to better meet the needs of the Chesapeake Bay watershed management stakeholders.

54 ENVIRONMENTAL SCIENCES↗