Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupling coefficient”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A Kinetic Model for Predicting Trace Gas Uptake and Reaction

A model is developed to describe trace gas uptake and reaction with applications to aerosols and microdroplets. Gas uptake by the liquid is formulated as a coupled equilibria that links gas, surface, and bulk regions of the droplet or solution. Previously, this framework was used in explicit stochastic reaction-diffusion simulations to predict the reactive uptake kinetics of ozone with droplets containing aqueous aconitic acid, maleic acid, and sodium nitrite. With the use of prior data and simulation results, a new equation for the uptake coefficient is derived, which accounts for both surface and bulk reactions. Lambert W functions are used to obtain closed form solutions to the integrated rate laws for the multiphase kinetics; similar to previous expressions that describe Michaelis-Menten enzyme kinetics. Together these equations couple interface and bulk processes over a wide range of conditions and do not require many of the limiting assumptions needed to apply resistor model formulations to explain trace gas uptake and reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1999

Development of Analysis Methods for Designing with Composites

The project involved the development of new analysis methods to achieve efficient design of composite structures. We developed a complex variational formulation to analyze the in-plane and bending coupling response of an unsymmetrically laminated plate with an elliptical cutout subjected to arbitrary edge loading as shown in Figure 1. This formulation utilizes four independent complex potentials that satisfy the coupled in-plane and bending equilibrium equations, thus eliminating the area integrals from the strain energy expression. The solution to a finite geometry laminate under arbitrary loading is obtained by minimizing the total potential energy function and solving for the unknown coefficients of the complex potentials. The validity of this approach is demonstrated by comparison with finite element analysis predictions for a laminate with an inclined elliptical cutout under bi-axial loading.The geometry and loading of this laminate with a lay-up of [-45/45] are shown in Figure 2. The deformed configuration shown in Figure 3 reflects the presence of bending-stretching coupling. The validity of the present method is established by comparing the out-of-plane deflections along the boundary of the elliptical cutout from the present approach with those of the finite element method. The comparison shown in Figure 4 indicates remarkable agreement. The details of this method are described in a manuscript by Madenci et al. (1998).

Madenci, E.↗
DOE OSTI2023

Development of an accelerator-based neutron source to prototype Mo-99 production, Part II: A liquid LBE loop under a high vacuum

To provide US domestic supply of Mo-99 without using high-enriched uranium (HEU), a subcritical uranium target assembly (UTA) is irradiated by an accelerator-based neutron source to create Mo-99 through fission. Part I of this work discusses the design of a liquid lead–bismuth eutectic (LBE) windowless target for an accelerator-based neutron source development. Part II discusses how to couple this windowless target to an accelerator operating at an ultra-high vacuum and the subcritical UTA cooled by water at room temperature. Due to the windowless design of the target, the liquid LBE flow shares an ultra-high vacuum (<1.3 × 10 -7 Pa or 10 -9 Torr) space with the accelerator. As a result of this shared vacuum space, the LBE system must operate at a high vacuum (10 -3 ~10 -6 Pa or 10 -5 ~10 -8 Torr). A magnetic rotary motion feedthrough unit utilizes magnetic fluid to allow rotation of the pump while maintaining a high vacuum environment. Prior to testing the LBE system under vacuum, a pump curve measurement is performed to estimate flowrate in the system. This measurement also generates data on orifice loss coefficients, which are compared to correlations in literature. The second experiment investigates vacuum level in the LBE system during operation. High vacuum is maintained (10 -3 ~10 -5 Pa or 10 -5 ~10 -7 Torr) during system operation, and a residual gas analyzer (RGA) scan shows that partial pressures of residual gases in the LBE system lower over the duration of LBE system operations. The third experiment investigates the gravity driven liquid LBE flowing out of the target chamber in the return line, which is partially full. If the liquid LBE is not drained quickly enough, flooding in the target chamber could occur. The coefficient n in the Manning equation is found to be around 0.008 s/m 1/3 . The last experiment performed is a demonstration that a vacuum jacket could provide sufficient thermal insulation to allow coupling between 300 °C LBE loop and a water tank at room temperature. In conclusion, the results from these experiments have influenced the development of the neutron source for the future commercial scale Mo-99 production system.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

