Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complex computing”

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,405 records · Page 78

A comparison of digital computer programs for the numerical solution of ordinary differential equations

Recently the determination of the best technique for numerically solving systems of ordinary differential equations on a digital computer has received much attention. The use of these formulas in conjunction with a stepsize control developed is explained, and one of the formulas is chosen for comparison with other integration techniques. This comparison of one of the best of Fehlberg's formulas with the different numerical techniques described in previous studies on a variety of test problems clearly shows the superiority of Fehlberg's formula. That is, on each of the test problems, the chosen Fehlberg formula is able to achieve a given accuracy in less computer time than any of the other techniques tested. Also, the computer program for the chosen Fehlberg formula is less complex and easier to use than the computer programs for most of the other techniques. To illustrate the use of the chosen Fehlberg formula, a computer listing of its application to several example problems is included along with a sample of the computer output from these applications.

Ingram, H. L.↗

Three dimensional grid generation for complex configurations: Recent progress

Some of the capabilities of the Computational Fluid Dynamics (CFD) community are surveyed for griding complex three dimensional configurations. The intent is to provide some insight as to the present state of grid generation for aircraft configurations in order to help assess whether this task presents a long term stumbling block to routine use of CFD in aerodynamic applications. A brief review is given of some of the techniques that are available for generating body conforming curvilinear grids. In order to assess capabilities in grid generation, colleagues at selected institutions were solicited to describe their experiences and difficulties in grid generation of complex configurations.

J.F. Thompson↗

Dislocation Patterning in Deforming Crystals: Theory, Computational Predictions and Validation (Final Technical Report)

This project was awarded for an initial period of three years, followed by a no-cost extension for one year and a funded one-year renewal. During the first four years, the project focused on investigating the role of dislocation reactions and dislocation correlation in dislocation patterning in FCC metals. During the fifth year, the scope of research was expanded to investigate the effects of composition inhomogeneity on the mesoscale plastic response of Body-Centered Cubic (BCC) alloys. In addition to its impact on metal hardening during deformation, dislocation patterning provides the microstructure information required to understand phenomena like fracture and recrystallization in metals. The Continuum Dislocation Dynamics (CDD) framework was the methodology used during the first four years, followed by the use of Discrete Dislocation Dynamics (DDD) during the fifth year. CDD is a density-based formalism of dislocation dynamics in which the plastic deformation of crystals is predicted together with the mesoscale dislocation patterns by tracking the space-time evolution of dislocations as driven by the applied stress, short-range and long-range interactions of dislocations, and cross slip. CDD is a crystal mechanics approach in which the plastic constitutive part is replaced with the equations of dislocation dynamics, which is driven by the internal stress via a mobility laws, while the evolution of the density gives the eigenstrain required to update the internal stress itself. The governing equations are thus those of crystal mechanics cast as an eigenstrain problem and those of dislocation transport and reactions. Our investigations during the first four years were driven by the hypothesis that that dislocation patterning is triggered by spatiotemporal dislocation density and internal stress fluctuations, and that such fluctuations influence the collective dislocation dynamic through their effects on the short-range reaction rates and the collective dislocation mobility. This hypothesis was tested by modeling the influence of the dislocation reactions and dislocation correlations on collective dislocation dynamics. The main research components during the first four years were: 1) Reformulation of the CDD framework to integrate the dislocation correlations and dislocation reactions. 2) Simulation of the dislocation correlations and quantifying their contribution to the long-range stress of complex dislocation systems. 3) Computational implementation of the updated CDD model for FCC single crystals and development of an efficient CDD code. 4) Investigation of dislocation patterning in FCC crystals based on the updated CDD model. Some key findings from these investigations were: 1) The patterning of dislocations is initiated by cross slip and internal stress fluctuations, and the refinement of the pattern results from junction formation. 2) Patterning is more prominent under high stress. 3) The correlation stress was found to be a significant part of the mean field stress in continuum representations of dislocation dynamics. During the last year of the project, we have established a method for performing dislocation dynamics in inhomogeneous alloys and demonstrated the impact of composition inhomogeneity on the characteristics of initial yielding and dislocation character in BCC alloys. We also tested this method using an irradiated ferritic alloy in which the composition and dislocation loops effects are both present. In doing so, we have considered two aspects of the role of inhomogeneity, the creation of internal coherency stress due to the dependence of the lattice parameter on the local composition, and the dependence of the dislocation mobility on the local composition. These two aspects resulted in an unexpected behavior of the dislocation system, and, in turn, the yielding behavior of BCC alloys. Key findings include: 1) The composition inhomogeneity in BCC alloys alters the well-known role of screw dislocations by forcing a level of waviness on the average dislocation character, thus destroying the screw character dominance in the yielding observed in pure or homogeneous BCC metals. 2) In the case of irradiated alloys, while composition inhomogeneity by itself and dislocation loops by themselves result in some hardening, the superposition of the two mechanisms does not result in the sum of the two contributions via any known rule. The latter result contradicts the classical works related hardening via multiple mechanisms, which state that various hardening mechanisms can be added linearly or using some Pythagorean additive rule. We were able to rationalize this unexpected result by a closer loop at the origin of the hardening mechanisms themselves, which both involved long-range elastic interactions with dislocations. We reached the conclusion that the hardening resulting from the total stress associated with these two mechanisms (current work) differs from the sum of the effects of the individual stress fields (classical literature). We have also found a major flaw in the classical theory of spinodal hardening. The latter theory never considered the impact of composition undulation on the solute hardening. We have built a new general theory accounting for the effect of solute friction together with the composition undulations on the dislocation configuration at Critical Resolved Shear Stress (CRSS), which is yielding a new definition of the spinodal hardening and new different results. This work is the first to recognize the lack of the role of solute in spinodal hardening theory of alloys.

