Search NASA⌕ Search

SEARCH · Search NASA

Results for “engineering, methods, limits, mathematics”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Adjoint-based sensitivity analysis for high-energy density radiative transfer using flux-limited diffusion

Uncertainty quantification and sensitivity analyses are a vital component for predictive modeling in the sciences and engineering. The adjoint approach to sensitivity analysis requires solving a primary system of equations and a mathematically related set of adjoint equations. The information contained in the equations can be combined to produce sensitivity information in a computationally efficient manner. Here in this work, sensitivity analyses are performed on systems described by flux-limited radiative diffusion using the adjoint approach. The sensitivities computed are shown to agree with standard perturbation theory and require significantly less computational time. The adjoint approach saves the computational cost of one forward solve per sensitivity, making the method attractive when multiple sensitivities are of interest.

42 ENGINEERING↗
NASA NTRS1997

Automated Reuse of Scientific Subroutine Libraries through Deductive Synthesis

Systematic software construction offers the potential of elevating software engineering from an art-form to an engineering discipline. The desired result is more predictable software development leading to better quality and more maintainable software. However, the overhead costs associated with the formalisms, mathematics, and methods of systematic software construction have largely precluded their adoption in real-world software development. In fact, many mainstream software development organizations, such as Microsoft, still maintain a predominantly oral culture for software development projects; which is far removed from a formalism-based culture for software development. An exception is the limited domain of safety-critical software, where the high-assuiance inherent in systematic software construction justifies the additional cost. We believe that systematic software construction will only be adopted by mainstream software development organization when the overhead costs have been greatly reduced. Two approaches to cost mitigation are reuse (amortizing costs over many applications) and automation. For the last four years, NASA Ames has funded the Amphion project, whose objective is to automate software reuse through techniques from systematic software construction. In particular, deductive program synthesis (i.e., program extraction from proofs) is used to derive a composition of software components (e.g., subroutines) that correctly implements a specification. The construction of reuse libraries of software components is the standard software engineering solution for improving software development productivity and quality.

Lowry, Michael R.↗
NASA NTRS1972

Users manual for Streamtube Curvature Analysis: Analytical method for predicting the pressure distribution about a nacelle at transonic speeds, volume 1

The computer program, Streamtube Curvature Analysis, is described for the engineering user and for the programmer. The user oriented documentation includes a description of the mathematical governing equations, their use in the solution, and the method of solution. The general logical flow of the program is outlined and detailed instructions for program usage and operation are explained. General procedures for program use and the program capabilities and limitations are described. From the standpoint of the grammar, the overlay structure of the program is described. The various storage tables are defined and their uses explained. The input and output are discussed in detail. The program listing includes numerous comments so that the logical flow within the program is easily followed. A test case showing input data and output format is included as well as an error printout description.

Keith, J. S.↗
DOE OSTI2022

A universal method to compare parts from STEP files

Abstract Model Based Definition (MBD) captures the complete specification of a part in digital form and leverages (at least) the universal “Standard for the Exchange of Product” (STEP) file format. MBD has revolutionized manufacturing due to time and cost savings associated with containing all engineering data within a single digital source. This work presents a novel method to transform digital definitions in any given STEP file into a tensor-like structure that is unique for each part and can be used to regenerate the original STEP file completely. Resulting STEP tensors are amenable to part comparison based on various part specifications in a general and straightforward manner. Here, part similarity is evaluated among sets of parts according to specific geometry, material composition, and design intent. Importantly, specification similarity can be quantified using only the tensors’ structure. As such, this approach is not limited to families of geometric shapes, part types, or fabrication methods; nor does it require any prior knowledge about the parts being compared.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Validation and bias quantification of criticality safety codes for SRNS operations

The validation and bias quantification methods employed by the Criticality Safety Engineering group at Savannah River Nuclear Solutions (SRNS) are discussed. Both SCALE and MCNP are validated against benchmark experiments, where these benchmark models have been divided into subgroups representative of the types of fissile material operations at SRNS. A quantification of the bias is made for each subgroup for each code and the three methods employed for quantification are discussed. The bias information is employed directly in the establishment of the Upper Safety Limit (USL) value of the neutron multiplication factor when using standardized documented SRNS methods. A simple example of validation and bias quantification is given and the procedure for setting the USL using the derived bias information is also presented. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Adaptive Power Flow Approximations With Second-Order Sensitivity Insights

