Engineering Papers⌕ Search

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.

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

At least 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

The Aqua-Planet Experiment (APE): CONTROL SST Simulation

Climate simulations by 16 atmospheric general circulation models (AGCMs) are compared on an aqua-planet, a water-covered Earth with prescribed sea surface temperature varying only in latitude. The idealised configuration is designed to expose differences in the circulation simulated by different models. Basic features of the aqua-planet climate are characterised by comparison with Earth. The models display a wide range of behaviour. The balanced component of the tropospheric mean flow, and mid-latitude eddy covariances subject to budget constraints, vary relatively little among the models. In contrast, differences in damping in the dynamical core strongly influence transient eddy amplitudes. Historical uncertainty in modelled lower stratospheric temperatures persists in APE.Aspects of the circulation generated more directly by interactions between the resolved fluid dynamics and parameterized moist processes vary greatly. The tropical Hadley circulation forms either a single or double inter-tropical convergence zone (ITCZ) at the equator, with large variations in mean precipitation. The equatorial wave spectrum shows a wide range of precipitation intensity and propagation characteristics. Kelvin mode-like eastward propagation with remarkably constant phase speed dominates in most models. Westward propagation, less dispersive than the equatorial Rossby modes, dominates in a few models or occurs within an eastward propagating envelope in others. The mean structure of the ITCZ is related to precipitation variability, consistent with previous studies.The aqua-planet global energy balance is unknown but the models produce a surprisingly large range of top of atmosphere global net flux, dominated by differences in shortwave reflection by clouds. A number of newly developed models, not optimised for Earth climate, contribute to this. Possible reasons for differences in the optimised models are discussed.The aqua-planet configuration is intended as one component of an experimental hierarchy used to evaluate AGCMs. This comparison does suggest that the range of model behaviour could be better understood and reduced in conjunction with Earth climate simulations. Controlled experimentation is required to explore individual model behavior and investigate convergence of the aqua-planet climate with increasing resolution.

AGCM↗
NASA NTRS2009

Mean Element Propagations Using Numerical Averaging

The long-term evolution characteristics (and stability) of an orbit are best characterized using a mean element propagation of the perturbed two body variational equations of motion. The averaging process eliminates short period terms leaving only secular and long period effects. In this study, a non-traditional approach is taken that averages the variational equations using adaptive numerical techniques and then numerically integrating the resulting EOMs. Doing this avoids the Fourier series expansions and truncations required by the traditional analytic methods. The resultant numerical techniques can be easily adapted to propagations at most solar system bodies.

orbit propagations↗
NASA NTRS1959

Longitudinal Laminar Flow Between Cylinders Arranged in Regular Array

The increasing complexity of heat transfer and process situations which involve fluid flow has demanded the frequent use of flow passages of unusual geometrical configuration. The present investigation is concerned with one such novel configuration, namely the longitudinal flow between solid cylindrical rods which are arranged in regular array. A schematic diagram of the situation under study. The rods may be located either in triangular or square array. The flow will be taken to be laminar and fully developed. The aim of this analysis is to determine the pressure drop, shear stress, and velocity-distribution characteristics of the system. The starting point of this study is the basic law of momentum conservation. The resulting differential equation has been solved in an approximate, but almost exact, manner by the use of truncated trigonometric series. Results are obtained over a wide range of porosity values for both the triangular and square arrays. Heat transfer has not been considered. The configuration under investigation has potential application in compact heat exchangers for nuclear reactors and other situations. Further the results should also be of interest in the theory of flow through unconsolidated porous beds (ia, 9a). The only related analytical work known to the authors is that of Emersleben (S), who considered only the square array. His rather involved solution, based on complex zeta functions, appears to be valid only at high porosities. Experiments covering a porosity range of 0.093 to 0.984 have been made by Sullivan (4) using parallel-oriented fibers, most of the tests being for fibers in random array. These previous investigations will be compared with the present theory in a later section.

Sparrow, E. M.↗
NASA NTRSDate not supplied

Comparison of Two Load Prediction Methods for Strain-Gage Balances

