Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MATHEMATICAL EVOLUTION”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Complex Dynamics of Air Traffic Flow

Air traffic in the United States has continued to grow at a steady pace since 1980, except for a dip immediately after the tragic events of September 11, 2001. There are different growth scenarios associated both with the magnitude and the composition of the future air traffic. The Terminal Area Forecast (TAF), prepared every year by the FAA, projects the growth of traffic in the United States. Both Boeing and Airbus publish market outlooks for air travel annually. Although predicting the future growth of traffic is difficult, there are two significant trends: heavily congested major airports continue to see an increase in traffic, and the emergence of regional jets and other smaller aircraft with fewer passengers operating directly between non-major airports. The interaction between air traffic demand and the ability of the system to provide the necessary airport and airspace resources can be modeled as a network. The size of the resulting network varies depending on the choice of its nodes. It would be useful to understand the properties of this network to guide future design and development. Many questions, such as the growth of delay with increasing traffic demand and impact of the en route weather on future air traffic, require a systematic understanding of the properties of the air traffic network. There has been a major advance in the understanding of the behavior of networks with a large number of components. Several theories have been advanced about the evolution of large biological and engineering networks by authors in diversified disciplines like physics, mathematics, biology and computer science. Several networks exhibit a scale-free property in the sense that the probabilistic distribution of their nodes as a function of connections decreases slower than an exponential. These networks are characterized by the fact that a small number of components have a disproportionate influence on the performance of the network. Scale-free networks are tolerant to random failure of components, but are vulnerable to selective attack on components. This paper examines two network representations for the baseline air traffic system. A network defined with the 40 major airports as nodes and with standard flight routes as links has a characteristic scale: all nodes have 60 or more links and no node has more than 460 links. Another network is defined with baseline aircraft routing structure exhibits an exponentially truncated scale-free behavior. Its degree ranges from 2 connections to 2900 connections, and 225 nodes have more than 250 connections. Furthermore, those high-degree nodes are homogeneously distributed in the airspace. A consequence of this scale-free behavior is that the random loss of a single node has little impact, but the loss of multiple high-degree nodes (such as occurs during major storms in busy airspace) can adversely impact the system. Two future scenarios of air traffic growth are used to predict the growth of air traffic in the United States. It is shown that a three-times growth in the overall traffic may result in a ten-times impact on the density of traffic in certain parts of the United States.

Scale-free Networks↗
NASA NTRS2001

Mathematical Modeling of Ni/H2 and Li-Ion Batteries

The modelling effort outlined in this viewgraph presentation encompasses the following topics: 1) Electrochemical Deposition of Nickel Hydroxide; 2) Deposition rates of thin films; 3) Impregnation of porous electrodes; 4) Experimental Characterization of Nickel Hydroxide; 5) Diffusion coefficients of protons; 6) Self-discharge rates (i.e., oxygen-evolution kinetics); 7) Hysteresis between charge and discharge; 8) Capacity loss on cycling; 9) Experimental Verification of the Ni/H2 Battery Model; 10) Mathematical Modeling Li-Ion Batteries; 11) Experimental Verification of the Li-Ion Battery Model; 11) Integrated Power System Models for Satellites; and 12) Experimental Verification of Integrated-Systems Model.

Weidner, John W.↗
DOE OSTI2024

Multiscale and Machine Learning Modeling for Process-informed Microstructure Prediction in Additively Manufactured Materials Using MALAMUTE

