Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Smooth approximation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 289 records · Page 16

Two-dimensional radiative transfer in cloudy atmospheres - The spherical harmonic spatial grid method

A new two-dimensional monochromatic method that computes the transfer of solar or thermal radiation through atmospheres with arbitrary optical properties is described. The model discretizes the radiative transfer equation by expanding the angular part of the radiance field in a spherical harmonic series and representing the spatial part with a discrete grid. The resulting sparse coupled system of equations is solved iteratively with the conjugate gradient method. A Monte Carlo model is used for extensive verification of outgoing flux and radiance values from both smooth and highly variable (multifractal) media. The spherical harmonic expansion naturally allows for different levels of approximation, but tests show that the 2D equivalent of the two-stream approximation is poor at approximating variations in the outgoing flux. The model developed here is shown to be highly efficient so that media with tens of thousands of grid points can be computed in minutes. The large improvement in efficiency will permit quick, accurate radiative transfer calculations of realistic cloud fields and improve our understanding of the effect of inhomogeneity on radiative transfer in cloudy atmospheres.

Evans, K. F.↗

An experimental study of light scattering by large, irregular particles

The intensity and polarization of light scattered by a variety of types of artificial partices large compared to the wavelength were measured as a function of phase angle. Shape, surface roughness, absorption coefficient, and internal scattering coefficient were varied systematically and their effects studied. Scattering by clear, smooth-surfaced spheres is in quantitative agreement with the predictions of the geometrical optics (ray theory) approximation to physical optics (Mie theory). The phase functions of almost all of the particles measured have both forward and backward scattering lobes. A two-parameter, double Henyey-Greenstein function generally provides reasonably good descriptions of the data, while keeping the number of free parameters to the minimum necessary. On a double Henyey- Greenstein parameter plot all of the particles fall into an L-shaped area of restricted size in which the location is characteristic of the particle type. Formalisms based on the equivalent slab model are also given for estimating the scattering efficiency of a large, irregular particle. For most dielectric particles the transmitted, forward scattered light is partially negatively polarized. It is this component that is respopnsible for the well-known maximum in the polarization curves of planetary regoliths at phase angles around 100 deg. For phase angles between about 30 deg and 70 deg the internally scattered light is found to be randomly polarized in the particles studied here, so that the only contribution to the second component of the Stokes vector is by Fresnel reflection from the particle surface. If this empirical result is general, measurement of the second Stokes vector of the light scattered from a regolith at these angles may provide a method of remotely measuring the mean refractive index.

Mcguire, Audrey F.↗

Empirical Equation for Turbulent Forced-Convection Heat Transfer for Prandtl Numbers from 0.001 to 1000

A review is made of some of the experimental data and analyses applicable to convective heat transfer in fully turbulent flow in smooth tubes with liquid metals and viscous Newtonian fluids. An empirical equation is evolved that closely approximates heat-transfer values obtained from selected analyses and experimental data for Prandtl numbers from 0.001 to 1000. The terms included in the equation are Reynolds number, Prandtl number, and an empirical diffusivity ratio between heat and momentum.

vonGlahn, Uwe H.↗

The MSFC Solar Activity Future Estimation (MSAFE) Model

The Natural Environments Branch of the Engineering Directorate at Marshall Space Flight Center (MSFC) provides solar cycle forecasts for NASA space flight programs and the aerospace community. These forecasts provide future statistical estimates of sunspot number, solar radio 10.7 cm flux (F10.7), and the geomagnetic planetary index, Ap, for input to various space environment models. For example, many thermosphere density computer models used in spacecraft operations, orbital lifetime analysis, and the planning of future spacecraft missions require as inputs the F10.7 and Ap. The solar forecast is updated each month by executing MSAFE using historical and the latest month's observed solar indices to provide estimates for the balance of the current solar cycle. The forecasted solar indices represent the 13-month smoothed values consisting of a best estimate value stated as a 50 percentile value along with approximate +/- 2 sigma values stated as 95 and 5 percentile statistical values. This presentation will give an overview of the MSAFE model and the forecast for the current solar cycle.

Suggs, Ron↗

The MSFC Solar Activity Future Estimation (MSAFE) Model