Data from a high-capacity semi-span balance is used to perform a detailed comparison of the load prediction accuracy of two strain-gage balance load prediction methods. Both methods independently obtain their load prediction equations from multivariate least squares fits of balance calibration data. The first method is called Non-Iterative Method. This approach directly uses regression models of the individual load components of a balance for the load prediction. The second method is called Iterative Method. This alternate approach uses a load iteration equation for the load prediction that is constructed from the regression models of the gage outputs of the balance. Basic characteristics of the two methods are reviewed. Afterwards, both methods are applied to calibration and check load data of the chosen balance. Finally, selected analysis results are compared. These comparisons confirmed that the load prediction accuracy of the two methods is the same for all practical purposes.

wind tunnel test↗
NASA NTRSDate not supplied

Comparison of Two Load Prediction Methods for Strain-Gage Balances

Data from a five-component semi-span balance is used to perform a systematic comparison of the load prediction accuracy of two load prediction methods. Both methods independently obtain the load prediction equations from multivariate least squares fits of balance calibration data. The first method is called the Non-Iterative Method. This approach directly uses regression models of the individual load components of a balance for the load prediction. The second method is called the Iterative Method. This alternate approach uses a load iteration equation for the load prediction that is constructed from the regression coefficients of the gage outputs of the balance. Basic characteristics of the two methods are reviewed. Afterwards, both methods are applied to calibration, check load, and wind tunnel test data of a five-component semi-span balance. Selected analysis results are compared. These comparisons confirm that the accuracy of the two methods is the same for all practical purposes.

strain-gage balance↗
NASA NTRS1987

Stability of a rigid rotor supported on oil-film journal bearings under dynamic load

Most published work relating to dynamically loaded journal bearings are directed to determining the minimum film thickness from the predicted journal trajectories. These do not give any information about the subsynchronous whirl stability of journal bearing systems since they do not consider the equations of motion. It is, however, necessary to know whether the bearing system operation is stable or not under such an operating condition. The stability characteristics of the system are analyzed. A linearized perturbation theory about the equilibrium point can predict the threshold of stability; however it does not indicate postwhirl orbit detail. The linearized method may indicate that a bearing is unstable for a given operating condition whereas the nonlinear analysis may indicate that it forms a stable limit cycle. For this reason, a nonlinear transient analysis of a rigid rotor supported on oil journal bearings under: (1) a unidirectional constant load, (2) a unidirectional periodic load, and (3) variable rotating load are performed. The hydrodynamic forces are calculated after solving the time-dependent Reynolds equation by a finite difference method with a successive overrelaxation scheme. Using these forces, equations of motion are solved by the fourth-order Runge-Kutta method to predict the transient behavior of the rotor. With the aid of a high-speed digital computer and graphics, the journal trajectories are obtained for several different operating conditions.

Majumdar, B. C.↗
NASA NTRS1986

Geostrophic adjustment in a shallow-water numerical model as it relates to thermospheric dynamics

The theory of geostrophic adjustment and its application to the dynamics of the high latitude thermosphere have been discussed in previous papers based on a linearized treatment of the fluid dynamical equations. However, a linearized treatment is only valid for small Rossby numbers given by Ro = V/fL, where V is the wind speed, f is the local value of the Coriolis parameter, and L is a characteristic horizontal scale for the flow. For typical values in the auroral zone, the approximation is not reasonable for wind speeds greater than 25 m/s or so. A shallow-water (one layer) model was developed that includes the spherical geometry and full nonlinear dynamics in the momentum equations in order to isolate the effects of the nonlinearities on the adjustment process. A belt of accelerated winds between 60 deg and 70 deg latitude was used as the initial condition. The adjustment process was found to proceed as expected from the linear formulation, but that an asymmetry between the response for an eastward and westward flow results from the nonlineawr curvature (centrifugal) terms. In general, the amplitude of an eastward flowing wind will be less after adjustment than a westward wind. For instance, if the initial wind velocity is 300 m/s, the linearized theory predicts a final wind speed of 240 m/s, regardless of the flow direction. However, the nonlinear curvature terms modify the response and produce a final wind speed of only 200 m/s for an initial eastward wind and a final wind speed of almost 300 m/s for an initial westward flow direction. Also, less gravity wave energy is produced by the adjustment of the westward flow than by the adjustment of the eastward flow. The implications are that the response of the thermosphere should be significantly different on the dawn and dusk sides of the auroral oval. Larger flow velocities would be expected on the dusk side since the plasma will accelerate the flow in a westward direction in that sector.