High-Temperature Thermodynamics of Cerium Silicates, A-Ce 2 Si 2 O 7 , and Ce 4.67 (SiO 4 ) 3 O

Lanthanide disilicates and oxyapatites have potential roles in high temperature applications as thermal (TBC) and environmental barrier coatings (EBC), or possible alteration phases in geological nuclear waste repositories. However, those Ce 3+ -bearing silicates have only been limitedly studied. In this work, we performed detailed structural and thermodynamic investigations on A-Ce 2 Si 2 O 7 (tetragonal, P4 1 ) and Ce 4.67 (SiO 4 ) 3 O (hexagonal, P6 3 /m). The high temperature structural behaviors and coefficients of thermal expansion were determined by in situ high temperature synchrotron X-ray diffraction (HT-XRD) implemented with Rietveld analysis and thermogravimetric analysis coupled with differential scanning calorimetry (TGA-DSC). A-Ce 2 Si 2 O 7 was found to be stable in N 2 and air up to ~1483 K with anisotropic thermal expansion along the a and c axes (α a = 12.3 × 10 -6 K -1 and α c = 12.4 × 10 -6 K -1 ). Ce 4.67 (SiO 4 ) 3 O had a slow partial oxidation between 533 K and 873 K to a new nonstoichiometric phase Ce 3+ 1.67- xCe4 + x Ce 3+ 3 (SiO 4 ) 3 O 1+0.5x , followed by a thermal decomposition to CeO 2 and SiO 2 at ~1000 K in air. By using high temperature oxide melt solution calorimetry at 973 K with lead borate as the solvent, the standard enthalpy of formation was determined for A-Ce 2 Si 2 O 7 (-3825.1 ± 6.0 kJ/mol) and Ce 4.67 (SiO 4 ) 3 O (-7391.3 ± 9.5 kJ/mol). Finally, these thermodynamic parameters were compared with those of CeO 2 , CeSiO 4 , and other silicate oxyapatites for examining their chemical stability in high temperature environments relevant for aeronautical applications, mineral formation, and nuclear fuel cycle.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electron transport in a collisional plasma with multiple ion species in the presence of a magnetic field

The classical work of Braginskii [Zh. Eksp. Teor. Fiz. 33, 459 (1957)] published almost 65 years ago was the first to provide a complete, closed fluid description of a weakly coupled, fully ionized, collisional plasma immersed in a magnetic field. While this fact is not widely known or appreciated, the Braginskii expressions for the electron transport coefficients can under- or overestimate the said coefficients by up to a factor of two for the electron Hall parameter of order unity (with the Hall parameter being proportional to the electron gyro-frequency over the electron–ion collision frequency) and can provide incorrect Hall-parameter scalings for its large values. Starting with the work of Epperlein and Haines [Phys. Fluids 29, 1029 (1986)], several papers attempted to correct the Braginskii electron results with varying degrees of success. Herein, we present our own effort with a hope to finally put this problem to rest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Foundations of magnetohydrodynamics

In this tutorial, a derivation of magnetohydrodynamics (MHD) valid beyond the usual ideal gas approximation is presented. Non-equilibrium thermodynamics is used to obtain conservation equations and linear constitutive relations. When coupled with Maxwell's equations, this provides closed fluid equations in terms of material properties of the plasma, described by the equation of state and transport coefficients. These properties are connected to microscopic dynamics using the Irving–Kirkwood procedure and Green–Kubo relations. Symmetry arguments and the Onsager–Casimir relations allow one to vastly simplify the number of independent coefficients. Importantly, expressions for current density, heat flux, and stress (conventionally Ohm's law, Fourier's law, and Newton's law) take different forms in systems with a non-ideal equation of state. The traditional form of the MHD equations, which is usually obtained from a Chapman–Enskog solution of the Boltzmann equation, corresponds to the ideal gas limit of the general equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Quantum simulation of Lindbladian dynamics via repeated interactions