The Natural Environments Branch of the Engineering Directorate at Marshall Space Flight Center (MSFC) provides solar cycle forecasts for NASA space flight programs and the aerospace community. These forecasts provide future statistical estimates of sunspot number, solar radio 10.7 cm flux (F10.7), and the geomagnetic planetary index, Ap, for input to various space environment models. For example, many thermosphere density computer models used in spacecraft operations, orbital lifetime analysis, and the planning of future spacecraft missions require as inputs the F10.7 and Ap. The solar forecast is updated each month by executing MSAFE using historical and the latest month's observed solar indices to provide estimates for the balance of the current solar cycle. The forecasted solar indices represent the 13-month smoothed values consisting of a best estimate value stated as a 50 percentile value along with approximate +/- 2 sigma values stated as 95 and 5 percentile statistical values. This presentation will give an overview of the MSAFE model and the forecast for the current solar cycle.

Suggs, Ron↗

The MSFC Solar Activity Future Estimation (MSAFE) Model

The Natural Environments Branch of the Engineering Directorate at Marshall Space Flight Center (MSFC) provides solar cycle forecasts for NASA space flight programs and the aerospace community. These forecasts provide future statistical estimates of sunspot number, solar radio 10.7 cm flux (F10.7), and the geomagnetic planetary index, Ap, for input to various space environment models. For example, many thermosphere density computer models used in spacecraft operations, orbital lifetime analysis, and the planning of future spacecraft missions require as inputs the F10.7 and Ap. The solar forecast is updated each month by executing MSAFE using historical and the latest month's observed solar indices to provide estimates for the balance of the current solar cycle. The forecasted solar indices represent the 13-month smoothed values consisting of a best estimate value stated as a 50 percentile value along with approximate +/- 2 sigma values stated as 95 and 5 percentile statistical values. This presentation will give an overview of the MSAFE model and the forecast for the current solar cycle.

Suggs, Ron↗

Implications of Viking color data for evolution of the Amenthes Region, Mars

In the Eastern Hemisphere of Mars, where the hemispheric dichotomy is not overlain by relatively recent flows from the Tharsis volcanoes, the Mars cratered terrain boundary (CTB) now exists as highly fractured terrain juxtaposed with the smooth plains of the Northern Hemisphere. Based on Mariner 9 photography, the southern cratered terrain in the Amenthes region (0 to 30 deg N, 225 to 270 deg W) is approximately 3 to 4 km above the northern plains, and the transition zone is marked by both broad plateaus and knobby terrain. In the Amenthes Southest quadrangle, in particular, the continuum between large detached plateaus, smaller smooth topped plateaus and knobby hills can be seen. To define the previous position and evolution of the boundary, the distribution of cratered terrain, characteristic landforms (knobby terrain and detached plateaus) have been mapped, as well as structural features across the boundary. Structural mapping indicates that geological formations parallel the boundary in the Eastern Hemisphere with the conclusion that CTB evolution involved normal faulting through extension perpendicular to the boundary.

Maxwell, Ted A.↗

Process for Smoothing an Si Substrate after Etching of SiO2

A reactive-ion etching (RIE) process for smoothing a silicon substrate has been devised. The process is especially useful for smoothing those silicon areas that have been exposed by etching a pattern of holes in a layer of silicon dioxide that covers the substrate. Applications in which one could utilize smooth silicon surfaces like those produced by this process include fabrication of optical waveguides, epitaxial deposition of silicon on selected areas of silicon substrates, and preparation of silicon substrates for deposition of adherent metal layers. During etching away of a layer of SiO2 that covers an Si substrate, a polymer becomes deposited on the substrate, and the substrate surface becomes rough (roughness height approximately equal to 50 nm) as a result of over-etching or of deposition of the polymer. While it is possible to smooth a silicon substrate by wet chemical etching, the undesired consequences of wet chemical etching can include compromising the integrity of the SiO2 sidewalls and undercutting of the adjacent areas of the silicon dioxide that are meant to be left intact. The present RIE process results in anisotropic etching that removes the polymer and reduces height of roughness of the silicon substrate to less than 10 nm while leaving the SiO2 sidewalls intact and vertical. Control over substrate versus sidewall etching (in particular, preferential etching of the substrate) is achieved through selection of process parameters, including gas flow, power, and pressure. Such control is not uniformly and repeatably achievable in wet chemical etching. The recipe for the present RIE process is the following: Etch 1 - A mixture of CF4 and O2 gases flowing at rates of 25 to 75 and 75 to 125 standard cubic centimeters per minute (stdcm3/min), respectively; power between 44 and 55 W; and pressure between 45 and 55 mtorr (between 6.0 and 7.3 Pa). The etch rate lies between approximately equal to 3 and approximately equal to 6 nm/minute. Etch 2 - O2 gas flowing at 75 to 125 stdcm3/min, power between 44 and 55 W, and pressure between 50 and 100 mtorr (between 6.7 and 13.3 Pa).

