SEARCH · Search NASA

Results for “EQUATIONS - DRAG, PROFILE - BICKNELL”

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 397 records · Page 22

All sources NASA NTRS DOE OSTI

NASA NTRS2019

SI IV Resonance Line Emission during Solar Flares: Non-LTE, Non-equilibrium, Radiation Transfer Simulations

The Interface Region Imaging Spectrograph routinely observes the Si IV resonance lines. When analyzing quiescent observations of these lines, it has typically been assumed that they form under optically thin conditions. This is likely valid for the quiescent Sun, but this assumption has also been applied to the more extreme flaring scenario. We used 36 electron-beam-driven radiation hydrodynamic solar flare simulations, computed using the RADYN code, to probe the validity of this assumption. Using these simulated atmospheres, we solved the radiation transfer equations to obtain the non-LTE, nonequilibrium populations, line profiles, and opacities for a model silicon atom, including charge exchange processes. This was achieved using the "minority species" version of RADYN. The inclusion of charge exchange resulted in a substantial fraction of Si IV at cooler temperatures than those predicted by ionization equilibrium. All simulations with an injected energy flux F> 5 x 10(exp 10) erg cm(exp -2) s(exp -1) resulted in optical depth effects on the Si IV emission, with differences in both intensity and line shape compared to the optically thin calculation. Weaker flares (down to F ~ 10(exp 9) erg cm(exp -2) s(exp -1)) also resulted in Si IV emission forming under optically thick conditions, depending on the other beam parameters. When opacity was significant, the atmospheres generally had column masses in excess of 5 x 10-6 g cm-2 over the temperature range 40-100 kK, and the Si IV formation temperatures were between 30 and 60 kK. We urge caution when analyzing Si IV flare observations, or when computing synthetic emission without performing a full radiation transfer calculation.

Kerr, Graham S.↗

DOE OSTI2024

Constraining the Dense Matter Equation of State with New NICER Mass–Radius Measurements and New Chiral Effective Field Theory Inputs

Pulse profile modeling of X-ray data from the Neutron Star Interior Composition Explorer is now enabling precision inference of neutron star mass and radius. Combined with nuclear physics constraints from chiral effective field theory (χEFT), and masses and tidal deformabilities inferred from gravitational-wave detections of binary neutron star mergers, this has led to a steady improvement in our understanding of the dense matter equation of state (EOS). Here, we consider the impact of several new results: the radius measurement for the 1.42 M ⊙ pulsar PSR J0437−4715 presented by Choudhury et al., updates to the masses and radii of PSR J0740+6620 and PSR J0030+0451, and new χEFT results for neutron star matter up to 1.5 times nuclear saturation density. Using two different high-density EOS extensions—a piecewise-polytropic (PP) model and a model based on the speed of sound in a neutron star (CS)—we find the radius of a 1.4 M ⊙ (2.0 M ⊙ ) neutron star to be constrained to the 95% credible ranges ${12.28}_{-0.76}^{+0.50}$ km (${12.33}_{-1.34}^{+0.70}$ km) for the PP model and ${12.01}_{-0.75}^{+0.56}$ km (${11.55}_{-1.09}^{+0.94}$ km) for the CS model. The maximum neutron star mass is predicted to be ${2.15}_{-0.16}^{+0.14}$M ⊙ and ${2.08}_{-0.16}^{+0.28}$M ⊙ for the PP and CS models, respectively. We explore the sensitivity of our results to different orders and different densities up to which χEFT is used, and show how the astrophysical observations provide constraints for the pressure at intermediate densities. Moreover, we investigate the difference R 2.0 − R 1.4 of the radius of 2 M ⊙ and 1.4 M ⊙ neutron stars within our EOS inference.

gravitational wave sources↗

NASA NTRS1991

Adjoint operator approach to shape design for internal incompressible flows