The power flow equations are fundamental to power system planning, analysis, and control. However, the inherent non-linearity and non-convexity of these equations present formidable obstacles in problem-solving processes. To mitigate these challenges, recent research has proposed adaptive power flow linearizations that aim to achieve accuracy over wide operating ranges. The accuracy of these approximations inherently depends on the curvature of the power flow equations within these ranges, which necessitates considering second-order sensitivities. In this paper, we leverage second-order sensitivities to both analyze and improve power flow approximations. We evaluate the curvature across broad operational ranges and subsequently utilize this information to inform the computation of various sample-based power flow approximation techniques. Additionally, we leverage second-order sensitivities to guide the development of rational approximations that yield linear constraints in optimization problems. In conclusion, this approach is extended to enhance accuracy beyond the limitations of linear functions across varied operational scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2004

Design and Experimental Verification of Deployable/Inflatable Ultra-Lightweight Structures

Because launch cost of a space structural system is often proportional to the launch volume and mass and there is no significant gravity in space, NASA's space exploration programs and various science missions have stimulated extensive use of ultra-lightweight deployable/inflatable structures. These structures are named here as Highly Flexible Structures (HFSs) because they are designed to undergo large displacements, rotations, and/or buckling without plastic deformation under normal operation conditions. Except recent applications to space structural systems, HFSs have been used in many mechanical systems, civil structures, aerospace vehicles, home appliances, and medical devices to satisfy space limitations, provide special mechanisms, and/or reduce structural weight. The extensive use of HFSs in today's structural engineering reveals the need of a design and analysis software and a database system with design guidelines for practicing engineers to perform computer-aided design and rapid prototyping of HFSs. Also to prepare engineering students for future structural engineering requires a new and easy-to- understand method of presenting the complex mathematics of the modeling and analysis of HFSs. However, because of the high flexibility of HFSs, many unique challenging problems in the modeling, design and analysis of HFSs need to be studied. The current state of research on HFSs needs advances in the following areas: (1) modeling of large rotations using appropriate strain measures, (2) modeling of cross-section warpings of structures, (3) how to account for both large rotations and cross- section warpings in 2D (two-dimensional) and 1D structural theories, (4) modeling of thickness thinning of membranes due to inflation pressure, pretension, and temperature change, (5) prediction of inflated shapes and wrinkles of inflatable structures, (6) development of efficient numerical methods for nonlinear static and dynamic analyses, and (7) filling the gap between geometrically exact elastic analysis and elastoplastic analysis. The objectives of this research project were: (1) to study the modeling, design, and analysis of deployable/inflatable ultra-lightweight structures, (2) to perform numerical and experimental studies on the static and dynamic characteristics and deployability of HFSs, (3) to derive guidelines for designing HFSs, (4) to develop a MATLAB toolbox for the design, analysis, and dynamic animation of HFSs, and (5) to perform experiments and establish an adequate database of post-buckling characteristics of HFSs.

Pai, P. Frank↗
DOE OSTI2026

Automated and Accelerated Continuum Model Development for Electrochemical Systems (Abbreviated Report)

Despite the availability of computational resources and advancements in numerical computing capabilities, the multiscale models core to understanding, predicting the behaviors of, and designing energy and environmental systems involving porous media are still 1.) developed through by-hand derivations and 2.) limited by many methodological assumptions employed during model derivation. As a result, the advancement of effective media models for engineering DOE mission-critical systems (e.g., batteries, flow batteries, electrolyzers, geothermal systems, subsurface chemical storage systems, etc.) is slow (i.e., it takes years for models to traverse from stages of “development” to “practical utilization”), hindering our ability to effectively optimize such systems and stay at the cutting-edge of the energy frontier. In this work, we aimed to address these limitations by 1.) automating and accelerating multiscale model derivation via symbolic computing and 2.) develop a novel multiscale modeling methodology for flow and transport through porous media that avoids the typical assumptions hindering previous models. As a result of our efforts, we 1.) developed a hybrid symbolic-numeric code called Fouriera for fully-automating the implementation of multiphysical and phase-field models via the Fourier spectral method for materials science research, and 2.) advanced a multiscale modeling methodology called The Method of Finite Averages that rigorously predicts the behaviors of flow and transport through heterogeneous porous media under the influence of non-local effects and strong advection. Ultimately, these deliverables provide strong foundations from which further efforts can advance multiscale modeling tools and capabilities that do not intrinsically rely on 1.) the speed and mathematical capabilities of humans, nor 2.) the methodological assumptions limiting current models.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Impact of particle arrays on phase separation composition patterns

