Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical calculations: variational”

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 289 records · Page 16

Time Dependence of Collision Probabilities During Satellite Conjunctions

The NASA Conjunction Assessment Risk Analysis (CARA) team has recently implemented updated software to calculate the probability of collision (P (sub c)) for Earth-orbiting satellites. The algorithm can employ complex dynamical models for orbital motion, and account for the effects of non-linear trajectories as well as both position and velocity uncertainties. This “3D P (sub c)” method entails computing a 3-dimensional numerical integral for each estimated probability. Our analysis indicates that the 3D method provides several new insights over the traditional “2D P (sub c)” method, even when approximating the orbital motion using the relatively simple Keplerian two-body dynamical model. First, the formulation provides the means to estimate variations in the time derivative of the collision probability, or the probability rate, R (sub c). For close-proximity satellites, such as those orbiting in formations or clusters, R (sub c) variations can show multiple peaks that repeat or blend with one another, providing insight into the ongoing temporal distribution of risk. For single, isolated conjunctions, R (sub c) analysis provides the means to identify and bound the times of peak collision risk. Additionally, analysis of multiple actual archived conjunctions demonstrates that the commonly used “2D P (sub c)” approximation can occasionally provide inaccurate estimates. These include cases in which the 2D method yields negligibly small probabilities (e.g., P (sub c)) is greater than 10 (sup -10)), but the 3D estimates are sufficiently large to prompt increased monitoring or collision mitigation (e.g., P (sub c) is greater than or equal to 10 (sup -5)). Finally, the archive analysis indicates that a relatively efficient calculation can be used to identify which conjunctions will have negligibly small probabilities. This small-P (sub c) screening test can significantly speed the overall risk analysis computation for large numbers of conjunctions.

Hall, Doyle T.↗

Large Array Channel Capacity in the Presence of Interference

We develop a model for a large array ground receiver system for use in deep-space communications, and analyze the resulting array channel capacity. The model includes effects of array geometry, time-dependent spacecraft orbital trajectory, point and extended interference sources, and elevation-dependent noise and tropospheric channel variations. Channel capacity is expressed as the ratio of determinants of covariance matrices characterizing source, interference, and additive noise, and then reduced to a simpler quadratic form more amenable to analysis and numerical computation. This formulation facilitates inclusion of array and channel characteristics into the model, as well as comparison of optimal, suboptimal, and equivalent single antenna configurations on achievable throughput. Realistic examples of ground array channel capacity calculations are presented, demonstrating the impact of array geometry, planetary interference sources, and array combining algorithm design upon the achievable data throughput.

Vilnrotter, V.↗

Digital Twin Model for Advanced Manufacture of a Rotating Detonation Engine Injector

