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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

The forward scattering of microwave solar radiation from a water surface

The forward scattering of microwave solar radiation from both smooth and rough water surfaces is computed. The smooth surface is assumed specular, and the rough surface is represented by a two-scale surface, for which two small-scale perturbation parameters, 0.10 and 0.25, are considered. The contribution of the scattered sunlight to the antenna temperature is found using the scalar approximation, and the results are compared with radiometer measurements. The overall agreement is good, but in some cases the smooth-surface measurements are higher than the computations. This discrepancy possibly indicates an absolute calibration error or a slight misalignment of the antennas' boresights. The computations for the two perturbation parameters bracket the rough-surface measurements except when the sun's mirror image is far removed from the boresight direction. The small disagreement in this case may be due to a peaked large-scale slope distribution.

Wentz, F. J.↗
DOE OSTI2023

Calculation of Dynamical Response Functions Using a Bound-State Method

Abstract We investigate a method to extract response functions (dynamical polarisabilities) directly from a bound-state approach applied to calculations of perturbation-induced reactions. The use of a square-integrable basis leads to a response in the form of a sum of $$\delta $$ δ functions. We integrate this over energy and fit a smooth function to the resulting stepwise-continuous one. Its derivative gives the final approximation to the physical response function. We show that the method reproduces analytical results where known, and analyse the details for a variety of models. We apply it to some simple models, using the stochastic variational method as the numerical method. Albeit we find that this approach, and other numerical techniques, have some difficulties with the threshold behavior in coupled-channel problems with multiple thresholds, its stochastic nature allows us to extract robust results even for such cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Gaussian process hydrodynamics

Abstract We present a Gaussian process (GP) approach, called Gaussian process hydrodynamics (GPH) for approximating the solution to the Euler and Navier-Stokes (NS) equations. Similar to smoothed particle hydrodynamics (SPH), GPH is a Lagrangian particle-based approach that involves the tracking of a finite number of particles transported by a flow. However, these particles do not represent mollified particles of matter but carry discrete/partial information about the continuous flow. Closure is achieved by placing a divergence-free GP prior ξ on the velocity field and conditioning it on the vorticity at the particle locations. Known physics (e.g., the Richardson cascade and velocity increment power laws) is incorporated into the GP prior by using physics-informed additive kernels. This is equivalent to expressing ξ as a sum of independent GPs ξ l , which we call modes, acting at different scales (each mode ξ l self-activates to represent the formation of eddies at the corresponding scales). This approach enables a quantitative analysis of the Richardson cascade through the analysis of the activation of these modes, and enables us to analyze coarse-grain turbulence statistically rather than deterministically. Because GPH is formulated by using the vorticity equations, it does not require solving a pressure equation. By enforcing incompressibility and fluid-structure boundary conditions through the selection of a kernel, GPH requires significantly fewer particles than SPH. Because GPH has a natural probabilistic interpretation, the numerical results come with uncertainty estimates, enabling their incorporation into an uncertainty quantification (UQ) pipeline and adding/removing particles (quanta of information) in an adapted manner. The proposed approach is suitable for analysis because it inherits the complexity of state-of-the-art solvers for dense kernel matrices and results in a natural definition of turbulence as information loss. Numerical experiments support the importance of selecting physics-informed kernels and illustrate the major impact of such kernels on the accuracy and stability. Because the proposed approach uses a Bayesian interpretation, it naturally enables data assimilation and predictions and estimations by mixing simulation data and experimental data.

Mathematics↗
DOE OSTI2021

Solar Transposition Modeling via Deep Neural Networks With Sky Images

This article presents a notable advance toward the development of a new method of increasing the single-axis tracking photovoltaic (PV) system power output by improving the determination and near-term prediction of the optimum module tilt angle. The tilt angle of the plane receiving the greatest total irradiance changes with Sun position and atmospheric conditions including cloud formation and movement, aerosols, and particulate loading, as well as varying albedo within a module's field of view. In this article, we present a multi-input convolutional neural network that can create a profile of plane-of-array irradiance versus surface tilt angle over a full 180° arc from horizon to horizon. As input, the neural network uses the calculated solar position and clear-sky irradiance values, along with sky images. The target irradiance values are provided by the multiplanar irradiance sensor (MPIS). In order to account for varying irradiance conditions, the MPIS signal is normalized by the theoretical clear-sky global horizontal irradiance. Using this information, the neural network outputs an N -dimensional vector, where N is the number of points to approximate the MPIS curve via Fourier resampling. The output vector of the model is smoothed with a Gaussian kernel to account for error in the downsamping and subsequent upsampling steps, as well as to smooth the unconstrained output of the model. Furthermore, these profiles may be used to perform near-term prediction of angular irradiance, which can then inform the movement of a PV tracker.

