SEARCH · Engineering Papers

Results for “characteristic equations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,225 records · Page 68

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Si 1-X Ge X /Si Heterojunction Internal Photoemission Long Wavelength Infrared Detector

Long wavelength Si 1-X Ge X /Si heterojunction internal photoemission (HIP) infrared detectors have been successfully demonstrated utilizing the growth of degenerately boron doped Si 1-X Ge X layers on Si. Recently, Si 0.7 GE 0.3 /SI HIP detectors with either a Si 1-X Ge X single layer or a Si 1-X Ge X /Si multilayer have been demonstrated with cutoff wavelengths out to 23 µm. Near-ideal thermionic emission dark current characteristics were measured and the electrical potential barriers were determined by the Richardson plot. A photoresponse model, similar to the modified Fowler Equation has been developed for the Si 1-X Ge X /Si heterojunction internal photoemission infrared detector at wavelengths corresponding to photon energies less than the Fermi energy. The optical potential barriers, the corresponding cutoff wavelengths, and the emission coefficients, C 1 , for the HIP detectors have been determined from the measured spectral responses using the photoresponse model. Similar optical and thermal potential barriers were obtained.

infrared↗

DOE OSTI2022

Advanced modeling and simulation of research reactors using dynamic mode decomposition

Full text of publication follows. Due to the ever-increasing safety requirements, the current trend of nuclear reactor analysis is shifting towards high-fidelity multi-physics models, which have a very high computational cost and modelling complexity. As the cost of even a single model run makes it impossible to analyse the behaviour and performance of these models on large-scale commercial plants, it has become even more significant to provide suitable benchmarks to validate and test them extensively. In this sense, research reactors offer a promising solution for the initial validation of high-fidelity models, as they are significantly smaller than commercial reactors and their characteristics are well known. In particular, the reactors of the TRIGA family have been used to assess and validate models and methods for Generation-IV designs, as they have some similar features (such as the dominance of natural convection as cooling mechanism and the difficulties in performing sub-channel analysis using standard codes). Still, the computational requirements of high-fidelity models make them unsuitable for real-time analysis, even following their assessment on research reactors. In this sense, Model Order Reduction (MOR) techniques give an additional strategy to reduce the computational cost of high-fidelity models (whilst preserving sufficient accuracy). In particular, this work focuses on Dynamic Mode Decomposition (DMD), a non-intrusive MOR technique that aims at representing models with explicit temporal dynamics by extracting the time-varying characteristics and the governing structures based only on a set of available data, thus without needing any underlying knowledge of the governing equations. In addition, DMD also computes a low-dimensional surrogate of the dynamic matrix of the system, making it suited for stability analysis and real-time evaluations. This work focuses on the application and validation of the DMD method on the Computational Fluid-Dynamics (CFD) model TRIGA Mark II reactor, also discussing in detail the potentiality of this algorithm as an advanced modelling tool for nuclear reactor analysis. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

Computational model for predicting the performance of fluid-powered longiudinal piston modules and motors.

A Matlab-based computational model line code has been developed for predicting the performance of fluid-powered longitudinal piston modules and motors. A system of equations has been developed based upon applying the conservation of mass to a control volume defined by a reciprocating piston. The method results in a differential equation that predicts the pressure gradient in the control volume as a function of rotor rotation angle. The model uses an integration scheme to predict the pressure state in the control volume. The orifice equation is used to predict the flowrate in and out of the control volume. From the resulting solution the torque and speed characteristics of the motor configuration can be predicted, optimized, or tailored for a specific requirement.SAND2019-13983 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Taylor, Michael↗

NASA NTRS2011

MASCOT - MATLAB Stability and Control Toolbox