Turner, Tasha↗

Electric Fields, Cloud Microphysics, and Reflectivity in Anvils of Florida Thunderstorms

A coordinated aircraft - radar project that investigated the electric fields, cloud microphysics and radar reflectivity of thunderstorm anvils near Kennedy Space Center is described. Measurements from two cases illustrate the extensive nature of the microphysics and electric field observations. As the aircraft flew from the edges of anvils into the interior, electric fields very frequently increased abruptly from approximately 1 to more than 10 kV m(exp -1) even though the particle concentration and radar reflectivity increased smoothly. The abrupt increase in field usually occurred when the aircraft entered regions with a reflectivity of 10 to 15 dBZ. It is suggested that the abrupt increase in electric field may be because the charge advection from the storm core did not occur across the entire breadth of the anvil and was not constant in time. Screening layers were not detected near the edges of the anvils. Some long-lived anvils showed subsequent enhancement of electric field and reflectivity and growth of particles, which if localized, might be a factor in explaining the abrupt change of field in some cases. Comparisons of electric field magnitude with particle concentration or reflectivity for a combined data set that included all anvil measurements showed a threshold behavior. When the average reflectivity, such as in a 3-km cube, was less than approximately 5 dBZ, the electric field magnitude was les than kV m(exp -1). Based on these findings, the Volume Averaged Height Integrated Radar Reflectivity (VAHIRR) is now being used by NASA, the Air Force and Federal Aviation Administration in new Lightning Launch Commit Criteria as a diagnostic for high electric fields in anvils.

Dye, J. E.↗

Intensity Oscillation of Low Energy Solar Particles Observed on the IMP 7 Satellite

Solar flare particles in the energy range from 125 to 160 KeV per charge were detected and analyzed by an electrostatic deflection vs. energy instrument on board the Explorer 47 (IMP 7) satellite. A harmonic analysis of the counting rates shows that frequently, the fluctuation of the particle flux can be approximately represented by N(t) = f(t)sin 2 pi nu t where f(t) is a smoothly varying function of t. The oscillation frequency nu varies from flare to flare and is of the order of 2 x .001/sec. The physical implications of the finding are discussed.

Fan, C. Y.↗

Minimal norm constrained interpolation

In computational fluid dynamics and in CAD/CAM, a physical boundary is usually known only discreetly and most often must be approximated. An acceptable approximation preserves the salient features of the data such as convexity and concavity. In this dissertation, a smooth interpolant which is locally concave where the data are concave and is locally convex where the data are convex is described. The interpolant is found by posing and solving a minimization problem whose solution is a piecewise cubic polynomial. The problem is solved indirectly by using the Peano Kernal theorem to recast it into an equivalent minimization problem having the second derivative of the interpolant as the solution. This approach leads to the solution of a nonlinear system of equations. It is shown that Newton's method is an exceptionally attractive and efficient method for solving the nonlinear system of equations. Examples of shape-preserving interpolants, as well as convergence results obtained by using Newton's method are also shown. A FORTRAN program to compute these interpolants is listed. The problem of computing the interpolant of minimal norm from a convex cone in a normal dual space is also discussed. An extension of de Boor's work on minimal norm unconstrained interpolation is presented.

Irvine, L. D.↗

Extending the Operational Envelope of a Turbofan Engine Simulation into the Sub-Idle Region

In many non-linear gas turbine simulations, operation in the sub-idle region can lead to model instability. This paper lays out a method for extending the operational envelope of a map based gas turbine simulation to include the sub-idle region. This method develops a multi-simulation solution where the baseline component maps are extrapolated below the idle level and an alternate model is developed to serve as a safety net when the baseline model becomes unstable or unreliable. Sub-idle model development takes place in two distinct operational areas, windmilling/shutdown and purge/cranking/ startup. These models are based on derived steady state operating points with transient values extrapolated between initial (known) and final (assumed) states. Model transitioning logic is developed to predict baseline model sub-idle instability, and transition smoothly and stably to the backup sub-idle model. Results from the simulation show a realistic approximation of sub-idle behavior as compared to generic sub-idle engine performance that allows the engine to operate continuously and stably from shutdown to full power.

