Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “strong coupling”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

An evaluation of three two-dimensional computational fluid dynamics codes including low Reynolds numbers and transonic Mach numbers

Full-potential, Euler, and Navier-Stokes computational fluid dynamics (CFD) codes were evaluated for use in analyzing the flow field about airfoils sections operating at Mach numbers from 0.20 to 0.60 and Reynolds numbers from 500,000 to 2,000,000. The potential code (LBAUER) includes weakly coupled integral boundary layer equations for laminar and turbulent flow with simple transition and separation models. The Navier-Stokes code (ARC2D) uses the thin-layer formulation of the Reynolds-averaged equations with an algebraic turbulence model. The Euler code (ISES) includes strongly coupled integral boundary layer equations and advanced transition and separation calculations with the capability to model laminar separation bubbles and limited zones of turbulent separation. The best experiment/CFD correlation was obtained with the Euler code because its boundary layer equations model the physics of the flow better than the other two codes. An unusual reversal of boundary layer separation with increasing angle of attack, following initial shock formation on the upper surface of the airfoil, was found in the experiment data. This phenomenon was not predicted by the CFD codes evaluated.

Hicks, Raymond M.↗
NASA NTRS2004

Neutron Star Structure In The Presence of Scalar Fields

Motivated by the possible presence of scalar fields on astrophysical scales, suggested by the recent measurement of the deceleration parameter by supernovae surveys, we present models of neutron star structure under the assumption that a scalar field makes a significant contribution to the stress energy momentum tensor, in addition to that made by the normal matter. To that end we solve the coupled Einstein - scalar field - hydrostatic balance equations to compute the effect of the presence of the scalar field on the neutron star structure. We find that the presence of the scalar field does change the structure of the neutron star, especially in cases of strong coupling between the scalar field and the matter density. We present the neutron star radius as a function of the matter-scalar field coupling constant for different values of the neutron star central density. The presence of the scalar field does affect both the maximum neutron star mass and its radius, the latter increasing with the value of the above coupling constant. Our results may be testable with the recent timing observations of accreting neutron stars.

Crawford, James P.↗
NASA NTRS2004

Neutron Star Structure in the Presence of Scalar Fields

Motivated by the possible presence of scalar fields on astrophysical scales, suggested by the apparent acceleration of the universe implied by the supernovae surveys, we present models of neutron star structure including the contribution of a (massless) scalar field to the stress energy momentum tensor, in addition to that made by the normal matter. To that end we solve the coupled Einstein -- scalar field -- hydrostatic balance equations to compute the effect of the presence of the scalar field on the neutron star structure. We find that the presence of the scalar field does change the structure of the neutron star, especially in cases of strong coupling between the scalar field and the matter density. We present the neutron star radius as a function of the matter--scalar field coupling constant for different values of the neutron star central density. The presence of the scalar field affects both the maximum neutron star mass and Its radius, the latter increasing with the value of the above coupling constant. We also compute particle and photon geodesics in the geometry of these neutron stars as well as to the geometry of black holes with different values of the scalar field. Our results may be testable with timing observations of accreting neutron stars.

Kazanas, Demosthenes↗
NASA NTRS1980

Thermal models for solar hard X-ray bursts

Thermal models for hard X-ray bursts consisting of a one-dimensional flux tube whose central electrons are heated to about 400 million K are examined. It is found that the evolution of a thermal X-ray source is a sensitive function of the electron-ion thermal coupling and the state of the plasma into which the source expands. When this coupling is weak, the heated electrons separate into a region of high temperature of about 400 million K and a region of lower temperature of about 100 million K, a process which leads to a power-law X-ray spectrum. In the case of strong coupling there is only one dominant temperature, about 200 million K, and the X-ray spectrum resembles a true thermal spectrum.

Smith, D. F.↗
NASA NTRS2000

An Analysis of the Torsion-Rotation-Vibration Rotational Spectrum of the Lowest In-Plane Bend and First Excited Torsional State of the C(3V) Internal Rotor C2H5CN