MASCOT software was created to provide the conceptual aircraft designer accurate predictions of air vehicle stability and control characteristics. The code takes as input mass property data in the form of an inertia tensor, aerodynamic loading data, and propulsion (i.e. thrust) loading data. Using fundamental non-linear equations of motion, MASCOT then calculates vehicle trim and static stability data for any desired flight condition. Common predefined flight conditions are included. The predefined flight conditions include six horizontal and six landing rotation conditions with varying options for engine out, crosswind and sideslip, plus three takeoff rotation conditions. Results are displayed through a unique graphical interface developed to provide stability and control information to the conceptual design engineers using a qualitative scale indicating whether the vehicle has acceptable, marginal, or unacceptable static stability characteristics. This software allows the user to prescribe the vehicle s CG location, mass, and inertia tensor so that any loading configuration between empty weight and maximum take-off weight can be analyzed. The required geometric and aerodynamic data as well as mass and inertia properties may be entered directly, passed through data files, or come from external programs such as Vehicle Sketch Pad (VSP). The current version of MASCOT has been tested with VSP used to compute the required data, which is then passed directly into the program. In VSP, the vehicle geometry is created and manipulated. The aerodynamic coefficients, stability and control derivatives, are calculated using VorLax, which is now available directly within VSP. MASCOT has been written exclusively using the technical computing language MATLAB . This innovation is able to bridge the gap between low-fidelity conceptual design and higher-fidelity stability and control analysis. This new tool enables the conceptual design engineer to include detailed static stability and trim constraints in the conceptual design loop. The unique graphical interface developed for this tool presents the stability data in a format that is understandable by the conceptual designer, yet also provides the detailed quantitative results if desired.

Kenny, Sean↗

DOE OSTI2023

Reactor physics characterization of triply periodic minimal surface-based nuclear fuel lattices

Triply periodic minimal surface (TPMS) lattices are receiving substantial attention in numerous engineering fields due to their impressive topology-driven physical characteristics. TPMS lattices are periodic structures of two distinct intertwined volume domains separated by an area-minimizing surface or wall. TPMS lattices have been observed in nature, such as biological membranes, skeletons, block copolymers, sea urchins, butterfly wings, and equipotential surfaces in crystals. Intriguingly, the topology of TPMS lattices can be easily parametrized via level-set equations and thus are heavily numerically and experimentally studied. Here, a significant research effort is currently applying TPMS lattices for heat exchangers and sinks. This paper extends TPMS lattice applications to nuclear reactor fuel designs, with a focus on identifying relevant TPMS geometric parameters controlling neutronics characteristics, such as reactivity, neutron spectrum, and heat removal properties. We found that fuel surface-area-to-volume ratios for TPMS lattices can be two orders of magnitude larger than current cylindrical fuel rods. Further, the selected TPMS lattice and its implicit equation, the unit cell pitch, wall thickness, and structure porosity are design parameters enabling neutronics optimization for both thermal and fast spectrum configurations, paving the way for exceptionally compact and dense nuclear core concepts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Cauchy problems for Einstein equations in three-dimensional spacetimes

Abstract We analyze existence and properties of solutions of two-dimensional general relativistic initial data sets with a negative cosmological constant, both on spacelike and characteristic surfaces. A new family of such vacuum spacelike data parameterised by poles at the conformal boundary at infinity is constructed. We review the notions of global Hamiltonian charges, emphasizing the difficulties arising in this dimension, both in a spacelike and characteristic setting. One or two, depending upon the topology, lower bounds for energy in terms of angular momentum, linear momentum, and center of mass are established.

Chruściel, Piotr T. (ORCID:0000000183627340)↗

NASA NTRS1970

Dynamic characteristics of a variable-mass flexible missile

The general motion of a variable mass flexible missile with internal flow and aerodynamic forces is considered. The resulting formulation comprises six ordinary differential equations for rigid body motion and three partial differential equations for elastic motion. The simultaneous differential equations are nonlinear and possess time-dependent coefficients. The differential equations are solved by a semi-analytical method leading to a set of purely ordinary differential equations which are then solved numerically. A computer program was developed for the numerical solution and results are presented for a given set of initial conditions.