aerothermodynamics↗

Extending the Operational Envelope of a Turbofan Engine Simulation into the Sub-Idle Region

In many non-linear gas turbine simulations, operation in the sub-idle region can lead to model instability. This paper lays out a method for extending the operational envelope of a map based gas turbine simulation to include the sub-idle region. This method develops a multi-simulation solution where the baseline component maps are extrapolated below the idle level and an alternate model is developed to serve as a safety net when the baseline model becomes unstable or unreliable. Sub-idle model development takes place in two distinct operational areas, windmilling/shutdown and purge/cranking/startup. These models are based on derived steady state operating points with transient values extrapolated between initial (known) and final (assumed) states. Model transitioning logic is developed to predict baseline model sub-idle instability, and transition smoothly and stably to the backup sub-idle model. Results from the simulation show a realistic approximation of sub-idle behavior as compared to generic sub-idle engine performance that allows the engine to operate continuously and stably from shutdown to full power.

gas path dynamics↗

Closed-Form Approximation of the Total Variation Proximal Operator

Total variation (TV) is a widely used function for regularizing imaging inverse problems that is particularly appropriate for images whose underlying structure is piecewise constant. TV regularized optimization problems are typically solved using proximal methods, but the way in which they are applied is constrained by the absence of a closed-form expression for the proximal operator of the TV function. A closed-form approximation of the TV proximal operator has previously been proposed, but its accuracy was not theoretically explored in detail. Here, we address this gap by making several new theoretical contributions, proving that the approximation leads to a proximal operator of some convex function, it is equivalent to a gradient descent step on a smoothed version of TV, and that its error can be fully characterized and controlled with its scaling parameter. We experimentally validate our theoretical results on image denoising and sparse-view computed tomography (CT) image reconstruction.

97 MATHEMATICS AND COMPUTING↗

Power Spectrum of Atmospheric Scintillation for the Deep Space Network Goldstone Ka-Band Downlink