42 ENGINEERING↗
NASA NTRS1989

Local error estimates for discontinuous solutions of nonlinear hyperbolic equations

Let u(x,t) be the possibly discontinuous entropy solution of a nonlinear scalar conservation law with smooth initial data. Suppose u sub epsilon(x,t) is the solution of an approximate viscosity regularization, where epsilon greater than 0 is the small viscosity amplitude. It is shown that by post-processing the small viscosity approximation u sub epsilon, pointwise values of u and its derivatives can be recovered with an error as close to epsilon as desired. The analysis relies on the adjoint problem of the forward error equation, which in this case amounts to a backward linear transport with discontinuous coefficients. The novelty of this approach is to use a (generalized) E-condition of the forward problem in order to deduce a W(exp 1,infinity) energy estimate for the discontinuous backward transport equation; this, in turn, leads one to an epsilon-uniform estimate on moments of the error u(sub epsilon) - u. This approach does not follow the characteristics and, therefore, applies mutatis mutandis to other approximate solutions such as E-difference schemes.

Tadmor, Eitan↗
NASA NTRS1991

High frequency scattering from corrugated stratified cylinders

Interest in applying radar remote sensing for the study of forested areas led to the development of a model for scattering from corrugated stratified dielectric cylinders. The model is used to investigate the effect of bark and its roughness on scattering from tree trunks and branches. The outer layer of the cylinder (bark) is assumed to be a low-loss dielectric material and to have a regular (periodic) corrugation pattern. The inner layers are treated as lossy dielectrics with smooth boundaries. A hybrid solution based on the moment method and the physical optics approximation is obtained. In the solution, the corrugations are replaced with polarization currents that are identical to those of the local tangential periodic corrugated surface, and the stratified cylinder is replaced with equivalent surface currents. New expressions for the equivalent physical-optics currents are used which are more convenient than the standard ones. It is shown that the bark layer and its roughness both reduce the radar cross-section. It is also demonstrated that the corrugations can be replaced by an equivalent anisotropic layer.

Sarabandi, Kamal↗
NASA NTRS1995

On High-Order Radiation Boundary Conditions

In this paper we develop the theory of high-order radiation boundary conditions for wave propagation problems. In particular, we study the convergence of sequences of time-local approximate conditions to the exact boundary condition, and subsequently estimate the error in the solutions obtained using these approximations. We show that for finite times the Pade approximants proposed by Engquist and Majda lead to exponential convergence if the solution is smooth, but that good long-time error estimates cannot hold for spatially local conditions. Applications in fluid dynamics are also discussed.

Hagstrom, Thomas↗
NASA NTRS1997

Dynamic Forms: Application to Aircraft Guidance - Part 2

The paper describes a method for guiding a dynamic system through a given set of points. The paradigm is a fully automatic aircraft subject to air traffic control (ATC). The ATC provides a sequence of waypoints through which the aircraft trajectory must pass. The waypoints typically specify time, position, and velocity. The guidance problem is to synthesize a system state trajectory that satisfies both the ATC and aircraft constraints. Complications arise because the controlled process is multidimensional, multiaxis, nonlinear, highly coupled, and the state space is not flat. In addition, there is a multitude of operating modes, which may number in the hundreds. Each such mode defines a distinct state space model of the process by specifying the state space coordinatization, the partition of the controls into active controls and configuration controls, and the output map. Furthermore, mode transitions are required to be smooth. The proposed guidance algorithm is based on the inversion of the pure feedback approximation, followed by correction for the effects of zero dynamics. The paper describes the structure and major modules of the algorithm, and the performance is illustrated by several example aircraft maneuvers.

Meyer, George↗
NASA NTRS2002

Tackling a Hot Paradox: Laminar Soot Processes-2 (LSP-2)

The last place you want to be in traffic is behind the bus or truck that is belching large clouds of soot onto your freshly washed car. Besides looking and smelling bad, soot is a health hazard. Particles range from big enough to see to microscopic and can accumulate in the lungs, potentially leading to debilitating or fatal lung diseases. Soot is wasted energy, and therein lies an interesting paradox: Soot forms in a flame's hottest regions where you would expect complete combustion and no waste. Soot enhances the emissions of other pollutants (carbon monoxide and polyaromatic hydrocarbons, etc.) from flames and radiates unwanted heat to combustion chambers (a candle's yellowish glow is soot radiating heat), among other effects. The mechanisms of soot formation are among the most important unresolved problems of combustion science because soot affects contemporary life in so many ways. Although we have used fire for centuries, many fundamental aspects of combustion remain elusive, in part because of limits imposed by the effects of gravity on Earth. Hot or warm air rises quickly and draws in fresh cold air behind it, thus giving flames the classical teardrop shape. Reactions occur in a very small zone, too fast for scientists to observe, in detail, what is happening inside the flame. The Laminar Soot Processes (LSP-2) experiments aboard STS-107 will use the microgravity environment of space to eliminate buoyancy effects and thus slow the reactions inside a flame so they can be more readily studied. 'Laminar' means a simple, smooth fuel jet burning in air, somewhat like a butane lighter. This classical flame approximates combustion in diesel engines, aircraft jet propulsion engines, and furnaces and other devices. LSP-2 will expand on surprising results developed from its first two flights in 1997. The data suggest the existence of a universal relationship, the soot paradigm, that, if proven, will be used to model and control combustion systems on Earth. STS-107 experiments also will help set the stage for extended combustion experiments aboard the International Space Station.