Advanced Materials and Manufacturing Technologies (AMMT) program under the Department of Energy Office of Nuclear Energy, aims to develop and qualify additively-manufactured materials for nuclear applications. The key challenges to these efforts are the microstructural variabilities observed on the AM products and their impact on the properties and performance of the material in extreme environments. AMMT is using a combination of high-through-put experimental and modeling techniques to accelerate the qualification efforts. Conventionally, in-situ and ex-situ characterizations and testing are performed to correlate different aspects of the AM process to the final product and its performance. However, adopting a trial-and-error approach to experimentally evaluate the vast range of process parameters required to capture the microstructural variabilities is cost-prohibitive. Modeling and simulation provide a comparatively inexpensive way to understand and correlate the microstructural evolution to the processing conditions. The modeling and simulation work-packages within the AMMT program aims to use physics-based and machine learning modeling capabilities to develop a digital twin for AM that can correlate the process conditions to the final product and establish a process-structure-property-performance (PSPP) correlation for AM materials. The melting and subsequent solidification that occurs during the AM process is a complex phenomenon that requires multiscale multiphysics analysis. Idaho National Laboratory’s (INL) Multiphysics Object-Oriented Simulation Environment (MOOSE), specifically the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE) software, provides an ideal platform for developing the multiphysics multiscale model to explore the intricacies of the microstructural evolution during the AM processes within a single framework. Furthermore, given that such full-fidelity simulations can be computationally intensive, reduced order models are necessary to explore the PSPP space for AM materials in an efficient, reliable, and cost-effective way. This work package focuses on understanding the role of process variabilities on the various microstructural characteristics of the AM materials. Microstructures unique to AM materials, such as compositional micro-heterogeneity and dislocation cells, are of particular interest here since they can influence the creep properties and radiation performance. In fiscal year (FY) 24, we significantly advanced upon our work in the last fiscal year, both on physics-based and ML models. The alloy solidification model available in MOOSE has been extended to incorporate the thermodynamic properties and free energy relevant to 316SS. The model demonstrates the Cr segregation that occurs during solidifcation. It is demonstrated that rate of solidification and solute segregation is primarily influence by the cooling rate dictating the level of freezing. This work captures the microstructural variabilities at the subgrain level that are often missing in the part-scale models. With an aim to connect the microstructural evolution model to realistic process conditions, a reduced order model is developed for predicting the thermal conditions around meltpool from high-fidelity process simulations. Furthermore, machine learning approach is used to accelerate the temperature prediction during the AM process. In the following years, MALAMUTE will be used to connect different aspects of the models and quantitatively predict the microstructural evolution. The developed ML-based surrogate model will consider the process conditions as the input to predict the microstructural features in a cost-effective way. The generated microstructures can be used by other work packages under AMMT to evaluate the properties and environmental response of the material at the mesoscale. Thus, this work help identify the key microstructural features at the subgrain level that are significant in property/performance prediction of the AM products. This work will provide inputs to the large-scale process variability models to reevaluate and validate assumptions/simplifications made in the part-scale models. Furthermore, through active learning this work will help identify the data need from both modeling and experimental sides for development of a robust digital twin for AM.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Niobium hydride formation in superconducting qubit thin films

The formation of nonsuperconducting hydrides in 160–170-nm-thick films of niobium is examined. We identify six elastically distinct orientation relationships ɛ−Nb 4 ⁢H 3 takes within the matrix, solid-solution 𝛼−NbH 𝑥 phase. We employ a phase field model to assess the impact of elastic energy induced by the strain of phase transformation on the morphology and transformation dynamics of ɛ precipitates within a thin film. We consider the dimensions of the thin film, crystallographic growth direction, and diffusion rates to predict the timescale of hydride evolution. Leveraging the finite element method, we predict two-dimensional and three-dimensional equilibrium shapes of ɛ-hydrides within a bulk sample and in a thin film that has a traction-free surface. Our results suggest that niobium hydrides migrate to the free surface of the film. Precipitates which reach the free surface coarsen, while precipitates within the film dissolve. Precipitates in both two dimensions and three dimensions experience a repulsive interaction force at the free surface, that is attractive in the bulk, shown in experiment and theory of previous studies.

36 - MATERIALS SCIENCE↗
NASA NTRS1987

Computational astrophysics

Astronomy is an area of applied physics in which unusually beautiful objects challenge the imagination to explain observed phenomena in terms of known laws of physics. It is a field that has stimulated the development of physical laws and of mathematical and computational methods. Current computational applications are discussed in terms of stellar and galactic evolution, galactic dynamics, and particle motions.

Miller, Richard H.↗
DOE OSTI2024

Is the exquisite specificity of lymphocytes generated by thymic selection or due to evolution?

We have previously argued that the antigen receptors of T and B lymphocytes evolved to be sufficiently specific to avoid massive deletion of clonotypes by negative selection. Their optimal ‘specificity’ level, i.e., probability of binding any particular epitope, was shown to be inversely related to the number of self-antigens that the cells have to be tolerant to. Experiments have demonstrated that T lymphocytes also become more specific during negative selection in the thymus, because cells expressing the most crossreactive receptors have the highest likelihood of binding a self-antigen, and hence to be tolerized (i.e., deleted, anergized, or diverted into a regulatory T cell phenotype). Thus, there are two —not mutually exclusive— explanations for the exquisite specificity of T cells, one involving evolution and the other thymic selection. To better understand the impact of both, we extend a previously developed mathematical model by allowing for T cells with very different binding probabilities in the pre-selection repertoire. We confirm that negative selection tends to tolerize the most crossreactive clonotypes. As a result, the average level of specificity in the functional post-selection repertoire depends on the number of self-antigens, even if there is no evolutionary optimization of binding probabilities. However, the evolutionary optimal range of binding probabilities in the pre-selection repertoire also depends on the number of self-antigens. Species with more self antigens need more specific pre-selection repertoires to avoid excessive loss of T cells during thymic selection, and hence mount protective immune responses. We conclude that both evolution and negative selection are responsible for the high level of specificity of lymphocytes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