Dynamic signal fluctuations due to atmospheric scintillations may impair the Ka-band (around 32-GHz) link sensitivities for a low-margin Deep Space Network (DSN) receiving system. The ranges of frequency and power of the fast fluctuating signals (time scale less than 1 min) are theoretically investigated using the spatial covariance and turbulence theory. Scintillation power spectrum solutions are derived for both a point receiver and a finite-aperture receiver. The aperture-smoothing frequency ((omega(sub s)), corner frequency ((omega(sub c)), and damping rate are introduced to define the shape of the spectrum for a finite-aperture antenna. The emphasis is put on quantitatively describing the aperture-smoothing effects and graphically estimating the corner frequency for a large aperture receiver. Power spectral shapes are analyzed parametrically in detail through both low- and high-frequency approximations. It is found that aperture-averaging effects become significant when the transverse correlation length of the scintillation is smaller than the antenna radius. The upper frequency or corner frequency for a finite-aperture receiver is controlled by both the Fresnel frequency and aperture-smoothing frequency. Above the aperture-smoothing frequency, the spectrum rolls off at a much faster rate of exp (-omega(sup 2)/omega(sup 2, sub s), rather than omega(sup -8/3), which is customary for a point receiver. However, a relatively higher receiver noise level can mask the fast falling-off shape and make it hard to be identified. We also predict that when the effective antenna radius a(sub r) less than or = 6 m, the corner frequency of its power spectrum becomes the same as that for a point receiver. The aperture-smoothing effects are not obvious. We have applied these solutions to the scenario of a DSN Goldstone 34-m-diameter antenna and predicted the power spectrum shape for the receiving station. The maximum corner frequency for the receiver (with omega(sub s) = 0.79 omega(sub 0) is found to be 0.44 Hz (or 1.0 omega(sub 0), while the fading rate (or fading slope) is about 0.06 dB/s.

Ho, C.↗

Collapse of primordial gas clouds and the formation of quasar black holes

The formation of quasar black holes during the hydrodynamic collapse of protogalactic gas clouds is discussed. The dissipational collapse and long-term dynamical evolution of these systems is analyzed using three-dimensional numerical simulations. The calculations focus on the final collapse stages of the inner baryonic component and therefore ignore the presence of dark matter. Two types of initial conditions are considered: uniformly rotating spherical clouds, and iirotational ellipsoidal clouds. In both cases the clouds are initially cold, homogeneous, and not far from rotational support (T/(absolute value of W) approximately equals 0.1). Although the details of the dynamical evolution depend sensitively on the initial conditions, the qualitative features of the final configurations do not. Most of the gas is found to fragment into small dense clumps, that eventually make up a spheroidal component resembling a galactic bulge. About 5% of the initial mass remains in the form of a smooth disk of gas supported by rotation in the gravitational potential potential well of the outer spheroid. If a central seed black hole of mass approximately greater than 10(exp 6) solar mass forms, it can grow by steady accretion from the disk and reach a typical quasar black hole mass approximately 10(exp 8) solar mass in less than 5 x 10(exp 8) yr. In the absence of a sufficiently massive seed, dynamical instabilities in a strongly self-gravitating inner region of the disk will inhibit steady accretion of gas and may prevent the immediate formation of quasar.

Loeb, Abraham↗

Can quasars photoionize the intergalactic medium at high redshift?

The reionization of the intergalactic medium (IGM) by quasar sources at high redshift are discussed. The integrated UV background from observed QSO's, taking into account the hydrogen opacity associated with intervening Ly-alpha clouds and Lyman limit systems are computed. It is noted that the published data appear to indicate a significant underdensity of absorption systems in the Ly-alpha forest with column densities N(sub HI) greater than 10(exp 15) cm(sup -2). This deficit results in a reduction of the opacity of the universe by a factor of 1.5-3 at z = 3-5 relative to previous estimates. The QSO contribution to the metagalactic flux at the Lyman edge may be as large as J(sub 912)(z) is approximately 6((1 + z)/4.5)(sup 0.5) x 10(exp -22) erg cm(sup -2) s(sup -1) Hz(sup -1) sr(sup -1) for q(sup o) = O, and slightly lower for q(sub o) = 1/2. For a density of the diffuse component of the IGM of omega(sub D)(h(sub 50)(sup 2)) less than 0.025, QSO's could photoionize a smooth IGM sufficiently to satisfy the constraints imposed by the Gunn-Peterson effect. The epoch of reionization could be as recent as z is approximately greater than 5. As a result, neutral patches of IGM would be detectable in the spectra of high redshift quasars. The patches would appear as absorption line systems with typical column densities of 10(exp 19) - 10(exp 20) cm(sup -2), and velocity widths of 100 - 1000 km s(sup -1).

Meiksin, Avery↗

TNet: A Model-Constrained Tikhonov Network Approach for Inverse Problems

Deep learning (DL), in particular deep neural networks, by default is purely data-driven and in general does not require physics. This is the strength of DL but also one of its key limitations when applied to science and engineering problems in which underlying physical properties—such as stability, conservation, and positivity—and accuracy are required. DL methods in their original forms are often not capable of respecting the underlying mathematical models or achieving desired accuracy even in big-data regimes. On the other hand, many data-driven science and engineering problems, such as inverse problems, typically have limited experimental or observational data, and DL would overfit the data in this case. Leveraging information encoded in the underlying mathematical models, we argue, not only compensates for missing information in low data regimes but also provides opportunities to equip DL methods with the underlying physics, hence promoting better generalization. This paper develops a model-constrained DL approach and its variant TNet—a Tikhonov neural network—which are capable of learning not only information hidden in the training data but also in the underlying mathematical models to solve inverse problems governed by partial differential equations in low data regimes. We provide the constructions and some theoretical results for the proposed approaches for both linear and nonlinear inverse problems. Since TNet is designed to learn inverse solutions with Tikhonov regularization, it is interpretable: in fact it recovers Tikhonov solutions for linear cases while potentially approximating Tikhonov solutions for nonlinear inverse problems. We also prove that data randomization can enhance not only the smoothness of the networks but also their generalizations. Comprehensive numerical results confirm the theoretical findings and show that with even as little as 1 training data sample for one-dimensional (1D) deconvolution, 5 for an inverse 2D heat conductivity problem, 100 for inverse initial conditions for a time-dependent 2D Burgers’s equation, and 50 for inverse initial conditions for 2D Navier–Stokes equations, TNet solutions can be as accurate as Tikhonov solutions while being several orders of magnitude faster. Furthermore, this is possible owing to the model-constrained term, replications, and randomization.

97 MATHEMATICS AND COMPUTING↗