Faeth, Gerard M.↗
NASA NTRS1995

Validation of Global Climatologies of Trace Gases Using NASA Global Tropospheric Experiment (GTE) Data

Methane is an important trace gas because it is a greenhouse gas that affects the oxidative capacity of the atmosphere. It is produced from biological and anthropogenic sources, and is increasing globally at a rate of approximately 0.6% per year [Climate Change 1992, IPCC]. By using National Oceanic and Atmospheric Administration/Climate Monitoring and Diagnostics Laboratory (NOAA/CMDL) ground station data, a global climatology of methane values was produced. Unfortunately, because the NOAA/CMDL ground stations are so sparse, the global climatology is low resolution. In order to compensate for this low resolution data, it was compared to in-situ flight data obtained from the NASA Global Tropospheric Experiment (GTE). The smoothed ground station data correlated well with the flight data. Thus, for the first time it is shown that the smoothing process used to make global contours of methane using the ground stations is a plausible way to approximate global atmospheric concentrations of the gas. These verified climatologies can be used for testing large-scale models of chemical production, destruction, and transport. This project develops the groundwork for further research in building global climatologies from sparse ground station data and studying the transport and distribution of trace gases.

Courchaine, Brian↗
NASA NTRS1997

Torque Reversal and Spin-Down of the Accretion-Powered Pulsar 4U 1626-67

Over 5 yr of hard X-ray (20-60 keV) monitoring of the 7.66 s accretion-powered pulsar 4U 1626-67 with the Compton Gamma Ray Observatory/BATSE large-area detectors has revealed that the neutron star is now steadily spinning down, in marked contrast to the steady spin-up and spin-down torques differ by only 15% with the neutron star spin changing on a timescale |v/v| approximately equals 5000 yr in both states. The current spin-down rate is itself decreasing on a timescale |v/v| approximately equals 26 yr. The long-term timing history shows small-amplitude variations on a 4000 day timescale, which are probably due to variations in the mass transfer rate. The pulsed 20-60 keV emission from 4U 1626-67 is well-fitted by a power-law spectrum with photon index gamma = 4.9 and a typical pulsed intensity of 1.5 x 10(exp -10) ergs cm (exp -2)s(exp -1). The low count rates with BATSE prohibited us from constraining the reported 42 minute binary orbit, but we can rule out long-period orbits in the range 2 days < or = P(orb) < or = 900 days. We compare the long-term torque behavior of 4U 1626-67 to other disk-fed accreting pulsars and discuss the implications of our results for the various theories of magnetic accretion torques. The abrupt change in the sign of the torque is difficult to reconcile with the extremely smooth spin-down now observed. The strength of the torque noise in 4U 1626-67, approximately 10(exp -22) Hz(exp 2)s(exp -2) Hz(exp -1), is the smallest ever measured for an accreting X-ray pulsar, and it is comparable to the timing noise seen in young radio pulsars. We close by pointing out that the core temperature and external torque (the two parameters potentially relevant to internal sources of timing noise) of an accreting neutron star are also comparable to those of young radio pulsars.

Chakrabarty, Deepto L.↗
DOE OSTI2024

Image-Based Digital Twin for Assessing the Coupled Electro-Chemo-Mechanical Behavior of Li-Ion Batteries

In energy storage materials, strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking during charge/discharge cycling, resulting in reduced performance and service life. In this work, a digital twin is created to investigate the performance of a Li-ion battery cathode and simulate degradation accumulation. Pixel-based model construction is used to represent the complex material geometries from microstructural images supplied by the National Renewable Energy Laboratory (NREL). Because of the expected large deformation and crack opening, the reproducing kernel particle method (RKPM), a meshfree method with discretization at the image pixels, is used to approximate the field variables: electrostatic potential, concentration, and displacement. An interface modified reproducing kernel (IM-RK) is constructed by scaling a smooth kernel function with an interface-distance function to achieve strategic discontinuity types (i.e. weak discontinuities for strain discontinuities and strong discontinuities for cracks) and alleviate Gibbs oscillations near these transition zones. The discrete nature of the images' pixel points is employed throughout the model and approximation construction. The mechanical model is verified using an image-based microstructure under tensile loading. A transient electro-chemo-mechanical coupled simulation is performed to evaluate the potential micro-cracking induced degradation of the battery cathode material subjected to charge/discharge cycling.

