Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Weight function”

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 127 records · Page 7

Quadrature imposition of compatibility conditions in Chebyshev methods

Often, in solving an elliptic equation with Neumann boundary conditions, a compatibility condition has to be imposed for well-posedness. This condition involves integrals of the forcing function. When pseudospectral Chebyshev methods are used to discretize the partial differential equation, these integrals have to be approximated by an appropriate quadrature formula. The Gauss-Chebyshev (or any variant of it, like the Gauss-Lobatto) formula can not be used here since the integrals under consideration do not include the weight function. A natural candidate to be used in approximating the integrals is the Clenshaw-Curtis formula, however it is shown that this is the wrong choice and it may lead to divergence if time dependent methods are used to march the solution to steady state. The correct quadrature formula is developed for these problems. This formula takes into account the degree of the polynomials involved. It is shown that this formula leads to a well conditioned Chebyshev approximation to the differential equations and that the compatibility condition is automatically satisfied.

Gottlieb, D.↗

Streamline-upwind/Petrov-Galerkin procedures for the vorticity-stream function form of the Navier-Stokes equations

The paper presents procedures for the solution of the Navier-Stokes equations in the vorticity-stream function form. The difficulties involved are related to the convection term in the vorticity transport equation and to the lack of boundary conditions for voritcity at no-slip surfaces. The use of a streamline-upwind/Petrov-Galerkin finite element formulation for the solution of the vorticity transport equation. In the present scheme, the weighting functions are dependent on both spatial and temporal discretizations. A proper numerical treatment of the boundary conditions leading to an implicit treatment of the vorticity at no-slip surfaces is presented. These procedures have successfully been employed to simulate various flows of engineering interest.

Tezduyar, T. E.↗

Finite element algorithms for compressible flow computation on a supercomputer

Methods of applying computational fluid dynamics techniques to engineering problems are discussed. The linking factors between theoretical and applied research are the shape and weight functions which produce the Galerkin finite element schemes. The flow is unsteady, viscous, heat conducting and compressible, and steady-state flows are the asymptotic limit of unsteady flows. All flows are described by the time-averaged Navier-Stokes equations (NSE) with mass and energy conservation. Local curvilinear intrinsic coordinates are applied to discretize the NSE in arbitrary geometric domains. A modified weighted residuals approach defines the discrete analogs of the physical systems modeled and a general interpolants method is used to derive families of numerical models, both implicit and explicit, finite difference and finite elemnt, from a single point of departure. Time-dependent element approximations are achieved with a progressive assembly of generalized elements method. Hyperbolic steady-state Euler algorithms and quasi-parabolic are spatial marching algorithms for solving the discretized equations. Application of the techniques is illustrated through calculation of the compressible flow around a hypersonic flight vehicle at 100,000 ft at zero angle of attack at speeds of Mach 4-10.

Spradley, L. W.↗

Adaptive Detector Arrays for Optical Communications Receivers

The structure of an optimal adaptive array receiver for ground-based optical communications is described and its performance investigated. Kolmogorov phase screen simulations are used to model the sample functions of the focal-plane signal distribution due to turbulence and to generate realistic spatial distributions of the received optical field. This novel array detector concept reduces interference from background radiation by effectively assigning higher confidence levels at each instant of time to those detector elements that contain significant signal energy and suppressing those that do not. A simpler suboptimum structure that replaces the continuous weighting function of the optimal receiver by a hard decision on the selection of the signal detector elements also is described and evaluated. Approximations and bounds to the error probability are derived and compared with the exact calculations and receiver simulation results. It is shown that, for photon-counting receivers observing Poisson-distributed signals, performance improvements of approximately 5 dB can be obtained over conventional single-detector photon-counting receivers, when operating in high background environments.

Vilnrotter, V.↗

A solution-adaptive mesh algorithm for dynamic/static refinement of two and three dimensional grids