The problem of determining the profile of a channel or duct that provides the maximum static pressure rise is solved. Incompressible, laminar flow governed by the steady state Navier-Stokes equations is assumed. Recent advances in computational resources and algorithms have made it possible to solve the direct problem of determining such a flow through a body of known geometry. It is possible to obtain a set of adjoint equations, the solution to which permits the calculation of the direction and relative magnitude of change in the diffuser profile that leads to a higher pressure rise. The solution to the adjoint problem can be shown to represent an artificially constructed flow. This interpretation provides a means to construct numerical solutions to the adjoint equations that do not compromise the fully viscous nature of the problem. The algorithmic and computational aspects of solving the adjoint equations are addressed. The form of these set of equations is similar but not identical to the Navier-Stokes equations. In particular some issues related to boundary conditions and stability are discussed.

Cabuk, H.↗

NASA NTRS1981

On the OLP prediction of the unstable modes of the flat plate turbulent boundary layer

The Orr et al. (OLP) stability analysis is reformulated to yield profiles from hot-wire velocity covariance data in the boundary layer of a flow over a flat plate. The proper orthogonal decomposition theorem (PODT) is used to maximize the mean square of the inner product of the velocity vector with a deterministic candidate vector function in Hilbert space. The OLP formulation couples the PODT for the eddy viscosity with a global extremum principle for the disturbance kinetic energy. Derivation of the governing equations is described analytically, as are the extraction of the mean velocity profile and the eddy viscosity. Finally, techniques for characterizing the dominant mode of the velocity field are introduced and applied to three-dimensional flow over a flat plate. Good agreement is found between the predictions, using laboratory data, and the measured energy and three-dimensional structures of the dominant eddies.

Wadia, A. R.↗

NASA NTRS1986

Solutions to bi-Maxwellian transport equations for SAR-arc conditions

The first subsonic solutions of the bi-Maxwellian-based 16-moment set of transport equations for stable auroral red (SAR) arc conditions are presented. These are compared with the solutions obtained from the Maxwellian-based 13-moment transport equations for the same boundary conditions. Close agreement between the 16-moment and 13-moment solutions was obtained for the drift velocity, total electron temperature, total proton heat flow, and total electron heat flow profiles. On the other hand, significant discrepancies were found. Thus, the 16-moment density profile falls off more rapidly with increasing altitude than that computed with the 13-moment equations; the total proton temperature is less in the 16-moment case than in the 13-moment case by several thousand degrees at most altitudes; and differences exist in the ratios of the proton and electron temperature anisotropies with the altitude. A simplified set of transport equations was obtained by dropping terms which remain relatively small at all altitudes.

Demars, H. G.↗

NASA NTRS2011

Evolution of the Equatorial Oscillation in Saturn's Stratosphere Between 2005 and 2010 from Cassini/CIRS Limb Data Analysis

We present an analysis of thermal infrared spectra acquired in limb viewing geometry by Cassini/CIRS in February 2010. We retrieve vertical profiles of Saturn's stratospheric temperature from 20 hPa to 10 (exp -2) hPa, at 9 latitudes between 20 deg N and 20 deg S. Using the gradient thermal wind equation, we derive a map of the zonal wind field. Both the temperature and the zonal wind vertical profiles exhibit an oscillation in the equatorial region. These results are compared to the temperature and zonal wind maps obtained from 2005-2006 CIRS limb data, when this oscillation was first reported. In both epochs, strong temperature anomalies at the equator (up to 20K) are consistent with adiabatic heating (cooling) due to a sinking (rising) motion at a speed of 0.1 - 0.2 mm/s. Finally, we show that the altitude of the maximum eastward wind has moved downwards by 1.3 scale heights in 4.2 years, hence with a 'phase' speed of approximately 0.5 mm/s. This rate is consistent with the estimated period of 14.7 years for the equatorial oscillation, and requires a local zonal acceleration of 1.1 x 10(exp -6) m.s(exp -2) at the 2.5 hPa pressure level. This downward propagation of the oscillation is consistent with it being driven by absorption of upwardly propagating waves.