The Lindblad equation generalizes the Schrödinger equation to quantum systems that undergo dissipative dynamics. The quantum simulation of Lindbladian dynamics is therefore non-unitary, preventing a naive application of state-of-the-art quantum algorithms. Here, we make use of an approximate correspondence between Lindbladian dynamics and evolution based on repeated interaction (RI) CPTP maps to write down a Hamiltonian formulation of the Lindblad dynamics and derive a rigorous error bound on the master equation. Specifically, we show that the number of interactions needed to simulate the Liouvillian within error e scales in most physical scenarios as . This is significant because the error in the Lindbladian approximation to the dynamics is not explicitly bounded in existing quantum algorithms for open system simulations. We then provide quantum algorithms to simulate RI maps using an iterative qubitization approach and Trotter–Suzuki formulas, and specifically show that for iterative qubitization the number of operations needed to simulate the dynamics (for a fixed value of ?) scales as in the limit where a0 (the coefficient 1-norm for the system and bath Hamiltonians) asymptotically dominates over the corresponding factor for the interaction Hamiltonian, which is often the case in weak coupling. This scaling would appear to be optimal if the complexity of ? is not considered, which underscores the importance of considering the error in the Liouvillian that we reveal in this work.

Quantum Computing↗
DOE OSTI2023

Cascade of vestigial orders in two-component superconductors: Nematic, ferromagnetic, s -wave charge- 4 e , and d -wave charge- 4 e states

Electronically ordered states that break multiple symmetries can melt in multiple stages, similarly to liquid crystals. In a partially melted phase, known as vestigial phase, a bilinear made out of combinations of the multiple components of the primary order parameter condenses. Multicomponent superconductors are thus natural candidates for vestigial order since they break both the U⁡(1)-gauge and also time-reversal or lattice symmetries. Here, we use group theory to classify all possible real-valued and complex-valued bilinears of a generic two-component superconductor on a tetragonal or hexagonal lattice. While the more widely investigated real-valued bilinears correspond to vestigial nematic or ferromagnetic order, the little explored complex-valued bilinears correspond to a vestigial charge-4e condensate, which itself can have an underlying s-wave, d x 2 –y 2 -wave, or d xy -wave symmetry. To properly describe the fluctuating regime of the superconducting Ginzburg-Landau action and thus access these competing vestigial phases, we employ both a large-N and a variational method. We show that while vestigial order can be understood as a weak-coupling effect in the large-N approach, it is akin to a moderate-coupling effect in the variational method. Despite these distinctions, both methods yield similar results in wide regions of the parameter space spanned by the quartic Landau coefficients. Specifically, we find that the nematic and ferromagnetic phases are the leading vestigial instabilities, whereas the various types of charge-4⁢e order are attractive albeit subleading vestigial channels. The only exception is for the hexagonal case, in which the nematic and s-wave charge-4e vestigial states are degenerate. Finally, we discuss the limitations of our approach, as well as the implications of our results for the realization of exotic charge-4⁢e states in material candidates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1976

Position error propagation in the simplex strapdown navigation system

The results of an analysis of the effects of deterministic error sources on position error in the simplex strapdown navigation system were documented. Improving the long term accuracy of the system was addressed in two phases: understanding and controlling the error within the system, and defining methods of damping the net system error through the use of an external reference velocity or position. Review of the flight and ground data revealed error containing the Schuler frequency as well as non-repeatable trends. The only unbounded terms are those involving gyro bias and azimuth error coupled with velocity. All forms of Schuler-periodic position error were found to be sufficiently large to require update or damping capability unless the source coefficients can be limited to values less than those used in this analysis for misalignment and gyro and accelerometer bias. The first-order effects of the deterministic error sources were determined with a simple error propagator which provided plots of error time functions in response to various source error values.

Source record↗
NASA NTRS1977

Low-speed power effects on advanced fighter configurations with two-dimensional deflected thrust

Wind-tunnel studies at the Langley Research Center have shown that significant increases in maximum lift coefficient and stability and decreases in drag due to lift are obtained when two-dimensional vectored thrust is used in conjuction with a close-coupled canard. The configuration tested was somewhat above the theoretical minimum drag due to lift because of the sharp leading edge on the biconvex airfoil used on the wing and canard. An effort to design a new configuration which will approach the minimum drag due to lift while maintaining high-lift configurations was completed. The resulting model will incorporate a realistic planform, airfoil section, and twist for a transonic maneuvering configuration.