C2H5CN (Propionitrile or ethyl cyanide) is a well known interstellar species abundantly observed in hot cores during the onset of star formation. The onset of star formation generally results in elevated temperature, which thermally populates may low lying vibrational states such as the 206/cm in-plane bend and the 212/cm first excited torsional state in C2H5CN. Unfortunately, these two states are strongly coupled through a complex series of torsion-vibration-rotation interactions, which dominate the spectrum. In order to understand the details of these interactions and develop models capable of predicting unmeasured transitions for astronomical observations in C2H5CN and similar molecules, several thousand rotational transitions in the lowest excited in-plane bend and first excited torsional state have been recorded, assigned and analyzed. The analysis reveals very strong a- and b-type Coriolis interactions and a number of other smaller interactions and has a number of important implications for other C3V torsion-rotation-vibration systems. The relative importance and the physical origins of the coupling among the rotational, vibrational and torsional motions will be presented along with a full spectroscopic analysis and supporting astronomical observations.

Pearson, J. C.↗
NASA NTRS1976

Ionospheric disturbances and gravity waves

The response of ionization to a gravity wave moving through the ionosphere is studied. Hydrodynamic equations are used, and local thermodynamic equilibrium is imposed for simplicity. The treatment involves a perturbation analysis, and the background medium is assumed to be time stationary, horizontally stratified, and known. It is shown that ionization may be locally resonant at each level for certain frequencies and directions, for which condition neutral and ionized particles are considered closely or critically coupled. The phase direction for this critical coupling is always downward in the absence of a magnetic field. A magnetic field results in two resonant directions for the same frequency, and these directions are mostly downward. Observed TID's associated with gravity waves may be indicative of such resonances. It is also noted that strong coupling may occur to neutral acoustic waves at high altitudes. Previous investigations restrict their use of momentum equations to the diffusion equation. The analysis also shows that such restrictions result in the neglect of terms arising from momentum transport due to any background ambipolar diffusion velocity and wave motion. These terms are mostly relevant at higher altitudes.

Eun, H.↗
NASA NTRS1987

Studies of pressure-velocity coupling schemes for analysis of incompressible and compressible flows

Two pressure-velocity coupling schemes, both of which solve the fully implicit discretized equations governing the flow of fluids were examined, and the capability of performing large Reynolds number, low Mach number compressible flow calculations were assessed. The semi-implicit iterative SIMPLE algorithm is extended to handle transient compressible flow calculations. This extension takes into account a strong coupling between the pressure and temperature through a correction procedure, based on the equations of state. Results obtained from the extended SIMPLE algorithm are then compared to similar results obtained from the non-iterative PISO algorithm. Both time-dependent and steady state calculations were performed using an axisymmetric 2:1 pipe expansion geometry and laminar flow conditions corresponding to Reynolds number of 1000 and Mach number of 2.0. For calculations simulating a time-dependent compression/expansion process, both schemes exhibit transient features in excellent agreement with each other, and moreover, the PISO method shows a significant computational time reduction of 60 percent over the SIMPLE scheme, regardless of the time step size or grid size employed. The effects of numerical diffusion are shown to be significant in these calculations. For steady state compressible flows, however, the SIMPLE algorithm displays increasing computational efficiency over the PISO method as the time step sizes employed to reach steady state conditions are decreased.

Burgreen, Gregory Wayne↗
NASA NTRS1989

Calculations of inlet distortion induced compressor flow field instability

Calculations of the onset of flow instability are carried out for low-speed multistage axial compressors operating with asymmetric inlet flow. The modeling of the fluid dynamic interaction between the spoiled and unspoiled sectors of the compressor is the most important feature of the calculation procedure. The calculations show that annulus average slope of the compressor pressure rise characteristic equal to zero is a useful approximate stability criterion for situations where the dynamics of the compressor flow field do not couple strongly to the compression system or the structure of the imposed distortion is not similar to that of the eigenmodes of the flow in the compressor annulus. This criterion is employed to investigate the relationship between the present model and the 'parallel compressor' model. Calculations are also presented for cases when compressor and compressor system are closely coupled, and situations in which the compressor is subjected to a rotating distortion. These first-of-a-kind computations, and the accompanying description of the physical mechanisms, show that the stability of the flow in the compressor can be adversely affected if the temporal or spatial structure of the distortion is such that resonant type responses can be evoked either from the compressor or from compressor/compression system interactions.

Chue, R.↗
NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗
NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗
NASA NTRS1994

Coupling of ELF/ULF energy from lightning and MeV particles to the middle atmosphere, inosphere, and global circuit