Guerlet, S.↗

DOE OSTI2022

Data for Catalyst and Reactor Design Considerations for Selective Production of Acids by Oxidative Cleavage of Alkenes with H2O2

Oxidative cleavage of alkenes and unsaturated fatty acids with hydrogen peroxide gives an efficient and sustainable process to obtain mono- and di-acids for polymers and lubricants with fewer safety risks and less environmental impact than processes that utilize ozone or other inorganic oxidizers (e.g., permanganate, dichromate, etc.). Guided by insight into the mechanisms for competing reaction pathways (i.e., epoxidation of alkene on W–(η2-O2) complexes vs. H2O2 decomposition) and the apparent kinetics derived from kinetic experiments, here, we postulate that W-based heterogeneous catalysts can provide high performance and stable operations at low H2O2 concentrations. Semi-batch reactors with continuous introduction of H2O2 solutions offer the means to maintain low H2O2 concentrations while providing sufficient quantities of H2O2 to satisfy the reaction stoichiometry. We derived simple kinetic model equations for the epoxidation, ring-opening, oxidative cleavage, and oxidation steps and fit theses equations to batch experimental data to obtain kinetic parameters. This kinetic model describes the concentration profiles of reactant, oxidant, and products well as shown by agreement with experimental data. Further predictions of the optimal H2O2 feed rate for semi-batch operation utilized by the proposed rate expressions and the reactor design equations suggest that low H2O2 feed rate increases selectivity towards oxidative cleavage products and selective use of H2O2 for oxidative cleavage pathway. Comparisons of oxidative cleavage of 4-octene in batch and semi-batch reactors show that semi-batch reactors with optimized molar feed rates of H2O2 increased oxidative cleavage product selectivities (76% to 99%; with an increase in butyric acid selectivity from 1% to 55%) and H2O2 selectivity (3% to 30%). In addition, semi-batch reaction conditions used avoid H2O2-mediated dissolution of W-atoms from the catalyst. Analysis of these findings suggest that solid oxide catalysts will be effective for continuous oxidative cleavage reactions if deployed within fixed-bed reactors that allow for distributed introduction of reactants and therefore low in situ concentrations of H2O2.

Catalysis↗

NASA NTRS2019

A NICER view of PSR J0030+0451: Implications for the Dense Matter Equation of State

Both the mass and radius of the millisecond pulsar PSRJ0030+0451 have been inferred via pulse-profile modeling of X-ray data obtained by NASA’s Neutron Star Interior Composition Explorer (NICER) mission. In this Letter we study the implications of the mass–radius inference reported for this source by Riley et al. for the dense matter equation of state (EoS), in the context of prior information from nuclear physics at low densities. Using a Bayesian framework we infer central densities and EoS properties for two choices of high-density extensions: a piecewise-polytropic model and a model based on assumptions of the speed of sound in dense matter. Around nuclear saturation density these extensions are matched to an EoS uncertainty band obtained from calculations based on chiral effective field theory interactions, which provide a realistic description of atomic nuclei as well as empirical nuclear matter properties within uncertainties. We further constrain EoS expectations with input from the current highest measured pulsar mass; together, these constraints offer a narrow Bayesian prior informed by theory as well as laboratory and astrophysical measurements. The NICER mass–radius likelihood function derived by Riley et al. using pulse-profile modeling is consistent with the highest-density region of this prior. The present relatively large uncertainties on mass and radius for PSR J0030+0451 offer, however, only a weak posterior information gain over the prior. We explore the sensitivity to the inferred geometry of the heated regions that give rise to the pulsed emission, and find a small increase in posterior gain for an alternative (but less preferred) model. Lastly, we investigate the hypothetical scenario of increasing the NICER exposure time for PSRJ0030+0451.

G Raaijmakers↗

NASA NTRS2019

The Impact of Neglecting Ice Phase on Cloud Optical Depth Retrievals from AERONET Cloud Mode Observations