Paulson, J. W., Jr.↗
NASA NTRS1986

The stability of coronal and prominence magnetic fields

The significance of the localized instabilities is not yet fully understood. The nonlinear coupling of these modes may give rise to an explosive instability, with the modes coupling to longer wavelengths. On the other hand, if the modes saturate early, than the main effect of the instability may be an enhancement of transport coefficients. Nonetheless, the equation provides a simple test for the stability of cylindrical magnetic fields.

Hood, A. W.↗
NASA NTRS1990

An empirical study of FORTRAN programs for parallelizing compilers

Some results are reported from an empirical study of program characteristics that are important in parallelizing compiler writers, especially in the area of data dependence analysis and program transformations. The state of the art in data dependence analysis and some parallel execution techniques are examined. The major findings are included. Many subscripts contain symbolic terms with unknown values. A few methods of determining their values at compile time are evaluated. Array references with coupled subscripts appear quite frequently; these subscripts must be handled simultaneously in a dependence test, rather than being handled separately as in current test algorithms. Nonzero coefficients of loop indexes in most subscripts are found to be simple: they are either 1 or -1. This allows an exact real-valued test to be as accurate as an exact integer-valued test for one-dimensional or two-dimensional arrays. Dependencies with uncertain distance are found to be rather common, and one of the main reasons is the frequent appearance of symbolic terms with unknown values.

Shen, Zhiyu↗
NASA NTRS1993

Thermal stresses in thick laminated composite shells

The paper provides an analytical formulation to investigate the thermomechanical behavior of thick composite shells subjected to a temperature distribution which varies arbitrarily in the radial direction. For illustrative purposes, shells under uniform temperature change are presented. It is found that thermal twist would occur even for symmetric laminated shells. Under uniform temperature rise, results for off-axis graphite/epoxy shells show that extensional-shear coupling can cause tensile radial stress throughout the shell and tensile hoop stress in the inner region. Laminated graphite/epoxy shells can exhibit negative effective thermal expansion coefficients in the longitudinal and transverse directions. Finally, the stacking sequence has a strong influence on the thermal stress distributions.

Yuan, F. G.↗
NASA NTRS2000

A Common Calibration Source Framework for Fully-Polarimetric and Interferometric Radiometers

Two types of microwave radiometry--synthetic thinned array radiometry (STAR) and fully-polarimetric (FP) radiometry--have received increasing attention during the last several years. STAR radiometers offer a technological solution to achieving high spatial resolution imaging from orbit without requiring a filled aperture or a moving antenna, and FP radiometers measure extra polarization state information upon which entirely new or more robust geophysical retrieval algorithms can be based. Radiometer configurations used for both STAR and FP instruments share one fundamental feature that distinguishes them from more 'standard' radiometers, namely, they measure correlations between pairs of microwave signals. The calibration requirements for correlation radiometers are broader than those for standard radiometers. Quantities of interest include total powers, complex correlation coefficients, various offsets, and possible nonlinearities. A candidate for an ideal calibration source would be one that injects test signals with precisely controllable correlation coefficients and absolute powers simultaneously into a pair of receivers, permitting all of these calibration quantities to be measured. The complex nature of correlation radiometer calibration, coupled with certain inherent similarities between STAR and FP instruments, suggests significant leverage in addressing both problems together. Recognizing this, a project was recently begun at NASA Goddard Space Flight Center to develop a compact low-power subsystem for spaceflight STAR or FP receiver calibration. We present a common theoretical framework for the design of signals for a controlled correlation calibration source. A statistical model is described, along with temporal and spectral constraints on such signals. Finally, a method for realizing these signals is demonstrated using a Matlab-based implementation.

Kim, Edward J.↗
NASA NTRS1981

Eye-head coordination during optokinetic stimulation in squirrel monkeys