battery degradation↗
DOE OSTI2021

Source term method for binary neutron stars initial data

The initial condition problem for a binary neutron star system requires a Poisson equation solver for the velocity potential with a Neumann-like boundary condition on the surface of the star. Difficulties that arise in this boundary value problem are: (a) the boundary is not known a priori, but constitutes part of the solution of the problem; (b) various terms become singular at the boundary. In this work, we present a new method to solve the fluid Poisson equation for irrotational/spinning binary neutron stars. The advantage of the new method is that it does not require complex fluid surface fitted coordinates and it can be implemented in a Cartesian grid, which is a standard choice in numerical relativity calculations. This is accomplished by employing the source term method proposed by Towers, where the boundary condition is treated as a jump condition and is incorporated as additional source terms in the Poisson equation, which is then solved iteratively. The issue of singular terms caused by vanishing density on the surface is resolved with an additional separation that shifts the computation boundary to the interior of the star. We present two-dimensional tests to show the convergence of the source term method, and we further apply this solver to a realistic three-dimensional binary neutron star problem. By comparing our solution with the one coming from the initial data solver cocal, we demonstrate agreement to approximately 1%. We report our method can be used in other problems with non-smooth solutions like in magnetized neutron stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Applications of the modified Hulthén-Kohn method for bound and scattering states

We adapt the Hulthén–Kohn method suggested by Efros [Phys. Rev. C 99, 034620 (2019)] for calculating various observables in the continuum and discrete spectrum using two-body interactions in single- and coupled channel systems. We explore the convergence of phase shifts and wave functions as well as the location of S-matrix poles which enables obtaining both resonance and bound state parameters. We find that employing a harmonic oscillator basis, together with an interaction smoothing scheme introduced by Gyarmati et al. [Nucl. Phys. A 326, 119 (1979)], and adopting approximate bound-state solutions for the short-range components of basis wave functions lead to good convergence even with restricted oscillator quanta accessible for modern no-core shell model codes. The adapted Efros method will facilitate ab initio many-body nuclear structure applications.

Nuclear reactions↗
DOE OSTI2022

Generalized Bootstrap AMG and AIR-AMG for coupled PDE systems with a focus on spacetime discretizatoins

The Pennsylvania State University (“Subcontractor”) has worked on the design of new algebraic, parallel, multilevel methods that obtain the full space and time solution of systems of PDEs. In particular, the PI and his collaborators explored semi-intrusive approaches based on algebraic multigrid (AMG). The focus of the research has been on the development of these techniques for the Euler equations in 1d and 2d. The overall research focused on the development of adaptive AIR (approximate ideal restriction) AMG solvers for these problems. The PI also explored the use of smoothed aggregation and root-node energy-based AMG solvers for these problems.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

Marching iterative methods for the parabolized and thin layer Navier-Stokes equations

Downstream marching iterative schemes for the solution of the Parabolized or Thin Layer (PNS or TL) Navier-Stokes equations are described. Modifications of the primitive equation global relaxation sweep procedure result in efficient second-order marching schemes. These schemes take full account of the reduced order of the approximate equations as they behave like the SLOR for a single elliptic equation. The improved smoothing properties permit the introduction of Multi-Grid acceleration. The proposed algorithm is essentially Reynolds number independent and therefore can be applied to the solution of the subsonic Euler equations. The convergence rates are similar to those obtained by the Multi-Grid solution of a single elliptic equation; the storage is also comparable as only the pressure has to be stored on all levels. Extensions to three-dimensional and compressible subsonic flows are discussed. Numerical results are presented.

Israeli, M.↗
NASA NTRS1985

Marching iterative methods for the Parabolized and Thin Layer Navier-Stokes equations

Downstream marching iterative schemes for the solution of the Parabolized or Thin Layer (PNS or TL) Navier-Stokes equations are described. Modifications of the primitive equation global relaxation sweep procedure result in efficient second-order marching schemes. These schemes take full account of the reduced order of the approximate equations as they behave like the SLOR for a single elliptic equation. The improved smoothing properties permit the introduction of Multi-Grid acceleration. The proposed algorithm is essentially Reynolds number independent and therefore can be applied to the solution of the subsonic Euler equations. The convergence rates are similar to those obtained by the Multi-Grid solution of a single elliptic equation; the storage is also comparable as only the pressure has to be stored on all levels. Extensions to three-dimensional and compressible subsonic flows are discussed. Numerical results are presented.

Israeli, M.↗