36 MATERIALS SCIENCE↗

Leveraging Machine Learning Capabilities for the Characterization of Irradiated Uranium: A Case Study of Analysis Methods for Nuclear Safeguards and Nuclear Forensics

Nondestructively determining the initial enrichment of irradiated uranium is a complex and laborious multivariable problem due to the presence of fission products. This work demonstrates the capabilities of machine learning to analyze gamma-ray spectral data to determine initial enrichment without knowledge of the decay time of the sample. The approach developed is agnostic to the particular scenario and is applicable to a wide variety of applications in nuclear forensics and nuclear safeguards. We irradiated 5 mg uranium standard reference materials at discrete enrichment values ranging from 0.02% to 97% 235 U (weight percent) in UT Austin’s Nuclear Engineering Teaching Laboratory TRIGA Mark II 1.1 MW research reactor, allowed each to decay for 8 hours, and then measured each sample via gamma-ray spectrometry for 50 hours post-irradiation yielding 1,400 individual gamma-ray spectra discretized into 8,192 energy bins. We then trained decision trees models to analyze individual gamma-ray spectra and estimate the associated initial enrichment without knowledge of the time since end of irradiation. We evaluated the performance of the models with a reserved test set not used for training or calibrating the model. A decision tree model constructed with this procedure achieved a mean absolute error in initial enrichment determination of 2.3% (weight percent 235 U). Next, we implemented a principal component analysis pre-processing routine of the gamma-ray spectrometry data to reduce the dimensionality of the dataset from 8,192 channels in the spectrum to 10 principal components while retaining over 99% of the inherent variance in the data. Decision tree models constructed with these data demonstrated decreased mean absolute error in enrichment determination, reduced computational time, and decreased complexity. A single decision tree model constructed with this procedure achieved a mean absolute error in initial enrichment determination of 0.05% (weight percent 235 U). Furthermore, we analyzed these models with learning curves to ensure that overfitting did not occur. The capabilities provided by these models can be naturally extended to other application-focused measurements in the fields of nuclear safeguards, nuclear forensics, and nuclear non-proliferation.

Drescher, Adam↗

BIPSPI+: Mining Type-Specific Datasets of Protein Complexes to Improve Protein Binding Site Prediction