Clouds present many challenges to climate modelling. To develop and verify the parameterisations needed to allow climate models to represent cloud structure and processes, there is a need for high-quality observations of cloud optical depth from locations around the world. Retrievals of cloud optical depth are obtainable from radiances measured by Aerosol Robotic Network (AERONET) radiometers in “cloud mode” using a two-wavelength retrieval method. However, the method is unable to detect cloud phase, and hence assumes that all of the cloud in a profile is liquid. This assumption has the potential to introduce errors into long term statistics of retrieved optical depth for clouds that also contain ice. Using a set of idealized cloud profiles we find that, for optical depths above 20, the fractional error in retrieved optical depth is a linear function of the fraction of the optical depth that is due to the presence of ice cloud (“ice fraction”). Clouds that are entirely ice have positive errors with magnitudes of the order of 55% to 70 %. We derive a simple linear equation that can be used as a correction at AERONET sites where ice fraction can be independently estimated. Using this linear equation, we estimate the magnitude of the error for a set of cloud profiles from five sites of the Atmospheric Radiation Measurement programme. The dataset contains separate retrievals of ice and liquid retrievals; hence ice fraction can be estimated. The magnitude of the error at each location was related to the relative frequencies of occurrence in thick frontal cloud at the mid-latitude sites and of deep convection at the tropical sites – that is, of deep cloud containing both ice and liquid particles. The long-term mean optical depth error at the five locations spans the range 2– 4, which we show to be small enough to allow calculation of top-of-atmosphere flux to within 10% and surface flux to about 15 %.

Shonk, Jonathan K. P.↗

DOE OSTI2025

Advanced Reactor Designs Security Analysis, Risk, and Recommendations: Risks, Consequences, and Possible by-Design Mitigation Approaches Associated with Select Advanced Reactors

Next-generation advanced reactors (ARs) incorporate enhanced safety systems, have smaller source terms, and feature compact modular designs, which should lessen their collective risk profiles. However, to fully evaluate risk, security needs to be a part of the equation. Without taking security into consideration, safety systems and components in the new ARs may be vulnerable to sabotage. These base attributes, coupled with enhanced security features specific to AR design through sound engineering and security-by-design (SeBD), should provide developers and operators with lower inherent security risk profiles. Building security early into the AR design may remove or passively secure potential critical targets from an adversary’s reach , thereby increasing overall safety and security. An integrated approach and diverse design team that includes engineering, operations, and security experts are fundamental to building security into the design without sacrificing fundamental operational efficiencies and principles. The objective of this project was to evaluate the security and safety interfaces for five classes of reactors, identify potential security vulnerabilities of structures, systems, and components (SSC), and underscore the need to consider security alongside safety in the design o f these concepts. The five reactor classes evaluated in this project and presented in this report are molten-salt reactors (MSR), high temperature gas reactors (HTGR), sodium-fast reactors (SFR), advanced light-water reactors (ALWR), and microreactors. These designs were selected because they reflect the concepts that are closest to market deployment and have received significant resource investments from the public and private sector. This project assesses the inherent security risks posed by common classes of ARs, provides a methodology and framework to assess security along with safety, and offers an analysis of potential mitigation strategies that could be incorporated. For each AR technology, the SSCs that relate to radionuclide source safety functions are discussed to understand the SSC contribution to safety and relative importance in the protective strategy for the design. The assumptions that went into evaluating each reactor concept originated from generic publicly available nonproprietary information and should not directly be used to qualify an absolute risk profile nor to rank specific AR designs. Instead, the purpose of the analysis is to understand and compare the generic inherent security risks of different AR technologies.

98 - NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL↗

NASA NTRS1974

A rapid numerical solution to subsonic flow over planar and axisymmetric profiles at an angle of attack of 0 deg