Larsen, M. F.↗
NASA NTRS1985

Dynamic response characteristics of a circulation control rotor model pneumatic system

Numerical and experimental simulation of unsteady airflow through the control valve and slotted air duct of a circulation control rotor is described. The numerical analysis involves the solution of the quasi-one-dimensional compressible fluid-dynamic equations in the blade air duct together with the coupled isentropic flow equations for flow into the blade through the valve and out of the blade through the Coanda slot. Numerical solutions are compared with basic experimental results obtained for a mockup of a circulation control rotor and its pneumatic valving system. The pneumodynamic phenomena that were observed are discussed with particular emphasis on the characteristic system time lags associated with the response of the flow variables to transient and periodic control valve inputs.

Watkins, C. B.↗
DOE OSTI2025

The absence of ray-effects in the discrete ordinate solution to the transport equation in spherical coordinates in multi-dimensions

The streaming operator of the transport equation is derived for spherical coordinates by starting from Newton’s second law for a free particle expressed in spherical coordinates. We shall show that the partial derivatives with respect to the velocity variables of the particle, which are absent in the Cartesian coordinate formulation of the transport equation, arise in the spherical coordinate formulation of the transport equation in response to the centrifugal force which prevents a free particle from ‘falling into the origin’ of the coordinate system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1972

Fixed base simulator study of an externally blown flap STOL transport airplane during approach and landing

A fixed-base simulator study was conducted to determine the flight characteristics of a representative STOL transport having a high wing and equipped with an external-flow jet flap in combination with four high-bypass-ratio fan-jet engines during the approach and landing. Real-time digital simulation techniques were used. The computer was programed with equations of motion for six degrees of freedom and the aerodynamic inputs were based on measured wind-tunnel data. A visual display of a STOL airport was provided for simulation of the flare and touchdown characteristics. The primary piloting task was an instrument approach to a breakout at a 200-ft ceiling with a visual landing.

Grantham, W. D.↗
DOE OSTI2025

Dynamics of heavy quarks in strongly coupled $\mathcal{N}$ = 4 SYM plasma

We calculate the probability distribution P(k) for a heavy quark with velocity v propagating through strongly coupled N = 4 SYM plasma in the ’t Hooft limit (N c → ∞, λ = g 2 N c → ∞) at a temperature T to acquire a momentum k due to interactions with the plasma. This distribution encodes the well-known drag coefficient η D and the transverse and longitudinal momentum diffusion coefficients κ T and κ L . The jet quenching parameter $\hat{q}$ can be extracted from P(k) for v = 1. Going beyond these known Gaussian characteristics of P(k), our calculation determines all of the higher order and mixed moments to leading order in 1/$\sqrt{λ}$ for the first time. These non-Gaussian features of P(k) include qualitatively novel correlations between longitudinal energy loss and transverse momentum broadening at nonzero v. We show that all higher moments scale characteristically with an effective temperature of the boosted plasma in the heavy quark rest frame, and we demonstrate that these non-Gaussian characteristics can be sizable in magnitude and even dominant in physically relevant situations. We use these results to derive a Kolmogorov equation for the evolution of the probability distribution for the total momentum of a heavy quark that propagates through strongly coupled plasma. This evolution equation accounts for all higher order correlations between transverse momentum broadening and longitudinal energy loss, which we have calculated from first principles. It reduces to a Fokker-Planck equation when truncated to only include the effects of η D , κ T and κ L . Remarkably, while heavy quarks do not reach kinetic equilibrium with the plasma if evolved with this Fokker-Planck equation, by showing that the Boltzmann distribution is a static solution of the all-order Kolmogorov equation that we have derived we demonstrate that heavy quarks do reach kinetic equilibrium if evolved with this equation. Our results thus provide a dynamically complete framework for understanding the thermalization of a heavy quark that may be initially far from equilibrium in the strongly coupled N = 4 SYM plasma — as well as new insight into heavy quark transport and equilibration in quark-gluon plasma.

Holography and Hydrodynamics↗
DOE OSTI2024

Degradation of performance in ICF implosions due to Rayleigh–Taylor instabilities: A Hamiltonian perspective