Computational approaches for predicting protein-protein interfaces are extremely useful for understanding and modelling the quaternary structure of protein assemblies. In particular, partner-specific binding site prediction methods allow delineating the specific residues that compose the interface of protein complexes. In recent years, new machine learning and other algorithmic approaches have been proposed to solve this problem. However, little effort has been made in finding better training datasets to improve the performance of these methods. With the aim of vindicating the importance of the training set compilation procedure, in this work we present BIPSPI+, a new version of our original server trained on carefully curated datasets that outperforms our original predictor. We show how prediction performance can be improved by selecting specific datasets that better describe particular types of protein interactions and interfaces (e.g. homo/hetero). In addition, our upgraded web server offers a new set of functionalities such as the sequence-structure prediction mode, hetero- or homo-complex specialization and the guided docking tool that allows to compute 3D quaternary structure poses using the predicted interfaces. BIPSPI+ is freely available at https://bipspi.cnb.csic.es.

59 BASIC BIOLOGICAL SCIENCES↗

A Bayesian Approach for Quantifying Data Scarcity when Modeling Human Behavior via Inverse Reinforcement Learning

Computational models that formalize complex human behaviors enable study and understanding of such behaviors. However, collecting behavior data required to estimate the parameters of such models is often tedious and resource intensive. Thus, estimating dataset size as part of data collection planning (also known as Sample Size Determination) is important to reduce the time and effort of behavior data collection while maintaining an accurate estimate of model parameters. In this paper, we present a sample size determination method based on Uncertainty Quantification (UQ) for a specific Inverse Reinforcement Learning (IRL) model of human behavior, in two cases: 1) pre-hoc experiment design—conducted in the planning stage before any data is collected, to guide the estimation of how many samples to collect; and 2) post-hoc dataset analysis—performed after data is collected, to decide if the existing dataset has sufficient samples and whether more data is needed. Here, we validate our approach in experiments with a realistic model of behaviors of people with Multiple Sclerosis (MS) and illustrate how to pick a reasonable sample size target. Our work enables model designers to perform a deeper, principled investigation of effects of dataset size on IRL.

97 MATHEMATICS AND COMPUTING↗

Time Series Dimension Reduction for Surrogate Models of Port Scanning Cyber Emulations

Surrogate model development is a key resource in the scientific modeling community for providing computational expedience when simulating complex systems without loss of great fidelity. The initial step to development of a surrogate model is identification of the primary governing components of the system. Principal component analysis (PCA) is a widely used data science technique that provides inspection of such driving factors, when the objective for modeling is to capture the greatest sources of variance inherent to a dataset. Although an efficient linear dimension reduction tool, PCA makes the fundamental assumption that the data is continuous and normally distributed. Thus, it provides ideal performance when these conditions are met. In the case for which cyber emulations provide realizations of a port scanning scenario, the data to be modeled follows a discrete time series function comprised of monotonically increasing piece-wise constant steps. The sources of variance are related to the timing and magnitude of these steps. Therefore, we consider using XPCA, an extension to PCA for continuous and discrete random variates. This report provides the documentation of the trade-offs between the PCA and XPCA linear dimension reduction algorithms, for the intended purpose to identify key components of greatest variance in our time series data. These components will ultimately provide the basis for future surrogate models of port scanning cyber emulations.

97 MATHEMATICS AND COMPUTING↗

Mixed-mode oscillations in a three-timescale coupled Morris–Lecar system

Mixed-mode oscillations (MMOs) are complex oscillatory behaviors of multiple-timescale dynamical systems in which there is an alternation of large-amplitude and small-amplitude oscillations. It is well known that MMOs in two-timescale systems can arise either from a canard mechanism associated with folded node singularities or a delayed Andronov–Hopf bifurcation (DHB) of the fast subsystem. While MMOs in two-timescale systems have been extensively studied, less is known regarding MMOs emerging in three-timescale systems. In this work, we examine the mechanisms of MMOs in coupled Morris–Lecar neurons with three distinct timescales. We investigate two kinds of MMOs occurring in the presence of a singularity known as canard-delayed-Hopf (CDH) and in cases where CDH is absent. In both cases, we examine how features and mechanisms of MMOs vary with respect to variations in timescales. Our analysis reveals that MMOs supported by CDH demonstrate significantly stronger robustness than those in its absence. Moreover, we show that the mere presence of CDH does not guarantee the occurrence of MMOs. This work yields important insights into conditions under which the two separate mechanisms in two-timescale context, canard and DHB, can interact in a three-timescale setting and produce more robust MMOs, particularly against timescale variations.