A numerical solution is presented for the incompressible flow over thin planar and axisymmetric profiles at an angle of attack of 0 degrees. The method uses a finite-difference field solution to the governing equation with a Gauss-Seidel successive overrelaxation scheme. However, the use of a simple Cartesian grid system restricts this method to slender profiles. Results are presented for a cambered airfoil, airfoil in wall effect (two-dimensional flowthrough inlet), body of revolution, and flowthrough nacelle. A computer program is presented which can be used for any of the previously mentioned cases with simple input changes. Results for compressible flow are available with the use of the appropriate two-dimensional or axisymmetric compressibility corrections. Computational time for a typical field calculation of 3000 grid points and 200 cycles through the field is less than 1 minute with less than 50,000 octal storage on the Control Data Corporation 6600 computing system.

Mascitti, V. R.↗

NASA NTRS2003

Latent Heating Retrieval from TRMM Observations Using a Simplified Thermodynamic Model

A procedure for the retrieval of hydrometeor latent heating from TRMM active and passive observations is presented. The procedure is based on current methods for estimating multiple-species hydrometeor profiles from TRMM observations. The species include: cloud water, cloud ice, rain, and graupel (or snow). A three-dimensional wind field is prescribed based on the retrieved hydrometeor profiles, and, assuming a steady-state, the sources and sinks in the hydrometeor conservation equations are determined. Then, the momentum and thermodynamic equations, in which the heating and cooling are derived from the hydrometeor sources and sinks, are integrated one step forward in time. The hydrometeor sources and sinks are reevaluated based on the new wind field, and the momentum and thermodynamic equations are integrated one more step. The reevalution-integration process is repeated until a steady state is reached. The procedure is tested using cloud model simulations. Cloud-model derived fields are used to synthesize TRMM observations, from which hydrometeor profiles are derived. The procedure is applied to the retrieved hydrometeor profiles, and the latent heating estimates are compared to the actual latent heating produced by the cloud model. Examples of procedure's applications to real TRMM data are also provided.

Grecu, Mircea↗

DOE OSTI2023

How machine learning can extend electroanalytical measurements beyond analytical interpretation

Electroanalytical measurements are routinely used to estimate material properties exhibiting current and voltage signatures. Analysis of such measurements relies on analytical expressions of material properties to describe the experiments. The need for analytical expressions limits the experiments that can be used to measure properties as well as the properties that can be estimated from a given experiment. Such analytical relations are essentially solutions of the physics-based differential equations (with properties as coefficients) describing the material behavior under certain specific conditions. In recent years, a new machine learning-based approach has been gaining popularity wherein the differential equations are numerically solved to interpret the electroanalytical experiments in terms of corresponding material properties. Since the physics-based differential equations are solved, one can additionally estimate underlying fields, e.g., concentration profile, using such an approach. To exemplify the characteristics of such a machine learning assisted interpretation of electroanalytical measurements, we use data from the Hebb–Wagner test on a magnesium spinel intercalation host. In conclusion, as compared to the traditional analytical expression-based interpretation, the emerging approach decreases experimental efforts to characterize relevant material properties as well as provides field information that was previously inaccessible.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Phase transitions in high-purity zirconium under dynamic compression

In this work, we present results from ramp compression experiments on high-purity Zr that show the α → ω, ω → β, as well as reverse β → ω phase transitions. Simulations with a multi-phase equation of state and phenomenological kinetic model match the experimental wave profiles well. While the dynamic α ! ! transition occurs ~ 9 GPa above the equilibrium phase boundary, the ω → β transition occurs within 0.9 GPa of equilibrium. We estimate that the dynamic compression path intersects the equilibrium ω – β line at P = 29:2 GPa, and T = 490 K. The thermodynamic path in the interior of the sample lies ~ 100 K above the isentrope at the point of the ω → β transition. Approximately half of this dissipative temperature rise is due to plastic work, and half is due to the non-equilibrium α → ω transition. The inferred rate of the α → ω transition is several orders of magnitude higher than that measured in dynamic diamond anvil cell (DDAC) experiments in an overlapping pressure range. We discuss a model for the influence of shear stress on the nucleation rate. We find that the shear stress s ji has the same effect on the nucleation rate as a pressure increase δP = cϵ ij s ji /(ΔV/V); where c is a geometric constant ~ 1 and, ϵ ij are the transformation shear strains. The small fractional volume change ΔV/V ≈ 0:1 at the α → ω transition amplifies the effect of shear stress, and we estimate that for this case δP is in the range of several GPa. Correcting our transition rate to a hydrostatic rate brings it approximately into line with the DDAC results, suggesting that shear stress plays a significant role in the transformation rate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1973