In an attempt to explain numerous atmospheric electrical phenomena, the elements of the global electrical circuit are reexamined. In addition to being a 'quasi-static 'DC' generator' and source of radiated energy at VLF and higher, the thunderstorm is found to be a pulse generator, with most of the external energy contained in ELF and ULF pulse currents to the ionosphere (and Earth). The pulse energy is found to deposit largely in the middle atmosphere above the thunderstorm. The VLF and above components are well understood, as are the ULF components due to the conductivity gradient. However, a previously poorly understood ELF component on the millsecond timescale, or 'slow tail,' contains a large fraction of the electrical energy. This component couples strongly to the ionosphere and also launches a unipolar transverse electromagnetic (TEM) wavelet in the radial Earth-ionosphere transmission line. The increase in charge with distance associated with such wavelets, and their ensemble sum at a point, may explain some large mesospheric 'DC' fields but there are still difficulties explaining other than rare occurrences, except for antipodal reconvergence. These millisecond duration unipolar wavelets also coupled to the ionosphere and may trigger other lightning at a distance. A schema is elucidated by which the charge of MeV particles deposited in the middle atmosphere persists for much longer than the local relaxation time. This also gives rise to unipolar waves of global extent which may explain lower-latitude field perturbations associated with solar/geomagnetic events.

Hale, Leslie C.↗
NASA NTRS1985

Thermonuclear reactions in cool accreting neutron stars and burst phenomena

The ignition of accreting materials on neutron stars is explored using strongly coupled plasma analytical techniques. The calculations cover the ignition temperature and density at the bottom of the accreted envelope of a neutron star. Emphasis is placed on low-temperature ignitions which take place at high densities. The investigation is extended to the accretion of material from a white dwarf in the form of pure He, C + O, or O + Ne + Mg. It is shown that electrons are strongly degenerate in low-temperature flashes, where the ignition is more dependent on density than on temperature. Precursor flashes 0.4-0.7 the intensity of the main burst will appear before the main bursts. The intensity relationship indicates that the appropriate model for an X-ray burst from a neutron star accreting from a white dwarf is a He shell flash in the presence of a hydrogen-rich atmosphere. The flash will have a maximum energy of 2 x 10 to the 43 ergs and could last as long as 40,000 sec.

Miyaji, S.↗
NASA NTRS1991

Interaction of a finite-length ion beam with a background plasma - Reflected ions at the quasi-parallel bow shock

The coupling of a finite-length, field-aligned, ion beam with a uniform background plasma is investigated using one-dimensional hybrid computer simulations. The finite-length beam is used to study the interaction between the incident solar wind and ions reflected from the earth's quasi-parallel bow shock, where the reflection process may vary with time. The coupling between the reflected ions and the solar wind is relevant to ion heating at the bow shock and possibly to the formation of hot, flow anomalies and re-formation of the shock itself. Consistent with linear theory, the waves which dominate the interaction are the electromagnetic right-hand polarized resonant and nonresonant modes. However, in addition to the instability growth rates, the length of time that the waves are in contact with the beam is also an important factor in determining which wave mode will dominate the interaction. It is found that interaction will result in strong coupling, where a significant fraction of the available free energy is converted into thermal energy in a short time, provided the beam is sufficiently dense or sufficiently long.

Onsager, T. G.↗
NASA NTRS1983

Coupling of radiation into thin film modes by means of localized plasma resonances

The interaction between the surface plasmon mode that propagates at a metal dielectric interface and the localized plasma resonances (LPR) is investigated experimentally in Ag-island films. A stair-stepped sample geometry comprising a glass substrate, a continuous 50-nm Ag film, an LiF spacer film of thickness d = 5-60 nm, and an Ag-island film of mass thickness 3 nm is used in near-normal-reflectivity and plasmon-propagation-constant (k) determinations. The results are presented graphically and discussed. The overall shape of the reflectivity curves is found to be characteristic of Ag films, but with a dip at about 400 nm (corresponding to the absorption resonance of the island film) which is most pronounced with d = 25 nm. It is inferred that the island resonances are strongly coupled to a continuous-film dissipative mechanism at this d value. This inference is supported by the fact that the variation in k, correctd for LiF effects and plotted as a function of d, is greatest at around d = 25 nm. The implications of this finding for broad-band coupling into a thin-film mode, LPR enhancement of waveguide nonlinear effects, and new surface-enhanced-Raman-scattering geometries are indicated.

Holland, W. R.↗
NASA NTRS1981

Parameterization of linear wave chemical transport in planetary atmospheres by eddy diffusion