The Rayleigh–Taylor instability (RTI) is an ubiquitous phenomenon that occurs in inertial-confinement-fusion (ICF) implosions and is recognized as an important limiting factor of ICF performance. To analytically understand the RTI dynamics and its impact on ICF capsule implosions, we develop a first-principle variational theory that describes an imploding spherical shell undergoing RTI. The model is based on a thin-shell approximation and includes the dynamical coupling between the imploding spherical shell and an adiabatically compressed fluid within its interior. Using a quasilinear analysis, we study the degradation trends of key ICF performance metrics (e.g., stagnation pressure, residual kinetic energy, and areal density) as functions of initial RTI parameters (e.g., the initial amplitude and Legendre mode), as well as the 1D implosion characteristics (e.g., the convergence ratio). We compare analytical results from the theory against nonlinear results obtained by numerically integrating the governing equations of this reduced model. Our findings emphasize the need to incorporate polar flows in the calculation of residual kinetic energy and demonstrate that higher convergence ratios in ICF implosions lead to significantly greater degradation of key performance metrics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Gravothermal evolution of dark matter halos with differential elastic scattering

Here, we study gravothermal evolution of dark matter halos in the presence of differential self-scattering that has strong velocity and angular dependencies. We design controlled N-body simulations to model Rutherford and Møller scatterings in the halo, and follow its evolution in both core-expansion and -collapse phases. The simulations show the commonly-used transfer cross section underestimates the effects of dark matter self-interactions, but the viscosity cross section provides an accurate approximation for modeling angular-dependent dark matter scattering. We investigate thermodynamic properties of the halo, and find that the three moments of the Boltzmann equation under the fluid approximation are satisfied. We further propose a constant effective cross section, which integrates over the halo's characteristic velocity dispersion with weighting kernels motivated by kinetic theory of heat conduction. The effective cross section provides a good approximation to differential self-scattering for most of the halo evolution. It indicates that we can map astrophysical constraints on a constant self-interacting cross section to an SIDM model with velocity- and angular-dependent scatterings.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Salt marsh soil organic carbon is regulated by drivers of microbial activity

Abstract Soil organic carbon is the foundation for soil health and a livable climate. Organic carbon is concentrated in coastal wetland soils, but dynamics that govern carbon persistence in coastal ecosystems remain incompletely understood. Whether microbial activity results in a gain or loss of carbon depends on environmental conditions that regulate microbial community attributes. We sought to identify which drivers of microbial activity have the greatest impact on organic carbon content in salt marsh soils. To address this question, we used the PRISMA (Preferred Reporting Items for Systematic reviews and Meta-analyses) statement to compile data on soil and ecosystem characteristics from 50 studies of over 60 salt marshes located around the world. We conducted a meta-analysis with structural equation modeling, including mediation and moderation analyses, to identify environmental drivers of salt marsh soil organic carbon content. High salinity, pH, nitrogen, and phosphorus were associated with increased microbial biomass carbon and soil organic carbon. Correlations between microbial biomass and organic carbon were strengthened by soil salinity and nitrogen, and weakened by soil water content. These results suggest that environmental conditions that control microbial growth and activity have potential to preserve or degrade organic carbon in salt marsh soils.

Erb, Hailey (ORCID:0000000337295634)↗
DOE OSTI2020

Impact of dynamical chiral symmetry breaking and dynamical diquark correlations on proton generalized parton distributions

We calculate the leading-twist, helicity-independent generalized parton distributions (GPDs) of the proton, at finite skewness, in the Nambu–Jona-Lasinio (NJL) model of quantum chomodynamics (QCD). The NJL model reproduces low-energy characteristics of QCD, including dynamical chiral symmetry breaking (DCSB). The proton bound-state amplitude is solved for using the Faddeev equation in a quark-diquark approximation, including both dynamical scalar and axial vector diquarks. GPDs are calculated using a dressed non-local correlator, consistent with DCSB, which is obtained by solving a Bethe-Salpeter equation. The model and approximations used observe Lorentz covariance, and as a consequence the GPDs obey polynomiality sum rules. Extractions of electromagnetic and gravitational form factors are performed. We find a D-term of -1.08 when the non-local correlator is properly dressed, and 0.85 when the bare correlator is used instead, suggesting that within this framework proton stability requires the constituent quarks to be dressed consistently with DCSB. In conclusion, we also find that the anomalous gravitomagnetic moment vanishes, as required by Poincare symmetry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Estimating the Effects of Soil Remediation on Children’s Blood Lead near a Former Lead Smelter in Omaha, Nebraska, USA