An adaptive grid algorithm has been developed in two and three dimensions that can be used dynamically with a solver or as part of a grid refinement process. The algorithm employs a transformation from the Cartesian coordinate system to a general coordinate space, which is defined as a parallelepiped in three dimensions. A weighting function, independent for each coordinate direction, is developed that will provide the desired refinement criteria in regions of high solution gradient. The adaptation is performed in the general coordinate space and the new grid locations are returned to the Cartesian space via a simple, one-step inverse mapping. The algorithm for relocation of the mesh points in the parametric space is based on the center of mass for distributed weights. Dynamic solution-adaptive results are presented for laminar flows in two and three dimensions.

Benson, Rusty A.↗

Effects of CELF Pretreatment Severity on Lignin Structure and the Lignin-Based Polyurethane Properties

Conversion of technical lignin into performance biopolymers such as polyurethane offers environmental and economic advantages when combined with production of biofuels from biomass sugars, presenting significant interest toward studying the role of pretreatment on lignin structure and functionality. Co-solvent enhanced lignocellulosic fractionation (CELF) pretreatment, employing acidic aqueous tetrahydrofuran (THF) mixtures, was developed to effectively break down the lignin-carbohydrate matrix and promote extraction of lignin from lignocellulosic biomass with desirable purity and yield. In this study, we report the effects of CELF pretreatment reaction severity on the molecular structure of CELF-extracted lignin and its impact toward the mechanical properties of the resulting lignin-based polyurethanes. Reaction temperature was found to play the most significant role, compared to reaction time and acidity, in manipulating structural features such as molecular weight, functionality and intra-polymer structure. At the severe reaction conditions at 180°C, the order of reactivity for primary lignin interlinkages characterized by semiquantitative HSQC NMR analysis were found to be β-ether > phenylcoumaran (β−5′) > resinol (β−β′) facilitating a high degree of depolymerization and yielding a high frequency of free phenolics and reduced aliphatic hydroxyl groups. All side-chain interlinkages were depleted converting guaiacyl subunits into condensed forms, while still retaining uncondensed syringyl subunits. Under the mild 150°C temperature reaction, CELF lignin had higher molecular weight and retained more β-ether interlinkages. The results from CELF lignin-based polyurethane synthesis indicated that the tensile properties depended on the miscibility of CELF lignin with other components and low molecular weight cuts improved the dispersion of lignin in the polyurethane network. Pre-mixing of CELF with poly(ethylene glycol) (PEG) reduced the brittleness and improved the ductility of the CELF lignin-PEG polyurethanes.

59 BASIC BIOLOGICAL SCIENCES↗

Optimizing X-Ray Optical Prescriptions for Wide-Field Applications