A digital twin material model (DTMM) of an additive manufacturing (AM) process was created to advance the state of the art in rotating detonation engine (RDE) injector design. Current RDE injectors are designed with large pressure drops, enabling a stable and repeatable combustion process. However, this comes at the cost of system efficiency. For the technology to transition to commercial fossil-based power generation, it is important to develop injectors with reduced flow losses. Low-loss injectors are difficult to design and manufacture with conventional manufacturing techniques. AM enables new design options, but the AM manufacturing process must be thoroughly understood to result in a robust design. A DTMM provides the necessary insight by defining the cause-effect relationships between process parameters, microstructure features, and properties. Therefore, a DTMM to support the design and manufacturing process was developed and applied to the design of a new additively manufactured low-loss injector. The injector combustion behavior was characterized through hot-fire tests, and mechanical performance was compared to the DTMM predictions. The two project goals were the successful development of the DTMM and the demonstration of an improved RDE injector design. The RDE injector design and DTMM developments occurred on parallel but dependent paths. The injector was designed to reduce pressure drop by increasing the cross-sectional flow area ratio between the injector air passages and the combustor annulus. This resulted in less structural material, raising the concern that thin members would be susceptible to high-cycle fatigue (HCF) under the periodic loading inherent to an RDE. It was most important for the DTMM to predict behavior in these features; therefore, the injector design concept guided the material thicknesses used in fatigue tests. The DTMM development started by manufacturing a series of coupons over the range of possible AM process variations. A design-of-experiment approach was used to select which process variable combinations gave the most efficient coverage relevant to the injector design space. The microstructure in each of these coupons was characterized, and then computational methods were used to create a numerical model of the correlation between process variables and microstructure. Next, a set of HCF samples were tested to calibrate existing models that map microstructure to HCF performance. Together, these two links formed the DTMM that calculates HCF behavior from AM process variables. Two injector prototypes were additively manufactured. The first injector design strategy aggressively pursued low-loss performance by substantially increasing the oxidizer flow area. The combination of manufacturing lead times and the fatigue testing schedule meant that the DTMM was not available when building this first prototype. Therefore, its process parameters were chosen based on a manual review of the available coupon data. This prototype was built successfully and evaluated in 58 combustion tests. Sustained detonation was achieved with remarkably reduced pressure loss, and some tests even displayed pressure loss characteristics similar to conventional gas turbine combustors. This achieved the project goal of improving RDE injector design. The second injector was manufactured according to the optimized parameters predicted by the DTMM. The flow area modifications of this injector were less aggressive than the first injector since demonstrating low pressure loss was not an objective of the second hot-fire test series. Rather, the test objective was to cause high cycle fatigue failure in the part due to periodic loading from the rotating detonation wave. The observed number of cycles to failure was to be compared to the number predicted by the DTMM and thereby assess the utility of the DTMM in component design. However, the required level of vibration was not obtained during combustion. Therefore, high cycle fatigue was not experienced in the hot-fire tests of the second injector. Fatigue data was obtained by further testing the second injector in a conventional HCF test apparatus. The injector demonstrated HCF strength above the DTMM prediction. In fact, it did not fail and testing was only discontinued due to reaching the end of the period of performance. This points to some success in the project’s primary goal of successfully developing and applying the DTMM to a component design. Implementing the DTMM recommendations for optimal processing parameters led to a part with acceptable properties. The DTMM was also shown to be an efficient correlator of data and to provide insight into the relationship between process settings, microstructure, and property performance. However, the failure of the DTMM prediction to match the experimental result of the injector fatigue test also points to the need to include significantly more data in the model development. In this project, coupons made with identical processing parameters exhibited drastically different properties from each other and from the injector part, which clearly influences the accuracy of a model that predicts performance based on parameters. Uncertainties in the build process must be quantified to develop more robust models. A denser and broader matrix of coupon process and geometry variations, several repeated builds of every point, more in-situ build process measurements, and direct observation of tensile and HCF sample microstructure (as opposed to separate microstructure specimens) are recommendations to improve future AM modeling efforts.

20 FOSSIL-FUELED POWER PLANTS↗

Impeller shroud to casing leakage flow simulations in the Space Shuttle Main Engine high pressure fuel pump

Quasi-three-dimensional Navier-Stokes calculations were carried out for the Space Shuttle Main Engine high-pressure fuel pump to simulate the impeller shroud to casing leakage flow. This flow geometry was modeled as an axisymmetric cavity flow with a stationary surface representing the casing, and a rotating surface denoting the impeller. A 63 x 81-node mesh provided sufficient resolution in the regions of greatest flow variations and reduced the effects of numerical diffusion. The turbulence field was closed with the high Reynolds number form of the k-epsilon model supplemented with wall functions in the vicinity of the walls. Finally, a parametric study quantified the effects of through mass flow changes on this leakage flow.

Sindir, Munir M.↗

The calculation of rotor/fuselage interaction for two-dimensional bodies