We examine the symmetry-breaking effect of fixed constellations of particles on the surface-directed spinodal decomposition of binary blends in the presence of particles whose surfaces have a preferential affinity for one of the components. Our phase-field simulations indicate that the phase separation morphology in the presence of particle arrays can be tuned to have a continuous, droplet, lamellar, or hybrid morphology depending on the interparticle spacing, blend composition, and time. In particular, when the interparticle spacing is large compared to the spinodal wavelength, a transient target pattern composed of alternate rings of preferred and non-preferred phases emerges at early times, tending to adopt the symmetry of the particle configuration. We reveal that such target patterns stabilize for certain characteristic length, time, and composition scales characteristic of the pure phase-separating mixture. To illustrate the general range of phenomena exhibited by mixture-particle systems, we simulate the effects of single-particle, multi-particle, and cluster–particle systems having multiple geometrical configurations of the particle characteristic of pattern substrates on phase separation. Our simulations show that tailoring the particle configuration, or substrate pattern configuration, a relative fluid-particle composition should allow the desirable control of the phase separation morphology as in block copolymer materials, but where the scales accessible to this approach of organizing phase-separated fluids usually are significantly larger. Limited experiments confirm the trends observed in our simulations, which should provide some guidance in engineering patterned blend and other mixtures of technological interest.

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

A testpart for interdisciplinary analyses in micro production engineering

In 2011, a round robin test was initiated within the group of CIRP Research Affiliates. The aim was to establish a platform for linking interdisciplinary research in order to share the expertise and experiences of participants all over the world. This paper introduces a testpart which has been designed to allow an analysis of different manufacturing technologies, simulation methods, machinery and metrology as well as process and production planning aspects. Current investigations are presented focusing on the machining and additive processes to produce the geometry, simulation approaches, machine analysis, and a comparison of measuring technologies. Challenges and limitations regarding the manufacturing and evaluation of the testpart features by the applied methods are discussed.

42 ENGINEERING↗
DOE OSTI2018

Special Analysis: Impact of Updated GSA Flow Model on E-Area Low-Level Waste Facility Groundwater Performance

In March 2018, the Savannah River National Laboratory (SRNL) notified Solid Waste Management (SWM) of new information that could potentially impact groundwater (GW) disposal limits and possibly require temporary protective measures (Crowley 2018). GW flow directions in the E-Area Low-Level Waste Facility (LLWF) Slit Trenches (STs), Engineered Trenches (ETs), and the Low Activity Waste Vault (LAWV) have notably changed in the new 2018 General Separations Area (GSA) flow model (Flach 2018a). The 2018 GSA flow model had been updated using 20 years of new hydrologic field data and model calibration methods that employed mathematical optimization software. It was later determined that the primary reason for the change in GW flow directions was due to the influence of low permeability caps that had been placed over the Old Burial Grounds and the Low Level Radioactive Waste Disposal Facility (LLRWDF). Updated flow directions in the model produce a higher degree of plume overlap for disposal units (DU’s) in the southeastern portion of E-Area than had been predicted in the 2008 E-Area LLWF Performance Assessment (PA). SWM and SRNL outlined the scope for this Special Analysis (SA) and proposed a set of interim measures (IM’s) to protect trench operations from exceeding Performance Objectives (PO’s) during the preparation of this SA, both of which were approved by the SWM Performance Assessment Review Committee (PARC) (Mooneyhan 2018).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Flowability of Crumbler Rotary Shear Size-Reduced Granular Biomass: An Experiment-Informed Modeling Study on the Angle of Repose