97 MATHEMATICS AND COMPUTING↗

How generic is iodide-tagging photoelectron spectroscopy: An extended investigation on the Gly· X - (Gly = glycine, X = Cl or Br) complexes

We report a joint negative ion photoelectron spectroscopy (NIPES) and quantum chemical computational study on glycine-chloride/bromide complexes (denoted Gly·X - , X = Cl/Br) in close comparison to the previously studied Gly·I - cluster ion. Combining experimental NIPE spectra and theoretical calculations, various Gly·X - complexes were found to adopt the same types of low-lying isomers, albeit with different relative energies. Despite more congested spectral profiles for Gly·Cl - and Gly·Br - , spectral assignments were accomplished with the guidance of the knowledge learned from Gly·I - , where a larger spin–orbit splitting of iodine afforded well-resolved, recognizable spectral peaks. Three canonical plus one zwitterionic isomer for Gly·Cl - and four canonical conformers for Gly·Br - were experimentally identified and characterized in contrast to the five canonical ones observed for Gly·I - under similar experimental conditions. Taken together, this study investigates both genericity and variations in binding patterns for the complexes composed of glycine and various halides, demonstrating that iodide-tagging is an effective spectroscopic means to unravel diverse ion-molecule binding motifs for cluster anions with congested spectral bands by substituting the respective ion with iodide.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Correlation methods

A method for processing complex radar data with a computer using correlation functions is reviewed. Parameters including data storage, data reduction, and real time operation are addressed. Since complex auto- and cross-correlation functions are calculated and stored, almost no information is lost. These also can be analyzed in terms of the full correlation analysis of the spaced-antenna-drifts technique. The proposed approach therefore appears to be very feasible to suit most Mesosphere-Stratosphere-Troposphere (MST) radar applications.

Rottger, J.↗

Two-Layer Dynamic Recycling Model (2L-DRM): Learning from Moisture Tracking Models of Different Complexity

Atmospheric moisture tracking models are used to identify and quantify sources and sinks of water in the atmospheric branch of the hydrologic cycle. These models are primarily used to investigate the origin of moisture resulting in precipitation for particular regions around the globe. Moisture tracking models vary widely in their level of complexity, depending on the number of physical processes represented. Complex models are comprehensive in their physical representation, but computationally much more expensive than simple models, which only focus on specific physical processes and use simplifying assumptions. We present the mathematical derivation of the new two-layer dynamical recycling model (2L-DRM), a simple analytical moisture tracking model that relaxes the vertically-integrated formulation of the original one-layer DRM (Dominguez et al., 2006). By comparing the simple DRM to a very complex moisture tracking model that uses water vapor tracers embedded within the Weather Research and Forecast regional climate model (WRF-WVT) for the North American Monsoon region, we pinpoint the absence of vertical wind shear as the main deficiency in the simple DRM. When comparing both simple models (DRM and 2L-DRM) to the WRF-WVT (which we treat as ‘truth’), the 2L-DRM better captures the spatial extent, the net amount and the temporal variability of precipitation that originates from oceanic and local terrestrial sources. The 2L-DRM is well suited to study the large-scale climatological sources of moisture, and for these applications, performs on-par with the much more complex and computationally demanding WRF-WVT model.

Dominguez, Francina↗

Rigorous computation of short-range order unifies its controversial effects in complex concentrated alloys

Direct experimental observations of chemical short-range order (SRO) in complex concentrated alloys (CCAs) have triggered high interest. However, the reported effects of SRO on yield stresses are controversial, and their atomic-scale mechanisms are elusive, which limits our ability to utilize SRO in alloy design. Here we tackle this challenge using an advanced computational approach that rigorously takes into account the critical lattice distortion in CCAs and further verify our theoretical predictions with experiments. We show that the CoCrNi model alloy has a narrow temperature window around 670 °C for SRO formation. This explains why the mechanical effect of SRO is observed in some experiments but not in others. Here, we propose an effective alloy-doping method to control SRO and reveal atomic-bonding types that dominate SRO formation for different alloys. The strategies and insights generally apply to a broad spectrum of alloys, laying the foundation for designing advanced alloys by manipulating their SRO.