The effects of introducing elasticity using different interpolation schemes to the grand potential phase field model

Introducing elastic energy in the phase field method has been shown to influence interfacial energy, depending on the elastic interpolation scheme. This study investigates the impact of the elastic energy when using a grand potential-based phase field method, comparing the result of Khachaturyan’s strain interpolation scheme (KHS) and Voight-Taylor’s elastic energy interpolation scheme (VTS). The KHS model leads to a decrease in the interfacial energy, while the VTS model leads to an increase. The change in interfacial energy is greater with the VTS model than the KHS model, which suggests that the KHS model is more appropriate to limit the artificial impact of the elastic energy on the interfacial energy. When the contribution at the interface is not negligible, it is shown that both the microstructure evolution kinetics and the equilibrium microstructure can be influenced by the choice of the elastic scheme being used. In addition, this paper shows that the grand potential model might not be appropriate when the system requires the introduction of a composition-dependent term in the elastic energy contribution. This limitation is due to the need for an explicit and invertible relation between the total potential and the composition.

36 MATERIALS SCIENCE↗
NASA NTRS1985

Vortex motion in wall-bounded viscous flow

A factor of general interest in a broad class of wall-bounded flows is the dynamic evolution of vortical structures through the flow. The structures are three-dimensional, and an overall mathematical description of such entities has not yet been formulated. One of the objectives of the present investigation is concerned with the establishment of a framework, based on first principles, which may form a basis for more detailed analytical studies. Another aim is related to the establishment of boundary and initial conditions in numerical experiments. The mathematical framework employed involves the method of matched asymptotic expansions, and an inner solution field is constructed which consists of a two-dimensional vortical structure. The outer solution field is taken to be an otherwise undisturbed laminar two-dimensional parallel or self-similar viscous flowfield.

Gatski, T. B.↗
NASA NTRS2003

Transition-Sensitized Turbulence Models for Compressible and Incompressible Flows

During the grant period from January 1,2002 to December 31,2002 work was carried out on three projects to extend the range of applicability of advanced turbulence models. First, a new transition-sensitized turbulence model was tested and refined. Second, the influence of compressibility on the pressure-strain rate correlation was studied. Third, the relationship between time-filtered large eddy simulation (TLES) and Reynolds-averaged Navier Stokes (RANS) modeling was investigated leading to submission of the article. The transition-sensitized turbulence model encompasses the early-stage transition and turbulent flow regimes describing the evolution of the ensemble mean disturbance energy and dissipation rate. It is founded on a consistent mathematical description of the laminar regime with its linear disturbances and the fully turbulent regime with its stochastic fluctuations. The unified description is provided by the ensemble viewpoint.

Thacker, William D.↗
NASA NTRS2010

Transforming Mean and Osculating Elements Using Numerical Methods

Mean element propagation of perturbed two body orbits has as its mathematical basis averaging theory of nonlinear dynamical systems. Averaged mean elements define the long-term evolution characteristics of an orbit. Using averaging theory, a near identity transformation can be found that transforms the mean elements back to the osculating elements that contain short period terms in addition to the secular and long period mean elements. The ability to perform the conversion is necessary so that orbit design conducted in mean elements can be converted back into osculating results. In the present work, this near identity transformation is found using the Fast Fourier Transform. An efficient method is found that is capable of recovering the osculating elements to first order

averaging↗
DOE OSTI2022

Q-POP-Thermo: A general-purpose thermodynamics solver for ferroelectric materials

We report that Q-POP-Thermo is a program designed to compute thermodynamic monodomain equilibrium states and their properties for ferroelectric single crystals and thin films based on the Landau-Ginzburg-Devonshire (LGD) Theory. Utilizing symbolic manipulation with the SymPy Library, the governing equations along with appropriate boundary conditions are solved for speedy minimization of the free energy of a crystal. Utilizing the popular Differential Evolution algorithm, with appropriate hybridization, multiple phase diagrams, such as the pressure-temperature phase diagram for bulk single crystals and the common strain-temperature phase diagram for monodomain thin-film systems can be readily generated. Furthermore, a variety of material properties of stable ferroelectric phases, including dielectric, piezoelectric, and electrocaloric properties, can simultaneously be calculated. Validation studies are presented for both thin-film and single crystal systems to test the effectiveness and capability of the open-source program.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Approaches for the Simulation of Coupled Processes in Evolving Fractured Porous Media Enabled by Exascale Computing