Biomass has potential as a carbon-neutral alternative to petroleum for chemical and energy products. However, complete replacement of fossil fuel is contingent upon efficient processes to eliminate undesirable characteristics of biomass, e.g., low bulk density, variability, and storage-induced quality problems. Mechanical size reduction via comminution is a processing operation to engineer favorable biomass flowability in handling. Crumbler rotary shear mill has been empirically demonstrated to produce more uniformly shaped particles with higher flowability than hammermilled biomass. This study combines modeling and experimentation to unveil fundamental understandings of the relation between granular particle characteristics and biomass flow behavior, which elucidate underlying mechanisms and guide selection of critical processing parameters. For this purpose, the impact of critical material attributes, including particle size (2–6 mm), particle shape (briquette, chip, clumped-sphere, cube, etc.), and surface roughness, on the angle of repose (AOR) of milled pine chips were investigated using discrete element method (DEM) simulations. Forest Concepts Crumbler rotary shear system is used to produce milled pine particles within the same size range considered in DEM simulations. AOR of different sets of these particles were measured experimentally to benchmark DEM results against experimental data. Specific energy consumption for the comminution of biomass with different particle size and moisture content are measured for technoeconomic analysis. Our results show that the smaller size (2 mm) of pine particle achieves better followability (i.e., smaller AOR) while the energy cost of comminution is significantly higher and bulk density is almost the same as the 6-mm pine particles. For the 2-mm particle size, Crumbles from veneer have better flow properties than Crumbles from chips. Contrarily, no significant difference was observed between the AOR of the two materials for the 6-mm particle size. Furthermore, from DEM simulations, mechanical interlocking between particles was found as a dominant factor in determining AOR of complex-shaped particles such as milled pine, which cannot be accurately captured by using simple particle shapes (e.g., mono-sphere) with a rolling resistance model. Conversely, clumped-sphere model alleviates this limitation without increasing computational cost significantly and can be used for accurate representation of biomass granular particles when simulating free-flow behavior.

09 BIOMASS FUELS↗
DOE OSTI2015

Implementation and Validation of Reference Sensitivity Profile Calculations in TRIPOLI 4 - 14247

In current neutron transport calculations, particularly on continuous-energy Monte-Carlo codes, the major source of uncertainty is coming from nuclear data. In order to improve nuclear data evaluations, integral experiments are used. Feedback is then precisely obtained by re-estimating cross-section levels and correlations according to their sensitivity. k{sub eff}-sensitivity coefficients are traditionally obtained using deterministic codes since hands-on methods were not yet implemented in continuous energy Monte Carlo codes. In this paper, an Iterated-Fission-Probability-based method newly developed in the continuous energy Monte Carlo code TRIPOLI-4 will be used to calculate sensitivity coefficient breakdown on isotopes and reactions and results will be compared to those of ERANOS deterministic code. Objective is twofold: to show the reliability of the new TRIPOLI-4 method and to identify where and why profiles are slightly different. Sensitivity profiles, forward and adjoint fluxes are thus meticulously compared using a fine 1968-group energy structure. An OECD international benchmark (Uncertainty Analysis for Criticality Safety Assessment - UACSA - Phase III) is taken as working example to enforce validation and to compare to independent results (e.g. MCNP-6). This benchmark is constituted of two bare spheres filled with homogeneous mixture of uranium and polyethylene (respectively 2% and 50% U-235 enrichments) and gives an overview of issues encountered in thermal reactor cores. Several energy-group structures (1968G, 172G, 33G) are also tested in ERANOS to identify multi-group and implicit effects. Lessons learnt are that continuous Monte Carlo methods seem to be mature and compare well with deterministic codes. Only damper the U-238 elastic - and on certain conditions the U-238 capture - sensitivity coefficient for which difference can reach 30% between TRIPOLI-4 and ERANOS. Indeed concerning U-238, implicit effects are emphasized and multi-group theory remains limited by hypothesis; even a 1968-group structure is not refined enough to properly calculate U-238 elastic scattering k{sub eff}-sensitivity. This reference 3D continuous-energy sensitivity capability in TRIPOLI4-IFP will be extensively used in the uncertainty automated tool RIB of the French Criticality-Safety package CRISTAL. (authors)