36 MATERIALS SCIENCE↗

Coupling a Computational Fluid Dynamics Model to a Spacecraft Thermal System Model for the DraMS Instrument Thermal Analysis

The Dragonfly Mass Spectrometer (DraMS) is an instrument on the Dragonfly mission, which will spend 7 years in deep space cruise before landing and operating on the surface of Titan. Vacuum thermal analyses are required for deep space cruise, and convection analyses are required for the Titan surface operations. Model exchanges across multiple thermal teams are needed for all phases of the mission. For DraMS, Thermal Desktop® (TD) has been the main thermal analytical tool of choice due to its capability in modeling complex thermal systems with relatively low computational power and for its availability across thermal teams. However, TD does not have computational fluid dynamics (CFD) capability and struggles to accurately capture complex convective behavior. DraMS has fans operating in tandem and gas flow behaviors are not easily predicted due to its complex flow paths. CFD software, such as Fluent, can model and predict such complex flow behaviors, but CFD models are computationally expensive, and its workflow processes are not tailored towards simulating large and complex systems. Therefore, a coupled modeling approach was chosen for DraMS: A TD model was used for simulating all the conductive, radiative, and source terms, while a Fluent CFD model was added on, as needed, to the TD model to provide the convective boundary conditions using the System Coupling software. The coupling software allows the TD and Fluent models to communicate data and arrive at a co-solved and co-converged solution. Furthermore, Thermal Iso-value Exchange (TIE) method was developed to facilitate and improve the TD-Fluent data exchange process. This paper will discuss the analytical studies that were done to verify the accuracy and usability of the coupled approach and the challenges associated, which lead to the development of the TIE approach. DraMS thermal design and co-solved analysis results will also be discussed.

Heat transfer↗

Coupling a Computational Fluid Dynamics (CFD) Model to a Spacecraft Thermal System Model for the DraMS Instrument Thermal Analysis

The Dragonfly Mass Spectrometer (DraMS) is an instrument on the Dragonfly mission, which will spend 7 years in deep space cruise before landing and operating on the surface of Titan. Vacuum thermal analyses are required for deep space cruise, and convection analyses are required for the Titan surface operations. Model exchanges across multiple thermal teams are needed for all phases of the mission. For DraMS, Thermal Desktop (TD) has been the main thermal analytical tool of choice due to its capability in modeling complex thermal systems with relatively low computational power and for its availability across thermal teams. However, TD does not have computational fluid dynamics (CFD) capability and struggles to accurately capture complex convective behavior. DraMS has fans operating in tandem and gas flow behaviors are not easily predicted due to its complex flow paths. CFD software, such as Fluent, can model and predict such complex flow behaviors, but CFD models are computationally expensive, and its workflow processes are not tailored towards simulating large and complex systems. Therefore, a coupled modeling approach was chosen for DraMS: A TD model was used for simulating all the conductive, radiative, and source terms, while a Fluent CFD model was added on, as needed, to the TD model to provide the convective boundary conditions using the System Coupling software. The coupling software allows the TD and Fluent models to communicate data and arrive at a co-solved and co-converged solution. Furthermore, Thermal Iso-value Exchange (TIE) method was developed to facilitate and improve the TD-Fluent data exchange process. This paper will discuss the analytical studies that were done to verify the accuracy and usability of the coupled approach and the challenges associated, which lead to the development of the TIE approach. DraMS thermal design and co-solved analysis results will also be discussed.

heat transfer↗

Investigation of Main Bearing Fatigue Estimate Sensitivity to Synthetic Turbulence Models Using a Novel Drivetrain Model Implemented in OpenFAST