Unsteady rotor wake interactions with the empennage, tail boom, and other aerodynamic surfaces have a significant influence on the aerodynamic performance of the helicopter, ride quality, and vibration. A Computational Fluid Dynamic (CFD) method for computing the aerodynamic interaction between an interacting vortex wake and the viscous flow about arbitrary 2-D bodies was developed to address this helicopter problem. The vorticity and flow field velocities are calculated on a body-fitted computational mesh using an uncoupled iterative solution. The interacting vortex wake is represented by an array of discrete vortices which, in turn, are represented by a finite core model. The evolution of the interacting vortex wake is calculated by Lagrangian techniques. The flow around circular and elliptic cylinders in the absence of an interacting vortex wake was calculated. These results compare very well with other numerical results and with results obtained from experiment and thereby demonstrate the accuracy of the viscous solution. The interaction of a simulated rotor wake with the flow about 2-D bodies, representing cross sections of fuselage components, was calculated to address the vortex interaction problem. The vortex interaction was calculated for the flow about a circular and an elliptic cylinder at 45 and 90 degrees incidence. The results demonstrate the significant variation in lift and drag on the 2-D bodies during the vortex interaction.

Stremel, Paul M.↗

Determination of elastic stresses in gas-turbine disks

A method is presented for the calculation of elastic stresses in symmetrical disks typical of those of a high-temperature gas turbine. The method is essentially a finite-difference solution of the equilibrium and compatibility equations for elastic stresses in a symmetrical disk. Account can be taken of point-to-point variations in disk thickness, in temperature, in elastic modulus, in coefficient of thermal expansion, in material density, and in Poisson's ratio. No numerical integration or trial-and-error procedures are involved and the computations can be performed in rapid and routine fashion by nontechnical computers with little engineering supervision. Checks on problems for which exact mathematical solutions are known indicate that the method yields results of high accuracy. Illustrative examples are presented to show the manner of treating solid disks, disks with central holes, and disks constructed either of a single material or two or more welded materials. The effect of shrink fitting is taken into account by a very simple device.

Manson, S S↗

Redox Interactions between Iron and Carbon in Planetary Mantles: Implications for Degassing and Melting Processes

Carbon stability in planetary mantles has been studied by numerous authors because it is thought to be the source of C-bearing atmospheres and of C-rich lavas observed at the planetary surface. In the Earth, carbonaceous peridotites and eclogites compositions have been experimentally studied at mantle conditions [1] [2] [3]. [4] showed that the fO2 variations observed in martian meteorites can be explained by polybaric graphite-CO-CO2 equilibria in the Martian mantle. Based on thermodynamic calculations [4] and [5] inferred that the stable form of carbon in the source regions of the Martian basalts should be graphite (and/or diamond), and equilibrium with melts would be a source of CO2 for the martian atmosphere. Considering the high content of iron in the Martian mantle (approx.18.0 wt% FeO; [6]), compared to Earth s mantle (8.0 wt% FeO; [7]) Fe/C redox interactions should be studied in more detail.

Martin, A.↗

Marcus–Hush–Chidsey kinetics at electrode–electrolyte interfaces

Electrochemical kinetics at electrode–electrolyte interfaces limit the performance of devices including fuel cells and batteries. While the importance of moving beyond Butler–Volmer kinetics and incorporating the effect of electronic density of states of the electrode has been recognized, a unified framework that incorporates these aspects directly into electrochemical performance models is still lacking. In this work, we explicitly account for the density functional theory-calculated density of states numerically in calculating electrochemical reaction rates for a variety of electrode–electrolyte interfaces. We first show the utility of this for two cases related to Li metal electrodeposition and stripping on a Li surface and a Cu surface (anode-free configuration). The deviation in reaction rates is minor for cases with flat densities of states such as Li, but is significant for Cu due to nondispersive d-bands creating large variation. Finally, we consider a semiconducting case of a solid-electrolyte interphase consisting of LiF and Li2CO3 and note the importance of the Fermi level at the interface pinned by the redox reaction occurring there. We identify the asymmetry in reaction rates as a function of discharge/charge naturally within this approach.

Kurchin, Rachel (ORCID:0000000221474809)↗

Duct liner optimization for turbomachinery noise sources