Two-dimensional transport by linear waves is parameterized by wave velocity and temperature correlations and nonconservative properties of the atmosphere in a residual Eulerian framework. Explicit expressions derived for the eddy tracer flux vector are written in eddy diffusion tensor form, but it is demonstrated that strong coupling of chemistry and temperature through the rate constants, as is the case for ozone, requires additional terms proportional to the mixing ratio and not to its gradient. A variety of limiting cases are discussed. In the residual Eulerian framework the K(yz), K(zy) components are nonzero for a nonconservative atmosphere and/or a reactive tracer. Tracers have common transport coefficients in the limit where their chemistry is negligible in comparison to transience. It is shown that the usual one-dimensional vertical eddy diffusion continuity equation is rigorously valid only when the temperature dependence of the rate constants is negligible and the tracer has strong vertical stratification relative to its horizontal variability.

Strobel, D. F.↗
NASA NTRS1983

Structural dynamics studies of rotating bladed-disk assemblies coupled with flexible shaft motions

In order to analyze the dynamic behavior of the first stage compressor/fan of the 'E3' turbofan engine, a classical structural dynamics approach is employed to couple the motions of a flexible bladed disk to a rotating flexible shaft. The analysis accounts for flexible disk displacements which are transverse to the plane of rotation, and radial as well as tangential, and also accounts for rigid disk translations along, and rotations about, axes normal to the undeformed shaft axes. In the case of a wide range of E3 engine shaft flexibilities and speeds, some of the one-diametral node frequencies are shown to be affected by shaft degrees of freedom whose stiffness values are in general range of design practice. Coriolis forces are also found to significantly affect natural frequencies where strong coupling between certain modes is present.

Loewy, R. G.↗
NASA NTRS1995

Thermal and chemical convection in planetary mantles

Melting of the upper mantle and extraction of melt result in the formation of a less dense depleted mantle. This paper describes series of two-dimensional models that investigate the effects of chemical buoyancy induced by these density variations. A tracer particles method has been set up to follow as closely as possible the chemical state of the mantle and to model the chemical buoyant force at each grid point. Each series of models provides the evolution with time of magma production, crustal thickness, surface heat flux, and thermal and chemical state of the mantle. First, models that do not take into account the displacement of plates at the surface of Earth demonstrate that chemical buoyancy has an important effect on the geometry of convection. Then models include horizontal motion of plates 5000 km wide. Recycling of crust is taken into account. For a sufficiently high plate velocity which depends on the thermal Rayleigh number, the cell's size is strongly coupled with the plate's size. Plate motion forces chemically buoyant material to sink into the mantle. Then the positive chemical buoyancy yields upwelling as depleted mantle reaches the interface between the upper and the lower mantle. This process is very efficient in mixing the depleted and undepleted mantle at the scale of the grid spacing since these zones of upwelling disrupt the large convective flow. At low spreading rates, zones of upwelling develop quickly, melting occurs, and the model predicts intraplate volcanism by melting of subducted crust. At fast spreading rates, depleted mantle also favors the formation of these zones of upwelling, but they are not strong enough to yield partial melting. Their rapid displacement toward the ridge contributes to faster large-scale homogenization.

Dupeyrat, L.↗
NASA NTRS2016

Precipitation and Carbon-Water Coupling Jointly Control the Interannual Variability of Global Land Gross Primary Production

Carbon uptake by terrestrial ecosystems is increasing along with the rising of atmospheric CO2 concentration. Embedded in this trend, recent studies suggested that the interannual variability (IAV) of global carbon fluxes may be dominated by semi-arid ecosystems, but the underlying mechanisms of this high variability in these specific regions are not well known. Here we derive an ensemble of gross primary production (GPP) estimates using the average of three data-driven models and eleven process-based models. These models are weighted by their spatial representativeness of the satellite-based solar-induced chlorophyll fluorescence (SIF). We then use this weighted GPP ensemble to investigate the GPP variability for different aridity regimes. We show that semi-arid regions contribute to 57% of the detrended IAV of global GPP. Moreover, in regions with higher GPP variability, GPP fluctuations are mostly controlled by precipitation and strongly coupled with evapotranspiration (ET). This higher GPP IAV in semi-arid regions is co-limited by supply (precipitation)-induced ET variability and GPP-ET coupling strength. Our results demonstrate the importance of semi-arid regions to the global terrestrial carbon cycle and posit that there will be larger GPP and ET variations in the future with changes in precipitation patterns and dryland expansion.

chlorophyll fluorescence↗