ABSTRACT A coupled medium‐fidelity drivetrain model is developed and implemented in OpenFAST for a 10‐MW land‐based reference turbine. The implementation is verified against a fully coupled multibody wind turbine model, including a detailed drivetrain. The new model can simultaneously and accurately estimate main bearing loads and represent elastic bending of the drivetrain. It has low computational cost and is useful for early design phases, sensitivity analyses and complex systems like wind farms (where computational expense must be expended elsewhere). Here, the model is implemented for a monopile offshore wind turbine and used to investigate the sensitivity of main bearing basic rating life to different synthetic turbulence models. Large‐eddy simulations (LES) targeting stable, neutral, and unstable atmospheric conditions at below‐, near‐ and above‐rated wind speeds are used as a reference. The turbulence models recommended by the International Electrotechnical Commission, the Mann spectral tensor model, and the Kaimal spectral model with exponential coherence are fitted to the LES data. Additionally, a constrained turbulence generator, PyConTurb (short for Python Constrained Turbulence ), based on LES data, is applied in the aero‐hydro‐servo‐elastic simulations. Taking PyConTurb as the baseline, the Kaimal model significantly underestimates fatigue of the downwind main bearing, with between 10% and 40% less damage. The Mann model also underestimates the downwind main bearing fatigue by up to 30%. The upwind main bearing damage is driven by mean loads, and differences between models are less significant, although the trends are similar. Reasons for these discrepancies are investigated and attributed to differences in spatial and temporal variations among the turbulence models.

17 WIND ENERGY↗

Negative complexity of formation: the compact dimensions strike back

We show that the vacuum-subtracted maximal volume, the proposed holographic dual to complexity of formation, can be negative when contributions from compact directions are included. We construct explicit solutions with arbitrarily negative complexity of formation in asymptotically AdS 4 × S 7 SUGRA. These examples rely critically on the compact directions, specifically the fact that the full eleven-dimensional spacetime is not asymptotically AdS 11 . While there is some ambiguity in the extension of the holographic complexity proposal to the compact directions, we show that the two natural candidates can both have arbitrarily negative complexity of formation in SUGRA solutions. We further find examples in which complexity can even decrease at late times, including cases of both single-sided geometries and two-sided wormholes. In particular, we construct a cosmological wormhole with simultaneously negative and decreasing complexity of formation (as computed by volume) at late times. We find a distinguished role for relevant primaries in these constructions and comment on possible interpretations.

2D Gravity↗

Embedded symmetric positive semi-definite machine-learned elements for reduced-order modeling in finite-element simulations with application to threaded fasteners

Here, we present a machine-learning strategy for finite element analysis of solid mechanics wherein we replace complex portions of a computational domain with a data-driven surrogate. In the proposed strategy, we decompose a computational domain into an “outer” coarse-scale domain that we resolve using a finite element method (FEM) and an “inner” fine-scale domain. We then develop a machine-learned (ML) model for the impact of the inner domain on the outer domain. In essence, for solid mechanics, our machine-learned surrogate performs static condensation of the inner domain degrees of freedom. This is achieved by learning the map from displacements on the inner-outer domain interface boundary to forces contributed by the inner domain to the outer domain on the same interface boundary. We consider two such mappings, one that directly maps from displacements to forces without constraints, and one that maps from displacements to forces by virtue of learning a symmetric positive semi-definite (SPSD) stiffness matrix. We demonstrate, in a simplified setting, that learning an SPSD stiffness matrix results in a coarse-scale problem that is well-posed with a unique solution. We present numerical experiments on several exemplars, ranging from finite deformations of a cube to finite deformations with contact of a fastener-bushing geometry. We demonstrate that enforcing an SPSD stiffness matrix drastically improves the robustness and accuracy of FEM–ML coupled simulations, and that the resulting methods can accurately characterize out-of-sample loading configurations with significant speedups over the standard FEM simulations.

97 MATHEMATICS AND COMPUTING↗

Clustering and Cliques in Preferential Attachment Random Graphs with Edge Insertion

In this paper, we investigate the global clustering coefficient (a.k.a transitivity) and clique number of graphs generated by a preferential attachment random graph model with an additional feature of allowing edge connections between existing vertices. Specifically, at each time step t, either a new vertex is added with probability f(t), or an edge is added between two existing vertices with probability 1 – f(t). We establish concentration inequalities for the global clustering and clique number of the resulting graphs under the assumption that f(t) is a regularly varying function at infinity with index of regular variation –$\gamma$, where $\gamma$ $\in$ [0, 1). Finally, we also demonstrate an inverse relation between these two statistics: the clique number is essentially the reciprocal of the global clustering coefficient.

97 MATHEMATICS AND COMPUTING↗