Models have historically represented fractured porous media with continuum descriptions that characterize the media using bulk parameters. The impact of small-scale features is not captured in these models, although they may be controlling the performance of subsurface applications. Pore-scale models can simulate processes in small-scale features by representing the pore space geometry explicitly but are computationally expensive for large domains. The alternative multiscale approach entails the combination of pore-scale and continuum-scale descriptions in a single framework. We use Chombo-Crunch, a computational capability that discretizes complex geometries with an adaptive, embedded boundary method to contrast these two approaches. Chombo-Crunch takes advantage of recent computational performance and memory bandwidth improvements resulting from the emergence of exascale computing resources. These combined improvements enable the efficient simulation of reactive transport in fractured media with a high degree of fidelity and the ability to capture the control small-scale processes exert on the overall medium evolution.

42 ENGINEERING↗
DOE OSTI2023

Solving MaxCut with quantum imaginary time evolution

We introduce a method to solve the MaxCut problem efficiently based on quantum imaginary time evolution (QITE). We employ a linear Ansatz for unitary updates and an initial state involving no entanglement, as well as an imaginary-time-dependent Hamiltonian interpolating between a given graph and a subgraph with two edges excised. We apply the method to thousands of randomly selected graphs with up to fifty vertices. We show that our algorithm exhibits a 93% and above performance converging to the maximum solution of the MaxCut problem for all considered graphs. Our results compare favorably with the performance of classical algorithms, such as the greedy and Goemans–Williamson algorithms. We also discuss the overlap of the final state of the QITE algorithm with the ground state as a performance metric, which is a quantum feature not shared by other classical algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Experimental and Theoretical Thermokinetic Analysis of the Na 3 VO 4 –CO 2 Reaction System under Different Physicochemical Conditions

Sodium vanadate (Na 3 VO 4 ) was synthesized, structural and microstructurally characterized as well as tested as possible carbon dioxide (CO 2 ) captor through thermodynamic calculations in addition to thermogravimetric dynamic and isothermal analyses. Ceramic structural characterization evidenced the formation of the Na 3 VO 4 crystal phase, while microstructural features evidenced dense agglomerates with a poor specific surface area, although synthesis temperature (600 °C) was not as high as that for other sodium ceramics. Na 3 VO 4 was investigated for the CO 2 capture process through dynamic and isothermal experiments, in the presence or absence of oxygen. All these experiments showed that Na 3 VO 4 possesses interesting CO 2 capture properties in a specific temperature range (520 and 600 °C). Isothermal products characterization allowed to elucidate the reaction pathway, implying the Na 4 V 2 O 7 formation, as possible reactive intermediate. In fact, reaction path evolution was supported with theoretical thermodynamic data. Moreover, the kinetic analysis, performed using the Avrami-Erofeev mathematical model and subsequent thermodynamic data obtention, probed the positive influence of the temperature, in both CO 2 superficial and bulk chemical sorption processes, while oxygen addition did not enhance the superficial process. Additionally, as it would be expected, kinetics were enhanced as a function of the CO 2 concentration.

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

Control-Oriented Modeling of Cycle-to-Cycle Combustion Variability at the Misfire Limit in SI Engines

A control-oriented model is presented that can capture the prior-cycle correlation of combustion cycles during conditions with high levels of exhaust gas recirculation (EGR). Combustion events are modeled in discrete time and the dynamic evolution is captured by the residual air, fuel, and inert gas trapped in the combustion chamber. The mathematical formulation of the model is presented together with the calibration procedure to emulate a particular engine operating condition. A cycle-to-cycle system identification methodology is described which allows regressing model parameters from experimental data. Simulations are presented and compared to real engine measurements to show the modeling potential for analysis and control of combustion events.

Maldonado Puente, Bryan↗
NASA NTRS2001

Model of statoconia accumulation in gravireceptors of mollusks

The kinetics of formation and accumulation of statoconia are different for Aplysia californica and Biomphalaria glabrata. In Aplysia californica, the fast growth of statoconia number occurs after the critical size (approximately 45 micrometers) of statocyst is reached; then the increase of statoconia number is proceeding with the nonmonotonic rate during the life of an animal. In Biomphalaria the growth of statoconia number occurs only in the initial phase. Then long-term evolution of statoconia in the absence of their generation is the result of their growth in the cyst lumen. In the case of Aplysia californica it is not clear whether a temporal change of the statoconia size distribution (SSD) is caused by statoconia growth in the cyst lumen similar to that in Biomphalaria (Model 1) or statoconia growth takes place in supporting cells until their release into the cyst lumen occurs. (Model 2). This problem is of practical importance because the majority of experiments related to the development of molluscan gravireceptors in altered gravity dealt with an initial phase of statoconia evolution in Aplysia californica and Biomphalaria glabrata. The purpose of the present work is the application of mathematical modeling to the analysis of mechanisms of statoconia formation by supporting cells.

Non-NASA Center↗