42 ENGINEERING↗
DOE OSTI2023

Learning continuous models for continuous physics

Abstract Dynamical systems that evolve continuously over time are ubiquitous throughout science and engineering. Machine learning (ML) provides data-driven approaches to model and predict the dynamics of such systems. A core issue with this approach is that ML models are typically trained on discrete data, using ML methodologies that are not aware of underlying continuity properties. This results in models that often do not capture any underlying continuous dynamics—either of the system of interest, or indeed of any related system. To address this challenge, we develop a convergence test based on numerical analysis theory. Our test verifies whether a model has learned a function that accurately approximates an underlying continuous dynamics. Models that fail this test fail to capture relevant dynamics, rendering them of limited utility for many scientific prediction tasks; while models that pass this test enable both better interpolation and better extrapolation in multiple ways. Our results illustrate how principled numerical analysis methods can be coupled with existing ML training/testing methodologies to validate models for science and engineering applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

An assessment of the utility of multirate time integration for the modeling of laser powder bed fusion

Finite element simulation of the laser powder bed fusion process is made challenging by the inherently multiscale nature of the process. When using the typical global time stepping techniques, slowly-evolving regions of the domain receive the same numerical treatment as the regions with the highest temperature rates. The current work details the implementation and evaluation of an implicit multirate method which is able to advance different regions of the domain with distinct time steps, depending on their current solution rate. Previous work indicates that this representation of the temporal scales of the problem can lead to significant reductions in wall clock run time, and it is shown herein that multirate time integration, when used with a uniform process-scale mesh, can result in speedups of approximately 19, 42, and 87 times, for domains with edge lengths of approximately 1 mm, 2 mm, and 5 mm, respectively. When used in conjunction with $\ h$-refinement (limited to two levels of refinement), the resulting speedups (taken relative to the uniform mesh) are around 28, 85, and 665 times. The method is demonstrated to converge as indicated by the literature, and its use with an AM-Bench domain is demonstrated.

42 ENGINEERING↗
NASA NTRS1991

A mathematical model of the maximum power density attainable in an alkaline hydrogen/oxygen fuel cell

A mathematical model of a hydrogen/oxygen alkaline fuel cell is presented that can be used to predict the polarization behavior under various power loads. The major limitations to achieving high power densities are indicated and methods to increase the maximum attainable power density are suggested. The alkaline fuel cell model describes the phenomena occurring in the solid, liquid, and gaseous phases of the anode, separator, and cathode regions based on porous electrode theory applied to three phases. Fundamental equations of chemical engineering that describe conservation of mass and charge, species transport, and kinetic phenomena are used to develop the model by treating all phases as a homogeneous continuum.

Kimble, Michael C.↗
DOE OSTI2017

Aeroelastic oscillations of a cantilever with structural nonlinearities: theory and numerical simulation.

Our study details the derivation of the nonlinear equations of motion for the axial, biaxial bending and torsional vibrations of an aeroelastic cantilever undergoing rigid body (pitch) rotation at the base. The primary attenstion is focussed on the geometric nonlinearities of the system, whereby the aeroelastic load is modeled by the theory of linear quasisteady aerodynamics. This modelling effort is intended to mimic the wind-tunnel experimental setup at the Royal Military College of Canada. While the derivation closely follows the work of Hodges and Dowell [1] for rotor blades, this aeroelastic system contains new inertial terms which stem from the fundamentally different kinematics than those exhibited by helicopter or wind turbine blades. Using the Hamilton’s principle, a set of coupled nonlinear partial differential equations (PDEs) and an ordinary differential equation (ODE) are derived which describes the coupled axial-bending-bending-torsion-pitch motion of the aeroelastic cantilever with the pitch rotation. The finite dimensional approximation of the coupled system of PDEs are obtained using the Galerkin projection, leading to a coupled system of ODEs. Subsequently, these nonlinear ODEs are solved numerically using the built-in MATLAB implicit ODE solver and the associated numerical results are compared with those obtained using Houbolt’s method. It is demonstrated that the system undergoes coalescence flutter, leading to a limit cycle oscillation (LCO) due to coupling between the rigid body pitching mode and teh flexible mode arising from the flapwise bending motion.

42 ENGINEERING↗