Meirovitch, L.↗

NASA NTRS1991

Methodology and critical performance characteristics of aerobrake maneuvers for earth return

An effective general methodology for aerobraking in medium and high energy earth returns either with transfer to LEO or return to the surface is presented. Equations are derived or reformulated to indicate performance bounds for a generic aerobrake maneuver whose primary goal is the safe capture of a returning spacecraft. A secondary goal of this maneuver is the ability to maximize the range of entry conditions while minimizing both the strain of thermal and structural loads and the propulsive effort required to achieve a particular destination orbit.

Goodrick, Thomas F.↗

DOE OSTI2023

Framework for phase transitions between the Maxwell and Gibbs constructions

By taking the nucleon-to-quark phase transition within a neutron star as an example, we present a thermodynamically consistent method to calculate the equation of state of ambient matter so that transitions that are intermediate to those of the familiar Maxwell and Gibbs constructions can be described. This method does not address the poorly known surface tension between the two phases microscopically (as, for example, in the calculation of the core pasta phases via the Wigner-Seitz approximation) but instead combines the local and global charge neutrality conditions characteristic of the Maxwell and Gibbs constructions, respectively. Overall charge neutrality is achieved by dividing the leptons to those that obey local charge neutrality (Maxwell) and those that maintain global charge neutrality (Gibbs). The equation of state is obtained by using equilibrium constraints derived from minimizing the total energy density. The results of this minimization are then used to calculate neutron star mass-radius curves, tidal deformabilities, equilibrium and adiabatic sound speeds, and nonradial g-mode oscillation frequencies for several intermediate constructions. Various quantities of interest transform smoothly from their Gibbs structures to those of Maxwell as the local-to-total electron ratio η, introduced to mimic the hadron-to-quark interface tension from 0 (Gibbs) to ∞ (Maxwell), is raised from 0 to 1. As a result, a notable exception is the g-mode frequency for the specific case of η = 1 for which a gap appears between the quark and hadronic branches.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1993

An automatic multigrid method for the solution of sparse linear systems

An automatic version of the multigrid method for the solution of linear systems arising from the discretization of elliptic PDE's is presented. This version is based on the structure of the algebraic system solely, and does not use the original partial differential operator. Numerical experiments show that for the Poisson equation the rate of convergence of our method is equal to that of classical multigrid methods. Moreover, the method is robust in the sense that its high rate of convergence is conserved for other classes of problems: non-symmetric, hyperbolic (even with closed characteristics) and problems on non-uniform grids. No double discretization or special treatment of sub-domains (e.g. boundaries) is needed. When supplemented with a vector extrapolation method, high rates of convergence are achieved also for anisotropic and discontinuous problems and also for indefinite Helmholtz equations. A new double discretization strategy is proposed for finite and spectral element schemes and is found better than known strategies.

Shapira, Yair↗

NASA NTRS1976

Trajectory Traces of Charged Particles in the Magnetosphere

The characteristic enhancements of ring current particles with energies of about 1 to 100 keV, associated with magnetospheric substorms, were observed by Explorer 45 around the plasmapause in the afternoon to midnight region, and showed the characteristic structure called a 'nose' in the proton spectrograms. The motion of these particles in the equatorial magnetosphere, under a recently proposed convection electric field and a dipole magnetic field is described. Approximate equations of a bounce period, a second adiabatic invariant, and a bounce-average azimuthal velocity are given with inaccuracies less than about 0.001 for all pitch angles. The complete set of flow patterns of 90 deg pitch angle particles is also presented by means of stagnation lines through which radial drifts and/or azimuthal drifts change their directions. The particle tracings in the magnetosphere give a basic concept to explain the observed nose characteristics.

Ejiri, M.↗

NASA NTRS1984

Applications of CONMIN to wing design optimization with vortex flow effect