An acoustical field theory for axisymmetric, multisectioned lined ducts with uniform flow profiles was combined with a numerical minimization algorithm to predict optimal liner configurations having one, two, and three sections. Source models studied include a point source located on the axis of the duct and rotor/outlet-stator viscous wake interaction effects for a typical research compressor operating at an axial flow Mach number of about 0.4. For this latter source, optimal liners for equipartition-of energy, zero-phase, and least-attenuated-mode source variations were also calculated and compared with exact results. It is found that the potential benefits of liner segmentation for the attenuation of turbomachinery noise is greater than would be predicted from point source results. Furthermore, effective liner design requires precise knowledge of the circumferential and radial modal distributions.

Lester, H. C.↗

Global Changes of the Water Cycle Intensity

In this study, we evaluate numerical simulations of the twentieth century climate, focusing on the changes in the intensity of the global water cycle. A new diagnostic of atmospheric water vapor cycling rate is developed and employed, that relies on constituent tracers predicted at the model time step. This diagnostic is compared to a simplified traditional calculation of cycling rate, based on monthly averages of precipitation and total water content. The mean sensitivity of both diagnostics to variations in climate forcing is comparable. However, the new diagnostic produces systematically larger values and more variability than the traditional average approach. Climate simulations were performed using SSTs of the early (1902-1921) and late (1979- 1998) twentieth century along with the appropriate C02 forcing. In general, the increase of global precipitation with the increases in SST that occurred between the early and late twentieth century is small. However, an increase of atmospheric temperature leads to a systematic increase in total precipitable water. As a result, the residence time of water in the atmosphere increased, indicating a reduction of the global cycling rate. This result was explored further using a number of 50-year climate simulations from different models forced with observed SST. The anomalies and trends in the cycling rate and hydrologic variables of different GCMs are remarkably similar. The global annual anomalies of precipitation show a significant upward trend related to the upward trend of surface temperature, during the latter half of the twentieth century. While this implies an increase in the hydrologic cycle intensity, a concomitant increase of total precipitable water again leads to a decrease in the calculated global cycling rate. An analysis of the land/sea differences shows that the simulated precipitation over land has a decreasing trend while the oceanic precipitation has an upward trend consistent with previous studies and the available observations. The decreasing continental trend in precipitation is located primarily over tropical land regions, with some other regions, such as North America experiencing an increasing trend. Precipitation trends are diagnosed further using the water tracers to delineate the precipitation that occurs because of continental evaporation, as opposed to oceanic evaporation. These diagnostics show that over global land areas, the recycling of continental moisture is decreasing in time. However, the recycling changes are not spatially uniform so that some regions, most notably over the United States, experience continental recycling of water that increases in time.

Bosilovich, Michael G.↗

Modeling of Melt Growth During Carbothermal Processing of Lunar Regolith

The carbothermal processing of lunar regolith has been proposed as a means to produce carbon monoxide and ultimately oxygen to support human exploration of the moon. In this process, gaseous methane is pyrolyzed as it flows over the hot surface of a molten zone of lunar regolith and is converted to carbon and hydrogen. Carbon gets deposited on the surface of the melt, and mixes and reacts with the metal oxides in it to produce carbon monoxide that bubbles out of the melt. Carbon monoxide is further processed in other reactors downstream to ultimately produce oxygen. The amount of oxygen produced crucially depends on the amount of regolith that is molten. In this paper we develop a model of the heat transfer in carbothermal processing. Regolith in a suitable container is heated by a heat flux at its surface such as by continuously shining a beam of solar energy or a laser on it. The regolith on the surface absorbs the energy and its temperature rises until it attains the melting point. The energy from the heat flux is then used for the latent heat necessary to change phase from solid to liquid, after which the temperature continues to rise. Thus a small melt pool appears under the heated zone shortly after the heat flux is turned on. As time progresses, the pool absorbs more heat and supplies the energy required to melt more of the regolith, and the size of the molten zone increases. Ultimately, a steady-state is achieved when the heat flux absorbed by the melt is balanced by radiative losses from the surface. In this paper, we model the melting and the growth of the melt zone with time in a bed of regolith when a portion of its surface is subjected to a constant heat flux. The heat flux is assumed to impinge on a circular area. Our model is based on an axisymmetric three-dimensional variation of the temperature field in the domain. Heat transfer occurs only by conduction, and effects of convective heat transport are assumed negligible. Radiative heat loss from the surface of the melt and the regolith to the surroundings is permitted. We perform numerical computations to determine the shape and the mass of the melt at steady state and its time evolution. We first neglect the volume change upon melting, and subsequently perform calculations including it. Predictions from our model are compared to test data to determine the effective thermal conductivities of the regolith and the melt that are compatible with the data