Modeling of turbulent transport in the surface layer

The turbulence equations as written by Donaldson using the method of invariant modeling have been applied to the following limiting cases of the surface or constant flux layer of the planetary boundary layer: (1) Neutrally stable; (2) stable (above influence of surface roughness); (3) nearly neutrally stable; and (4) very unstable (free convection). For the neutrally stable case, the equations are shown to admit as a solution the familiar logarithmic profile. By use of this result, boundary conditions suitable for the surface layer are defined and are simple to apply to rough surfaces.

Smith, G. L.↗

NASA NTRS1975

Turbulent wall jet in a coflowing stream

A theoretical investigation was undertaken to develop a relatively simple model of a two-dimensional, turbulent wall jet in a coflowing stream. The incompressible jet flow was modeled by using an integral method which includes turbulent shear stress, entrainment, and heat transfer. The method solves the conservation equations for the average jet flow properties and uses the velocity profile suggested by Escudier and Nicoll to obtain detailed characteristics of the jet on a flat plate. The analytical results compare favorably with experimental data for a range of injection velocities, which demonstrates the usefulness of the theory for estimating jet growth, velocity decay, and wall skin friction. The theory, which was applied to a Coanda jet on a circular cylinder, provided estimates of suction pressures aft of the jet exit that were in close agreement with experimental values.

Campbell, J. F.↗

NASA NTRS1976

Thermospheric winds and exospheric temperatures from incoherent scatter radar measurements in four seasons

The time-dependent equations of neutral air motion are solved subject to three constraints: two of them are the usual upper and lower boundary conditions and the third is the value of the wind-induced ion drift at any given height. Using incoherent radar data, this procedure leads to a fast, direct numerical integration of the two coupled differential equations describing the horizontal wind components and yields time dependent wind profiles and meridional exospheric neutral temperature gradients. The diurnal behavior of the neutral wind system and of the exospheric temperature is presented for two solstice and two equinox days. The data used were obtained by the St. Santin and the Millstone Hill incoherent scatter radars. The derived geographic distributions of the exospheric temperatures are compared with those predicted by the OGO-6 empirical thermospheric model.

Antoniadis, D. A.↗

NASA NTRS1978

Shear flow beneath oceanic plates - Local nonsimilarity boundary layers for olivine rheology

The principle of local similarity, which has been used to model the two-dimensional boundary layers in the oceanic upper mantle, permits calculation of the temperature, velocity, and stress fields with essentially analytic techniques. Finite difference numerical methods are hard pressed to resolve the detail required by the large variation of viscosity between the lithosphere and the asthenosphere. In this paper the local similarity approximation has been justified by quantitatively evaluating the effect of nonsimilarity due to viscous heating, nonlinear temperature- and pressure-dependent rheology, buoyancy, adiabatic cooling, etc. Nonsimilar effects produce only small modifications of the locally similar boundary layers; important geophysical observables such as surface heat flux and ocean floor topography are given to better than 10 percent by the locally similar solution. A posteriori evaluations of the terms neglected in the boundary layer simplification of the complete equations have been conducted on the locally similar temperature and velocity profiles close to the spreading ridge. The boundary layer models are valid to depths of 100 km at 3 m.y. and 10 km at 0.3 m.y.

Yuen, D. A.↗