Slender wings on supersonic cruise configurations are expected to be thin and highly swept. As a result, edge-separated vortex flow is inevitable and must be accounted for in aerodynamic analysis and design. The present method is based on the method of suction analogy to calculate the total aerodynamic characteristics. The method requires the solution of the attached flow problem, the latter being solved by a low-order panel method in subsonic and supersonic flow. In essence, the lifting pressure is calculated by using a pressure-doublet distribution satisfying the Prandtl-Glauert equation. From the pressure distribution, the leading-edge suction is calculated. The latter is assumed to be the vortex lift through the method of suction analogy. For a cambered wing, the location of vortex-lift action point is important in predicting the aerodynamic characteristics. It is also seen that the effect of camber shape appears nonlinearly in all aerodynamic expressions. To design the camber shape, the camber slope is represented by a cosine Fourier series at each of several spanwise stations. The Fourier coefficients are the design variables. To design a leading-edge flap in the vortex flow (i.e., a vortex flap), the coordinates of corner points and the deflection angle are the design variables. The process of wing design is to determine the camber shape and twist distribution such that an objective function, typically the drag, is minimized, subject to various constraints.

Lan, C. E.↗

NASA NTRS2013

Combined Experimental and Numerical Simulations of Thermal Barrier Coated Turbine Blades Erosion

A combined experimental and computational study was conducted to investigate the erosion of thermal barrier coated (TBC) blade surfaces by alumina particles ingestion in a single stage turbine. In the experimental investigation, tests of particle surface interactions were performed in specially designed tunnels to determine the erosion rates and particle restitution characteristics under different impact conditions. The experimental results show that the erosion rates increase with increased impingement angle, impact velocity and temperature. In the computational simulations, an Euler-Lagrangian two stage approach is used in obtaining numerical solutions to the three-dimensional compressible Reynolds Averaged Navier-Stokes equations and the particles equations of motion in each blade passage reference frame. User defined functions (UDF) were developed to represent experimentally-based correlations for particle surface interaction models which were employed in the three-dimensional particle trajectory simulations to determine the particle rebound characteristics after each surface impact. The experimentally based erosion UDF model was used to predict the TBC erosion rates on the turbine blade surfaces based on the computed statistical data of the particles impact locations, velocities and angles relative to the blade surface. Computational results are presented for the predicted TBC blade erosion in a single stage commercial APU turbine, for a NASA designed automotive turbine, and for the NASA turbine scaled for modern rotorcraft operating conditions. The erosion patterns in the turbines are discussed for uniform particle ingestion and for particle ingestion concentrated in the inner and outer 5 percent of the stator blade span representing the flow cooling the combustor liner.

Hamed, Awate↗

DOE OSTI2025

Determining the profitability of energy storage over its life cycle using levelized cost of storage

Levelized cost of storage (LCOS) can be a simple, intuitive, and useful metric for determining whether a new energy storage plant would be profitable over its life cycle and to compare the cost of different energy storage technologies. However, researchers and industry decision makers still use conflicting definitions of LCOS. For example, some include charging cost, while others only include round trip efficiency (RTE) losses. Additionally, inputs to the existing formulations are not specific enough to generate repeatable results across studies, which reduces trust in the metric. To push for standardization in economic assessment of batteries and other energy storage devices, the authors review existing definitions of LCOS and identify the desired characteristics for a standard. They then propose a new definition and demonstrate that it fits these characteristics very well relative to other prominent options. Unit analysis is applied to this proposed definition to provide a deeper understanding of the equations and to demonstrate its effectiveness. Finally, the sensitivity of LCOS to different input parameters is investigated to help users understand how to compare analyses from literature to their own. The authors also provide a spreadsheet and a Python script to streamline adoption of the proposed definition.

25 ENERGY STORAGE↗

DOE OSTI2025

Acausality-driven instabilities in relativistic viscous hydrodynamics