Balasubramaniam, R.↗

Assessing the numerical stability of physics models to equilibrium variation through database comparisons on DIII-D

High fidelity kinetic equilibria are crucial for tokamak modeling and analysis. Manual workflows for constructing kinetic equilibria are time consuming and subject to user error, motivating development of automated equilibrium reconstruction tools to provide accurate and consistent reconstructions for downstream physics analysis. These automated tools also provide access to kinetic equilibria at large database scales, which enables the quantification of general uncertainties arising from equilibrium reconstruction techniques. In this paper, we compare a large database of DIII-D kinetic equilibria generated manually by physics experts to equilibria from automated kinetic reconstruction tools, assessing the impact of reconstruction method on equilibrium parameters and resulting magnetohydrodynamic stability calculations. We find agreement among scalar parameters, whereas profile quantities, such as the bootstrap current, show larger disagreements. We analyze ideal kink and classical tearing stability with DCON and STRIDE respectively, finding that the kink stability calculation is generally more robust than the tearing index Δ' calculation. We find that in 90% of cases, both kink stability classifications are unchanged between the manual expert and automated kinetic equilibria.

CAKE↗

Collapse, equilibrium, and fragmentation of rotating, adiabatic clouds

Numerical calculations of the collapse of adiabatic clouds from uniform density and rotation initial conditions show that when restricted to axisymmetry, the clouds form either near-equilibrium spheroids or rings. Rings form in the collapse of low thermal energy clouds and have a ratio of rotational kinetic energy to the absolute value of gravitational potential energy greater than approximately 0.43. When the axisymmetric constraint is removed and an initial m = 2 density variation is introduced, clouds either collapse to form near-equilibrium ellipsoids or else fragment into binary systems through a bar phase. Ellipsoids form in the collapse of high thermal energy clouds and have a rotational kinetic energy/absolute value of gravitational potential energy ratio less than approximately 0.27. The results are consistent with the critical values of the rotational kinetic energy/absolute value of gravitational potential energy ratio for instabilities in Maclaurin spheroids, and suggest that protostellar clouds may undergo a dynamic fragmentation in the nonisothermal collapse regime.

Boss, A. P.↗

A numerical study of the effect of geometry variation, turbulence models, and dissipation on the flow past control surfaces

The steady 3D viscous flow past the ONERA M6 wing and a slender delta wing-body with trailing edge control surfaces has been computed. A cell-centered finite-volume Navier-Stokes patched zonal method has been used for the numerical simulation. Both diagonalized and LUSGS schemes have been implemented. Besides the standard nonplanar zonal interfacing techniques, a new virtual zone capability has been employed. For code validation, the transonic flow past the ONERA M5 wing is calculated for angles-of-attack of 3.06 deg and 5.06 deg and compared with the available experiments. The wing-body computational results are compared with experimental data for both trailing-edge flaps deflected. The experimental flow conditions are M subinfinity = 0.4, a turbulent Reynolds number of 5.41 million based on a mean aerodynamic chord of 25.959 inches, adiabatic wall, and angles-of-attack varying from 0 deg to 23.85 deg. The computational results are presented for the 23.85 deg angle-of-attack case. The effects of the base flow due to a model sting, the varying second and fourth order numerical dissipation, and the turbulence model are all considered.

Chaussee, Denny S.↗

Solar Energy Deposition Rates in the Mesosphere Derived from Airglow Measurements: Implications for the Ozone Model Deficit Problem