Lead exposures from legacy sources threaten children’s health. Soil in Omaha, Nebraska, was contaminated by emissions from a lead smelter and refinery. The U.S. Environmental Protection Agency excavated and replaced contaminated soil at the Omaha Lead Superfund Site between 1999 and 2016. The goal of this study was to assess the association of soil lead level (SLL) and soil remediation status with blood lead levels (BLLs) in children living near or on the site. We linked information on SLL at residential properties with children’s BLLs and assigned remediation status to children’s BLL measurements based on whether their measurements occurred during residence at remediated or unremediated properties. We examined the association of SLL and remediation status with elevated BLL (EBLL). We distinguished the roles of temporal trend and the intervention with time-by-intervention-status interaction contrasts. All analyses estimated odds ratios (ORs) with a generalized estimating equations approach to ensure robustness under the complex correlations among BLL measurements. All analyses controlled for relevant covariates including children’s characteristics. EBLL (> 5 μg/dL ) was associated with both residential SLL [e.g., OR=2.00; 95% confidence interval (CI): 1.83, 2.19; > 400–800 vs. ≤ 200 ppm] and neighborhood SLL [e.g., OR=1.85 (95% CI: 1.62, 2.11; > 400–800 vs. ≤ 200 ppm)] before remediation but only with neighborhood SLL after remediation. The odds of EBLL were higher before remediation [OR 1.52 (95% CI: 1.34, 1.72)]. Similarly, EBLL was positively associated with preremediation status in our interaction analysis [interaction OR=1.18 (95%CI: 1.02, 1.37)]. Residential and neighborhood SLLs were important predictors of EBLLs in children residing near or on this Superfund site. Neighborhood SLL remained a strong predictor following remediation. Our data analyses showed the benefit of soil remediation. Results from the interaction analyses should be interpreted cautiously due to imperfect correspondence of remediation times between remediation and comparison groups. https://doi.org/10.1289/EHP8657

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Context-Aware Learning for Inverse Design in Photovoltaics

This document describes progress in the ARPA-E DIFFERENTIATE project titled “Context-Aware Learning for Inverse Design in Photovoltaics” during the period of May 2019 to May 2022. This project is being performed at Iowa State University, New York University, Stanford University, and National Renewable Energy Laboratory. The project aims to develop a new machine learning (ML) framework to significantly accelerate the design of organic microstructures for improved organic photovoltaic performance. In this project, we had developed an inverse design framework using Deep Learning called InvNets for generating microstructures with desired physics-driven properties. As a preliminary product, in Milestone 3, we demonstrated how InvNets show 20% improvement in the performance of the microstructures and over 100X speedup in the performance compared to traditional processes for physics-driven inverse design. Later, in Milestone 6, we demonstrated that InvNets work for more complex physics properties, specifically, generating microstructures for organic photovoltaic cells with desired current-voltage characteristics. Further, in Milestone 4, we explored the idea of using physics-aware surrogates for obtaining solutions of partial differential equations(PDE) called as DiffNets(now called as NeuFENets to avoid ambiguity of names). The connection between both frameworks is that DiffNet surrogates form the physics-aware surrogate in the InvNet framework. Finally in Milestone 8, we extend our framework for other physics domains. Specifically, we explore building geometry-aware NeuFENets by developing physics surrogates that exploit ideas from traditional immersed boundary finite element methods. With these updates, we are able to achieve all the Milestones.

36 MATERIALS SCIENCE↗
NASA NTRS1975

Polarization tracking of a partially coherent signal using a double loop

We consider the problem of simultaneously estimating the phase and polarization angles of a linearly polarized electromagnetic wave. The solution of this problem has applications to space communications and ionospheric probing. The nonlinear stochastic characteristics of a cross-coupled double-loop implementation for the tracker are evaluated using the multidimensional Fokker-Planck equation. The exact stationary probability density is found for a special case of equal signal-to-noise ratio (SNR) in the two loops. The stationary probability density is obtained for a general class of double loops when the two loop bandwidths differ greatly.

Ohlson, J. E.↗