X-ray telescopes with spatial resolution optimized over the field of view (FOV) are of special interest for missions, such as WFXT, focused on moderately deep and deep surveys of the x-ray sky, and for solar x-ray observations. Here we report on the present status of an on-going study of the properties of Wolter I and polynominal grazing incidence designs with a view to gain a deeper insight into their properties and simply the design process. With these goals in mind, we present some results in the complementary topics of (1) properties of Wolter I x-ray optics and polynominal x-ray optic ray tracing. Of crucial importance for the design of wide-field x-ray optics is the optimization criteria. Here we have adopted the minimization of a merit function, M, which measures the spatial resolution averaged over the FOV: M= ((integral of d phi) between the limits of 0 and 2 pi) (integral of d theta theta w(theta) sigma square (theta,phi) between the limits of 0 and theta(sub FOV)) (integral of d phi between the limits of 0 and phi/4) (Integral of d theta theta w(theta) between the limits of 0 and theta(sub FOV) where w(theta(sub 1) is a weighting function and Merit function: sigma-square (theta, phi) = summation of (x,y,z) [ -<(x,y,z)> (exp 2)] is the spatial variance for a point source on the sky at polar and azimuthal off-axis angles (theta,phi).

Elsner, R. F.↗

Computation of averaged monthly zonal albedo utilizing the solar zenith angle, properties of clear and cloudy atmospheres

The zonal temporal averages of albedos at the top of the atmosphere were considered as a function of the length of the day. The length of the day were used to determine the average daily values of mu sub 0(=Cos of the solar zenith angle, theta sub 0). Polynominal fits of the slope and intercept functions of A sub s (cloud-free albedo) and A sub c(cloud albedo) as function of Cos theta sub 0 were obtained by using the sample values of albedo corresponding to solar zenith angles from 0 to 90 deg with interval of 5 deg. The daily zonal values of mu sub 0 and the surface albedos were used to compute the daily zonal values of albedos at the top of the clear and cloudy atmosphere. The monthly zonal cloud fractions were used to compute planetary albedo A at the top of the atmosphere. The global values of monthly albedos A sub s, A sub c and A were computed by using the weighting function defined as the difference of the sins of zonal values of latitudes. The computer program implementation is also described.

Dhuria, H.↗

Spatial and Temporal Variability of Carbon Dioxide Using Structure Functions in Urban Areas: Insights for Future Active Remote CO2 Sensors

High resolution in-situ CO2 measurements were recorded onboard the NASA P-3B during the DISCOVER-AQ (Deriving Information on Surface Conditions from Column and Vertically Resolved Observations Relevant to Air Quality) Field Campaigns during July 2011 over Washington DC/Baltimore, MD; January-February 2013 over the San Joaquin Valley, CA; September 2013 over Houston, TX; and July-August 2014 over Denver, CO. Each of these campaigns have approximately two hundred vertical soundings of CO2 within the lower troposphere (surface to about 5 kilometers) at 6-8 different sites in each of the urban areas. In this study, we used structure function analysis, which is a useful way to quantify spatial and temporal variability, by displaying differences with average observations, to evaluate the variability of CO2 in the 0-2 kilometers range (representative of the planetary boundary layer). These results can then be used to provide guidance in the development of science requirements for the future ASCENDS (Active Sensing of CO2 Emissions over Nights, Days, and Seasons) mission to measure near-surface CO2 variability in different urban areas. We also compare the observed in-situ CO2 variability with the variability of the CO2 column-averaged optical depths in the 0-1 kilometer and 0-3.5 kilometers altitude ranges in the four geographically different urban areas, using vertical weighting functions for potential future ASCENDS lidar CO2 sensors operating in the 1.57 and 2.05 millimeter measurement regions. In addition to determining the natural variability of CO2 near the surface and in the column, radiocarbon method using continuous CO2 and CO measurements are used to examine the variation of emission quantification between anthropogenic and biogenic sources in the DC/Maryland urban site.

Choi, Yonghoon↗

Progress in modelling fast-ion D-alpha spectra and neutral particle analyzer fluxes using FIDASIM

FIDASIM is a code that models signals produced by charge-exchange reactions between neutrals and ions (both fast and thermal) in magnetically confined plasmas. With the ion distribution function as input, the code predicts the efflux to a neutral particle analyzer diagnostic and the photon radiance of Balmer-alpha light to a fast-ion D α diagnostic, in addition to many other related quantities. A new, parallelized version of the Monte Carlo code FIDASIM has been developed in Fortran90 that is substantially faster than the original interactive data language version. Modified algorithms include more accurate treatments of the time dependent collisional-radiative equations that describe neutral energy levels, of the cloud of ‘halo’ neutrals that surround the injected neutral beam, and of finite Larmor radius effects. Enhanced physics capabilities include modelling ‘passive’ signals from cold edge neutrals, the ability to treat general three-dimensional magnetic confinement configurations, and calculations of diagnostic-specific weight functions that enable tomographic reconstructions of the fast-ion distribution function. Neutral beam attenuation, beam emission, and fast-ion birth profiles are also modelled. Finally, the new algorithms have been successfully validated against experimental data and new features have been tested through benchmarks between two independently developed versions of the code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evaluation of Vertically Resolved Water Winds from AIRS using Hurricane Katrina

The knowledge of wind velocity as a function of altitude is key to weather forecast improvements. The ability of hyperspectral sounders in principle to measure vertically resolved water winds, which has long been recognized, has been tested with Atmospheric Infrared Sounder (AIRS) data. AIRS retrievals of total column water above 300 mb have been correlated with the radiosonde upper-tropospheric wind velocity and moisture data. The excellent correlation is illustrated with results obtained from hurricane Katrina and from the western United States. AIRS is a hyperspectral infrared sounder in low Earth orbit. It was launched in May 2002. We illustrate the use of AIRS data for the measurement of upper tropospheric water by using the 2387/cm CO2 R-branch channel and the 1551/cm water vapor channel. The 2387/cm channel measures the temperature at 300 mb totally independent of water vapor. The weighting function of the 1551/cm channel peaks at 300 mb only under moist conditions; the peak shifts downward (higher temperature) for less water and upward (lower temperature) for more water. The difference between the brightness temperatures bt2387 and bt1551 cancels the local several degree weather related variability of the temperature and measures the component due to the water vapor at 300 mb.

Atmospheric Infrared Sounder (AIRS)↗

Electrical Cable Design for Urban Air Mobility Aircraft

Urban Air Mobility (UAM) describes a new type of aviation focused on efficient flight within urban areas for moving people and goods. There are many different configurations of UAM vehicles, but they generally use an electric motor driving a propeller or ducted fan powered by batteries or a hybrid electric power generation system. Transmission cables are used to move energy from the storage or generation system to the electric motors. Though terrestrial power transmission cables are well established technology, aviation applications bring a whole host of new design challenges that are not typical considerations in terrestrial applications. Aircraft power transmission cable designs must compromise between resistance-per-length, weight-per-length, volume constraints, and other essential qualities. In this paper we use a multidisciplinary design optimization to explore the sensitivity of these qualities to a representative tiltwing turboelectric UAM aircraft concept. This is performed by coupling propulsion and thermal models for a given mission criteria. Results presented indicate that decreasing cable weight at the expense of increasing cable volume or cooling demand is effective at minimizing maximum takeoff weight (MTO). These findings indicate that subsystem designers should update their modeling approach in order to contribute to system-level optimality for highly-coupled novel aircraft. Mobility (UAM) vehicles have the potential to change urban and intra-urban transport in new and interesting ways. In a series of two papers Johnson et al.1 and Silva et al.2 presented four reference vehicle configurations that could service different niches in the UAM aviation category. Of those, this paper focuses on the Vertical Take-off and Landing (VTOL) tiltwing configuration shown in Figure 1. This configuration uses a turboelectric power system, feeding power from a turbo-generator through a system of transmission cables to four motors spinning large propellers on the wings. Previous work on electric cable subsystems leaves much yet to be explored, especially in the realm of subsystem coupling. Several aircraft optimization studies1, 3, 4 only considered aircraft electrical cable weight and ignored thermal effects. Electric and hybrid-electric aircraft studies by Mueller et al.5 and Hoelzen et al.6 selected a cable material but did not investigate alternative materials. Advanced cable materials have been examined by a number of authors: Alvarenga7 examined carbon nanotube (CNT) conductors for low-power applications. De Groh8, 9 examined CNT conductors for motor winding applications. Behabtu et al.,10 and Zhao et al.11 examined CNT conductors for a general applications. There were some studies that examined the thermal effects of cables but they did not allow the cable material to change; El-Kady12 optimized ground-cable insulation and cooling subject constraints. Vratny13 selected cable material based on vehicle power demand, and required resulting cable heat to be dissipated by the Thermal Management System (TMS). None of these previous studies allowed for the selection of the cable material based on a system level optimization goal. Instead, they focused on sub-system optimality such as minimum weight, which comes at the expense of incurring additional costs for other subsystems. Dama14 selected overhead transmission line materials using a weighting function and thermal constraints. However, that work was not coupled with any aircraft subsystems like a TMS. The traditional aircraft design approach, which relies on assembling groups of optimal subsystems, breaks down when considering novel aircraft concepts like the tiltwing vehicle. In a large part, this is because novel concepts have a much higher degree of interaction or coupling between subsystems. For example, when a cable creates heat, this heat needs to be dissipated by the TMS, which needs power supplied by the turbine, and delivering the power creates more heat. The cable, the TMS, and the turbine are all coupled. A change to one subsystem will affect all the other subsystems, much to the consternation of subsystem design experts. Multidisciplinary optimization is the design approach that can address these challenges. However, to fully take advantage of this, we must change the way we think about subsystem design. Specifically, we must move away from point design, and focus on creating solution spaces. The work presented in this paper uses the multidisciplinary optimization approach with aircraft level models to study the system-level sensitivity of cable traits: weight-per-length and resistance-per-length. Additionally, we examined the effects of vehicle imposed volume constraints on these traits. This is useful for three purposes: (1) to demonstrate a framework that can perform a coupled analysis between the aircraft thermal and propulsion systems, (2) to provide a method by which future cable designs can be evaluated against each other given a system-level design goal, (3) to provide insight into what cable properties may be promising for future research. This last element is explored given the caveat that the models contained in this analysis do not represent high-fidelity systems. Thus, while we can demonstrate coupling in between systems, the exact system-level sensitivity to a given parameter may change if a subsystem model or the assumptions governing that model change. The organization of this paper is as follows, in Sec II we outline a method to combine the VTOL vehicle design and cable information in order to produce cables sensitivity studies. Results analysis and discussion are contained in Sec III. Conclusions are presented in Sec IV.

Aretskin-Hariton, Eliot D.↗

Instability and treatments of the coupled discrete element and lattice Boltzmann method by the immersed moving boundary scheme

The immersed moving boundary (IMB) scheme has been extensively used to couple the discrete element method (DEM) with the lattice Boltzmann method (LBM). In the literature, only the formulation of IMB for lattice nodal cells covered by a single-solid particle was given. The treatment of situations where a nodal cell is covered by two or more solid particles is seldom discussed. It is found that some numerical instability can occur for such situations due to an inappropriate computation of the weighting function in the IMB formulation. This work presents an enhanced treatment that can resolve the issue and validates it using some benchmark tests. Furthermore, to avoid the extra costs associated with the treatment and simplify the complicated procedure introduced, a simplified IMB scheme is proposed. The accuracy of both enhanced and simplified IMB schemes are validated by test cases including single-particle sedimentation, two-particle drafting-kissing-tumbling phenomenon, and multiple-particle sedimentation. Then, the robustness of both schemes is examined and discussed using a specially designed flow past cylinders test. Overall, the simplified IMB scheme is proved to be robust and sufficiently accurate and simpler and more effective than the enhanced scheme.

42 ENGINEERING↗

Non-uniform active learning for Gaussian process models with applications to trajectory informed aerodynamic databases

The ability to non-uniformly weight the input space is desirable for many applications, and has been explored for space-filling approaches. Increased interests in linking models, such as in a digital twinning framework, increases the need for sampling emulators where they are most likely to be evaluated. In particular, here we apply non-uniform sampling methods for the construction of aerodynamic databases. This paper combines non-uniform weighting with active learning for Gaussian Processes (GPs) to develop a closed-form solution to a non-uniform active learning criterion. We accomplish this by utilizing a kernel density estimator as the weight function. We demonstrate the need and efficacy of this approach with an atmospheric entry example that accounts for both model uncertainty as well as the practical state space of the vehicle, as determined by forward modeling within the active learning loop.

42 ENGINEERING↗

A mass-transfer particle-tracking method for simulating transport with discontinuous diffusion coefficients

The problem of a spatially discontinuous diffusion coefficient (D(x)) is one that may be encountered in hydrogeologic systems due to natural geological features or as a consequence of numerical discretization of flow properties. To date, mass-transfer particle-tracking (MTPT) methods, a family of Lagrangian methods in which diffusion is jointly simulated by random walk and diffusive mass transfers, have been unable to solve this problem. This manuscript presents a new mass-transfer (MT) algorithm that enables MTPT methods to accurately solve the problem of discontinuous D(x). To achieve this, we derive a semi-analytical solution to the discontinuous D(x) problem by employing a predictor-corrector approach, and we use this semi-analytical solution as the weighting function in a reformulated MT algorithm. Furthermore, this semi-analytical solution is generalized for cases with multiple 1D interfaces as well as for 2D cases, including a 2 × 2 tiling of 4 subdomains that corresponds to a numerically-generated diffusion field. The solutions generated by this new mass-transfer algorithm closely agree with an analytical 1D solution or, in more complicated cases, trusted numerical results, demonstrating the success of our proposed approach.

42 ENGINEERING↗

Pan-tropical daily L-band microwave land surface emissivity retrieval from GNSS-R observations

The uncertainties in microwave land surface emissivity (MLSE) measurements have long limited the use of spaceborne microwave radiometer data. As an emerging observation method, Global Navigation Satellite System Reflectometry (GNSS-R) has demonstrated great potential in several land and ocean applications. In this study, a method for obtaining daily MLSE dataset in the pan-tropical region from Cyclone GNSS (CYGNSS) observations is presented and evaluated. The CYGNSS observations are first aggregated into the Equal-Area-Scalable-Earth (EASE) 2.0 36 km grid by a combined weight function of distance, time, and signal-to-noise ratio variance. Then, the method employs a pixel-by-pixel regression algorithm to conduct the daily MLSE retrieval using reference emissivity derived from the Soil Moisture Active Passive (SMAP) brightness temperature. The CYGNSS MLSE shows good agreement with SMAP MLSE, delivering an overall root-mean-square error (RMSE) of 0.022 and 0.017 for horizontal and vertical polarization, respectively, during the training set spanning the whole year of 2018. Furthermore, on the test set from January 2019 to May 2019, the RMSE values amounted to 0.030 and 0.023 for horizontal and vertical polarization, respectively. Temperature records from the International Soil Moisture Network are employed to calculate the emissivity and for in-situ validation, which yield an RMSE of 0.034 and 0.026 for the two polarizations, respectively. The proposed algorithm provides an encouraging approach to obtain accurate daily MLSE dataset for microwave remote sensing. Compared to the SMAP MLSE, the CYGNSS MLSE has a remarkable improvement of 86% in temporal resolution, greatly complementing the existing microwave emissivity datasets.

54 ENVIRONMENTAL SCIENCES↗

Models of bars − II. Exponential profiles

ABSTRACT We present a new model for galactic bars with exponentially falling major axis luminosity profiles and Gaussian cross-sections. This is based on the linear superposition of Gaussian potential–density pairs with an exponential weight function, using an extension of the method originally introduced by Long & Murali. We compute the density, potential, and forces, using Gaussian quadrature. These quantities are given as explicit functions of position. There are three independent scaled bar parameters that can be varied continuously to produce bespoke bars of a given mass and shape. We categorize the effective potential by splitting a reduced parameter space into six regions. Unusually, we find bars with three stable Lagrange points on the major axis are possible. Our model reveals a variety of unexpected orbital structure, including a bifurcating x1 orbit coexisting with a stable x4 orbit. Propeller orbits are found to play a dominant role in the orbital structure, and we find striking similarities between our bar configuration and the model of Kaufmann & Contopoulos. We find a candidate orbital family, sired from the propeller orbits, that may be responsible for the observed high-velocity peaks in the Milky Way’s bar. As a cross-check, we inspect, for the first time, the proper motions of stars in the high-velocity peaks, which also match our suggested orbital family well. This work adds to the increasing body of evidence that real galactic bars may be supported at least partly by propeller orbits rather than solely by elliptical-like orbits of the x1 family.

McGough, D. P.↗

Mitigating the effects of undersampling in weak lensing shear estimation with metacalibration

ABSTRACT metacalibration is a state-of-the-art technique for measuring weak gravitational lensing shear from well-sampled galaxy images. We investigate the accuracy of shear measured with metacalibration from fitting elliptical Gaussians to undersampled galaxy images. In this case, metacalibration introduces aliasing effects leading to an ensemble multiplicative shear bias about 0.01 for Euclid and even larger for the Roman Space Telescope, well exceeding the missions’ requirements. We find that this aliasing bias can be mitigated by computing shapes from weighted moments with wider Gaussians as weight functions, thereby trading bias for a slight increase in variance of the measurements. We show that this approach is robust to the point-spread function in consideration and meets the stringent requirements of Euclid for galaxies with moderate to high signal-to-noise ratios. We therefore advocate metacalibration as a viable shear measurement option for weak lensing from upcoming space missions.

Kannawadi, Arun↗