We derive rates of energy deposition in the mesosphere due to the absorption of solar ultraviolet radiation by ozone. The rates are derived directly from measurements of the 1.27-microns oxygen dayglow emission, independent of knowledge of the ozone abundance, the ozone absorption cross sections, and the ultraviolet solar irradiance in the ozone Hartley band. Fifty-six months of airglow data taken between 1982 and 1986 by the near-infrared spectrometer on the Solar-Mesosphere Explorer satellite are analyzed. The energy deposition rates exhibit altitude-dependent annual and semi-annual variations. We also find a positive correlation between temperatures and energy deposition rates near 90 km at low latitudes. This correlation is largely due to the semiannual oscillation in temperature and ozone and is consistent with model calculations. There is also a suggestion of possible tidal enhancement of this correlation based on recent theoretical and observational analyses. The airglow-derived rates of energy deposition are then compared with those computed by multidimensional numerical models. The observed and modeled deposition rates typically agree to within 20%. This agreement in energy deposition rates implies the same agreement exists between measured and modeled ozone volume mixing ratios in the mesosphere. Only in the upper mesosphere at midlatitudes during winter do we derive energy deposition rates (and hence ozone mixing ratios) consistently and significantly larger than the model calculations. This result is contrary to previous studies that have shown a large model deficit in the ozone abundance throughout the mesosphere. The climatology of solar energy deposition and heating presented in this paper is available to the community at the Middle Atmosphere Energy Budget Project web site at http://heat-budget.gats-inc.com.

Mlynczak, Martin G.↗

A finite-step method for estimating the spanwise lift distribution of wings in symmetric, yawed, and rotary flight at low speeds

The finite-step method was programmed for computing the span loading and stability derivatives of trapezoidal shaped wings in symmetric, yawed, and rotary flight. Calculations were made for a series of different wing planforms and the results compared with several available methods for estimating these derivatives in the linear angle of attack range. The agreement shown was generally good except in a few cases. An attempt was made to estimate the nonlinear variation of lift with angle of attack in the high alpha range by introducing the measured airfoil section data into the finite-step method. The numerical procedure was found to be stable only at low angles of attack.

Krenkel, A. R.↗

Intensification and refraction of acoustical signals in partially choked converging ducts

A computer code based on the wave-envelope technique is used to perform detailed numerical calculations for the intensification and refraction of sound in converging hard walled and lined circular ducts carrying high mean Mach number flows. The results show that converging ducts produce substantial refractions toward the duct center for waves propagating against near choked flows. As expected, the magnitude of the refraction decreases as the real part of the admittance increases. The pressure wave pattern is that of interference among the different modes, and hence the variation of the magnitude of pressure refraction with frequency is not monotonic.

Nayfeh, A. H.↗

Configuration and impurity quantification of AmLi sources using radiography and gamma spectroscopy

This paper presents the non-destructive assay of two AmLi sources used for detector characterization and active interrogation at PNNL, identified as MRC-101 and MRC-103. First, a detailed description of the internal configuration of 2724-BT encapsulated MRC sources is established. X-ray radiographs of the AmLi sources are reported for the first time to verify the dimensions and orientation of the encapsulation. Notably, a density variation was observed as a contrast change across the height of the source cavities. Contrary to commonly made assumptions in literature, it is indirectly observed that the AmLi source/target material only occupies a portion of the source cavity. The intensity of the prompt gamma-rays resulting from the inelastic scatter of alpha particles on 7Li was used to calculate the mass ratios of LiH:AmO2. The gamma-ray spectra revealed the presence of numerous impurities in both source and target material (243Am, 237Np, 154Eu, 23Na, and 9Be). Neutron-gamma production ratios were applied to the intensity of the Doppler-broadened peaks from alpha-induced reactions with 7Li, 9Be, and 23Na to calculate neutron rates for both sources. A neutron/alpha production yield (derived in this work) was then related to the atomic concentration of these isotopes to find the approximate masses within each source. Ultimately, the AmLi source density was estimated with prompt gamma analysis based on partial source volumes as 0.45±0.08 g/cm^3.

DTRA, AmLi, neutron↗