We investigate non-linear instabilities stemming from superluminal propagation of information in Israel-Stewart-like models of relativistic viscous fluid dynamics. In relativity, the characteristic speed of propagation of information, $w$, and the speed of the fluid, $v$, allow us to differentiate between regimes of the hydrodynamic equations that are acausal but stable ($w>1$), unstable ($v^{2} w^{2} \geq 1$), and covariantly ill-posed ($w^{2} \leq 0$). As an analytical benchmark, we present a new solution that illustrates these distinct regimes. We compare this analytical solution to the result of a numerical relativistic viscous fluid dynamics solver, and confirm that the analytical result can be recovered numerically in the stable regime, whether causal or acausal. The onset of numerical instabilities is further found to occur in the regime predicted by the analytical solution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Microscopic activated dynamics theory of the shear rheology and stress overshoot in ultradense glass-forming fluids and colloidal suspensions

We formulate a particle and force level, activated dynamics-based statistical mechanical theory for the continuous startup nonlinear shear rheology of ultradense glass-forming hard sphere fluids and colloidal suspensions in the context of the elastically collective nonlinear Langevin equation approach and a generalized Maxwell model constitutive equation. Activated structural relaxation is described as a coupled local-nonlocal event involving caging and longer range collective elasticity which controls the characteristic stress relaxation time. Theoretical predictions for the deformation-induced enhancement of mobility, the onset of relaxation acceleration at remarkably low values of stress, strain, or shear rate, apparent power law thinning of the steady-state structural relaxation time and viscosity, a nonvanishing activation barrier in the shear thinning regime, an apparent Herschel–Buckley form of the shear rate dependence of the steady-state shear stress, exponential growth of different measures of a yield or flow stress with packing fraction, and reduced fragility and dynamic heterogeneity under deformation were previously shown to be in good agreement with experiments. The central new question we address here is the defining feature of the transient response—the stress overshoot. In contrast to the steady-state flow regime, understanding the transient response requires an explicit treatment of the coupled nonequilibrium evolution of structure, elastic modulus, and stress relaxation time. We formulate a new quantitative model for this aspect in a physically motivated and computationally tractable manner. Theoretical predictions for the stress overshoot are shown to be in good agreement with experimental observations in the metastable ultradense regime of hard sphere colloidal suspensions as a function of shear rate and packing fraction, and accounting for deformation-assisted activated motion appears to be crucial for both the transient and steady-state responses.

Mechanics↗

NASA NTRS1969

An Analog Computer Study of the Low-Frequency Flow Dynamics of Two Nuclear-Rocket Cold-Flow Engine Systems

An analog computer study is made of the low-frequency flow dynamics observed with liquid hydrogen flowing in two different nuclear-rocket cold-flow (unfueled) engine systems. The prediction of the low-frequency response depends largely upon the generation of the proper mass-flow storage in the system. This is achieved by relating, in the model lumping, the average values of fluid enthalpy and a certain thermodynamic fluid property. A two-phase state equation is derived for thermodynamic equilibrium between the gas and liquid phases. Parametric studies are made which demonstrate the effect of the pump characteristics on the flow damping.

Fox, Harry W., Jr.↗

NASA NTRS1973

Investigation of gaseous propellant combustion and associated injector/chamber design guidelines

Injector design criteria are provided for gaseous hydrogen-gaseous oxygen propellants. Design equations and procedures are presented which will allow an injector-chamber designer to a priori estimate of the performance, compatibility and stability characteristics of prototype injectors. The effects of chamber length, element geometry, thrust per element, mixture ratio, impingement angle, and element spacing were evaluated for four element concepts and their derivatives. The data from this series of tests were reduced to a single valued mixing function that describes the mixing potential of the various elements. Performance, heat transfer and stability data were generated for various mixture ratios, propellant temperatures, chamber pressures, contraction ratios, and chamber lengths. Applications of the models resulted in the design of procedures, whereby the performance and chamber heat flux can be calculated directly, and the injector stability estimated in conjunction with existing models.

Calhoon, D. F.↗