Head and eye movements in the yaw plane were recorded during and after optokinetic stimulation in squirrel monkeys. 1) Phasic or tonic head deviations to the side of the ocular quick phase occurred in 94% of total recordings (n = 50) during the perstimulus period, and in 75% of recordings (n = 49) during the poststimulus period. Magnitude of mean head deviation was significantly different between perstimulus and poststimulus periods. 2) Head nystagmus associated with eye nystagmus was consistently observed in seven of nine squirrel monkeys during optokinetic stimulation. Squirrel monkeys are thereby less prone to display head nystagmus than either guinea pigs, pigeons or chickens. 3) Slow phase speeds of coupled head and eye nystagmus were subjected to statistical analysis. A highly significant negative correlation was found between slow phase head and eye speeds. The correlation coefficient was - 0.81 at 60 degrees / sec stimulus (n = 119) and -0.72 at 100 degrees / sec stimulus (n = 131). The gaze speed, calculated by summing the head and eye speeds, was 59.1 plus or minus 6.8 / sec at 60 degrees / sec and 92.2 plus or minus 11.4 at 100 degrees / sec stimulus. There was no significant difference between the gaze speed in a free head condition and the eye speed when the head was fixed.

Eye Movements↗
NASA NTRS2009

Data Parallel Line Relaxation (DPLR) Code User Manual: Acadia - Version 4.01.1

Data-Parallel Line Relaxation (DPLR) code is a computational fluid dynamic (CFD) solver that was developed at NASA Ames Research Center to help mission support teams generate high-value predictive solutions for hypersonic flow field problems. The DPLR Code Package is an MPI-based, parallel, full three-dimensional Navier-Stokes CFD solver with generalized models for finite-rate reaction kinetics, thermal and chemical non-equilibrium, accurate high-temperature transport coefficients, and ionized flow physics incorporated into the code. DPLR also includes a large selection of generalized realistic surface boundary conditions and links to enable loose coupling with external thermal protection system (TPS) material response and shock layer radiation codes.

Wright, Michael J.↗
NASA NTRS2022

Satellite Soil Moisture Data Assimilation Impacts on Modeling Weather Variables and Ozone in the Southeastern US – Part 2: Sensitivity to Dry-Deposition Parameterizations

Ozone (O3) dry deposition is a major O3 sink. As a follow-up study of Huang et al. (2021), we quantify the impact of satellite soil moisture (SM) on model representations of this process when different dry-deposition parameterizations are implemented, based on which the implications for interpreting O3 air pollution levels and assessing the O3 impacts on human and ecosystem health are provided. The SM data from NASA's Soil Moisture Active Passive mission are assimilated into the Noah-Multiparameterization (Noah-MP) land surface model within the NASA Land Information System framework, semicoupled with Weather Research and Forecasting model with online Chemistry (WRF-Chem) regional-scale simulations covering the southeastern US. Major changes in the modeling system used include enabling the dynamic vegetation option, adding the irrigation process, and updating the scheme for the surface exchange coefficient. Two dry-deposition schemes are implemented, i.e., the Wesely scheme and a “dynamic” scheme, in the latter of which dry-deposition parameterization is coupled with photosynthesis and vegetation dynamics. It is demonstrated that, when the dynamic scheme is applied, the simulated O3 dry-deposition velocities vd and their stomatal and cuticular portions, as well as the total O3 fluxes Ft, are larger overall; vd and Ft are 2–3 times more sensitive to the SM changes due to the data assimilation (DA). Further, through case studies at two forested sites with different soil types and hydrological regimes, we highlight that, applying the Community Land Model type of SM factor controlling stomatal resistance (i.e., β factor) scheme in replacement of the Noah-type β factor scheme reduced the v(d) sensitivity to SM changes by ∼75 % at one site, while it doubled this sensitivity at the other site. Referring to multiple evaluation datasets, which may be associated with variable extents of uncertainty, the model performance of vegetation, surface fluxes, weather, and surface O3 concentrations shows mixed responses to the DA, some of which display land cover dependency. Finally, using model-derived concentration- and flux-based policy-relevant O3 metrics as well as their matching exposure–response functions, the relative biomass/crop yield losses for several types of vegetation/crops are estimated to be within a wide range of 1 %–17 %. Their sensitivities to the model's dry-deposition scheme and the implementation of SM DA are discussed.

Min Huang↗