Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistent scheme”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Spectral consequences of photoreceptor sampling in the rhesus retina

Optical transforms were used to compute the power spectra of rhesus cones treated as arrays of image sampling points. Spectra were obtained for the central fovea, parafovea, periphery, and far periphery. All were consistent with a novel spatial sampling principle that introduces minimal noise for spatial frequencies below the Nyquist limits implied by local receptor densities, while frequencies above the nominal Nyquist limits are not converted into conspicuous moiree patterns, but instead are scattered into broadband noise. This sampling scheme allows the visual system to escape aliasing distortion despite a large mismatch between retinal image bandwidth and the Nyquist limits implied by extrafoveal cone densities.

Yellott, J. I., Jr.↗
NASA NTRS1995

AC motor controller with 180 degree conductive switches

An ac motor controller is operated by a modified time-switching scheme where the switches of the inverter are on for electrical-phase-and-rotation intervals of 180.degree. as opposed to the conventional 120.degree.. The motor is provided with three-phase drive windings, a power inverter for power supplied from a dc power source consisting of six switches, and a motor controller which controls the current controlled switches in voltage-fed mode. During full power, each switch is gated continuously for three successive intervals of 60.degree. and modulated for only one of said intervals. Thus, during each 60.degree. interval, the two switches with like signs are on continuously and the switch with the opposite sign is modulated.

Oximberg, Carol A.↗
NASA NTRS2013

A Hybrid Fiber/Solid-State Regenerative Amplifier with Tunable Pulse Widths for Satellite Laser Ranging

A fiber/solid-state hybrid seeded regenerative amplifier, capable of achieving high output energy with tunable pulse widths, has been developed for satellite laser ranging applications. The regenerative amplifier cavity uses a pair of Nd:YAG zigzag slabs oriented orthogonally to one another in order to make thermal lensing effects symmetrical and simplify optical correction schemes. The seed laser used is a fiber-coupled 1,064-nm narrowband (<0.02 nm) diode laser that is discretely driven in a new short-pulsed mode, enabling continuously tunable seed pulse widths in the 0.2-to-0.4-ns range. The amplifier gain unit consists of a pair of Brewster-cut 6-bounce zigzag Nd:YAG laser slabs, oriented 90deg relative to each other in the amplifier head. This arrangement creates a net-symmetrical thermal lens effect (an opposing singleaxis effect in each slab), and makes thermo-optical corrections simple by optimizing the curvature of the nearest cavity mirror. Each slab is pumped by a single 120-W, pulsed 808-nm laser diode array. In this configuration, the average pump beam distribution in the slabs had a 1-D Gaussian shape, which matches the estimated cavity mode size. A half-wave plate between the slabs reduces losses from Fresnel reflections due to the orthogonal slabs Brewster-cut end faces. Successful "temporal" seeding of the regenerative amplifier cavity results in a cavity Q-switch pulse envelope segmenting into shorter pulses, each having the width of the input seed, and having a uniform temporal separation corresponding to the cavity round-trip time of approx. =10 ns. The pulse energy is allowed to build on successive passes in the regenerative amplifier cavity until a maximum is reached, (when cavity gains and losses are equal), after which the pulse is electro- optically switched out on the next round trip The overall gain of the amplifier is approx. =82 dB (or a factor of 1.26 million). After directing the amplified output through a LBO frequency doubling crystal, approx. = 2.1 W of 532-nm output (>1 mJ) was measured. This corresponds to a nonlinear conversion efficiency of >60%. Furthermore, by pulse pumping this system, a single pulse per laser shot can be created for the SLR (satellite laser ranging) measurement, and this can be ejected into the instrument. This is operated at the precise frequency needed by the measurement, as opposed to commercial short-pulsed, mode-locked systems that need to operate in a continuous fashion, or CW (continuous wave), and create pulses at many MHz. Therefore, this design does not need to throw away or dump 99% of the laser energy to produce what is required; this system can be far smaller, more efficient, cheaper, and readily deployed in the field when packaged efficiently. Finally, by producing custom diode seed pulses electronically, two major advantages over commercial systems are realized: First, this pulse shape is customizable and not affected by the cavity length or gain of the amplifier cavity, and second, it can produce adjustable (selectable) pulse widths by simply adding multiple seed diodes and coupling each into commercial, low-cost fiber-optic combiners.

Coyle, Barry↗
DOE OSTI2025

Robust Implicit Adaptive Low Rank Time-Stepping Methods for Matrix Differential Equations

In this work, we develop implicit rank-adaptive schemes for time-dependent matrix differential equations. The dynamic low rank approximation (DLRA) is a well-known technique to capture the dynamic low rank structure based on Dirac–Frenkel time-dependent variational principle. In recent years, it has attracted a lot of attention due to its wide applicability. Our schemes are inspired by the three-step procedure used in the rank adaptive version of the unconventional robust integrator (the so called BUG integrator) (Ceruti et al. in BIT Numer Math 62(4):1149–1174, 2022) for DLRA. First, a prediction (basis update) step is made computing the approximate column and row spaces at the next time level. Second, a Galerkin evolution step is invoked using an implicit solves for the small core matrix. Finally, a truncation is made according to a prescribed error threshold. Since the DLRA is evolving the differential equation projected on to the tangent space of the low rank manifold, the error estimate of the BUG integrator contains the tangent projection (modeling) error which cannot be easily controlled by mesh refinement. This can cause convergence issue for equations with cross terms. To address this issue, we propose a simple modification, consisting of merging the row and column spaces from the explicit step truncation method together with the BUG spaces in the prediction step. In addition, we propose an adaptive strategy where the BUG spaces are only computed if the residual for the solution obtained from the prediction space by explicit step truncation method, is too large. Here, we prove stability and estimate the local truncation error of the schemes under assumptions. We benchmark the schemes in several tests, such as anisotropic diffusion, solid body rotation and the combination of the two, to show robust convergence properties.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The Completed SDSS-IV Extended Baryon Oscillation Spectroscopic Survey: N-body Mock Challenge for Galaxy Clustering Measurements

We develop a series of N-body data challenges, functional to the final analysis of the extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 16 (DR16) galaxy sample. The challenges are primarily based on high-fidelity catalogues constructed from the Outer Rim simulation - a large box size realization (3h(-1) Gpc) characterized by an unprecedented combination of volume and mass resolution, down to 1.85 x 10(9) h(-1)M(circle dot). We generate synthetic galaxy mocks by populating Outer Rim haloes with a variety of halo occupation distribution (HOD) schemes of increasing complexity, spanning different redshift intervals. We then assess the performance of three complementary redshift space distortion (RSD) models in configuration and Fourier space, adopted for the analysis of the complete DR16 eBOSS sample of Luminous Red Galaxies (LRG5). We find all the methods mutually consistent, with comparable systematic errors on the Alcock-Paczynski parameters and the growth of structure, and robust to different HOD prescriptions - thus validating the robustness of the models and the pipelines used for the baryon acoustic oscillation (BAO) and full shape clustering analysis. In particular, all the techniques are able to recover and alpha(11) to within 0.9 per cent, and f sigma(8) to within 1.5 per cent. As a by-product of our work, we are also able to gain interesting insights on the galaxy-halo connection. Our study is relevant for the final eBOSS DR16 'consensus cosmology', as the systematic error budget is informed by testing the results of analyses against these high-resolution mocks. In addition, it is also useful for future large-volume surveys, since similar mock-making techniques and systematic corrections can be readily extended to model for instance the Dark Energy Spectroscopic Instrument (DESI) galaxy sample.

cosmology: theory, large-scale structure of Univer↗
NASA NTRS1991

Hybrid solution for high-speed target acquisition and identification systems

A typical hierarchy for a general object recognition problem consists of object detection, classification and identification. This paper establishes necessary building blocks required for high-speed object recognition applications. An architecture that combines digital and optical processing, exploiting current image processing techniques for detection and classification, and optical processing hardware is described. An optical processing scheme is suggested for the identification aspect. Numerical results of each proposed concept are presented.

Udomkesmalee, Suraphol↗
NASA NTRS1995

Sensitivity analysis of frictional contact response of axisymmetric composite structures

A computational procedure is presented for evaluating the sensitivity coefficients of the static frictional contact response of axisymmetric composite structures. The structures are assumed to consist of an arbitrary number of perfectly bonded homogeneous anisotropic layers. The material of each layer is assumed to be hyperelastic, and the effect of geometric nonlinearity is included. The sensitivity coefficients measure the sensitivity of the response variations in different material, lamination and geometric parameters of the structure. A displacement finite element model is used for the discretization. The normal contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of nodal displacements, and Lagrange multipliers associated with the contact conditions. The Lagrange multipliers are allowed to be discontinuous at interelement boundaries. Tangential contact conditions are incorporated by using a penalty method in conjunction with the classical Coulomb's friction model. The Newton-Raphson iterative scheme is used for the solution of the resulting nonlinear algebraic equations, and for the determination of the contact region, contact conditions (sliding or sticking), and the contact pressures. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented for the frictional contact of a composite spherical cap pressed against a rigid plate.

Karaoglan, I.↗
NASA NTRS1992

A numerical study of three-dimensional surface tension driven convection with fre surface deformation

The steady three-dimensional thermocapillary motion with a deformable free surface is studied numerically in both normal and zero gravity environments. Flow configurations consist of a square cavity heated from the side. In the analysis, the free surface is allowed to deform and the grid distribution is adapted to the surface deformation. The divergence-free condition is satisfied by using a dual time-stepping approach in the numerical scheme. Convective flux derivatives are evaluated using a third-order accurate upwind-biased flux-split differencing technique. The numerical solutions at the midplane of the square cavity are compared with the results from two-dimensional calculations. In addition, numerial results for cases under zero and normal gravity conditions are compared. Significantly different flow structures and surface deformation have been observed. The comparison of calculated results will be compared with experimental data in the updated version of this paper.

Hsieh, Kwang-Chung↗
NASA NTRS1995

Calculations of Diffuser Flows with an Anisotropic K-Epsilon Model

A newly developed anisotropic K-epsilon model is applied to calculate three axisymmetric diffuser flows with or without separation. The new model uses a quadratic stress-strain relation and satisfies the realizability conditions, i.e., it ensures both the positivity of the turbulent normal stresses and the Schwarz' inequality between any fluctuating velocities. Calculations are carried out with a finite-element method. A second-order accurate, bounded convection scheme and sufficiently fine grids are used to ensure numerical credibility of the solutions. The standard K-epsilon model is also used in order to highlight the performance of the new model. Comparison with the experimental data shows that the anisotropic K-epsilon model performs consistently better than does the standard K-epsilon model in all of the three test cases.

Zhu, J.↗
NASA NTRS2001

A Coupled Plasma and Sheath Model for High Density Reactors

We present a coupled plasma and collisionless; sheath model for the simulation of high density plasma processing reactors. Due to inefficiencies in numerical schemes and the resulting computational burden, a coupled multidimensional plasma and sheath simulation has not been possible model for gas mixtures and high density reactors of practical interest. In this work we demonstrate that with a fully implicit algorithm and a refined computational mesh, a self-consistent plasma and sheath simulation is feasible. We discuss the details of the model equations, the importance of ion inertia, and the resulting sheath profiles for argon and chlorine plasmas. We find that at low operating pressures (10-30 mTorr), the charge separation occurs only within a 0.5 mm layer near the surface in a 300 mm inductively coupled plasma etch reactor. A unified model eliminates the use of off-line or loosely coupled sheath models with simplifying assumptions which generally lead to uncertainties in ion flux and sheath electrical properties.

Deepak, Bose↗
NASA NTRS2001

Formation of Electrostatic Potential Drops in the Auroral Zone

In order to examine the self-consistent formation of large-scale quasi-static parallel electric fields in the auroral zone on a micro/meso scale, a particle in cell simulation has been developed. The code resolves electron Debye length scales so that electron micro-processes are included and a variable grid scheme is used such that the overall length scale of the simulation is of the order of an Earth radii along the magnetic field. The simulation is electrostatic and includes the magnetic mirror force, as well as two types of plasmas, a cold dense ionospheric plasma and a warm tenuous magnetospheric plasma. In order to study the formation of parallel electric fields in the auroral zone, different magnetospheric ion and electron inflow boundary conditions are used to drive the system. It has been found that for conditions in the primary (upward) current region an upward directed quasi-static electric field can form across the system due to magnetic mirroring of the magnetospheric ions and electrons at different altitudes. For conditions in the return (downward) current region it is shown that a quasi-static parallel electric field in the opposite sense of that in the primary current region is formed, i.e., the parallel electric field is directed earthward. The conditions for how these different electric fields can be formed are discussed using satellite observations and numerical simulations.

Schriver, D.↗
NASA NTRS2017

Geometrically Flexible and Efficient Flow Analysis of High Speed Vehicles Via Domain Decomposition, Part 1: Unstructured-Grid Solver for High Speed Flows

The ability to solve the equations governing the hypersonic turbulent flow of a real gas on unstructured grids using a spatially-elliptic, 2nd-order accurate, cell-centered, finite-volume method has been recently implemented in the VULCAN-CFD code. This paper describes the key numerical methods and techniques that were found to be required to robustly obtain accurate solutions to hypersonic flows on non-hex-dominant unstructured grids. The methods and techniques described include: an augmented stencil, weighted linear least squares, cell-average gradient method, a robust multidimensional cell-average gradient-limiter process that is consistent with the augmented stencil of the cell-average gradient method and a cell-face gradient method that contains a cell skewness sensitive damping term derived using hyperbolic diffusion based concepts. A data-parallel matrix-based symmetric Gauss-Seidel point-implicit scheme, used to solve the governing equations, is described and shown to be more robust and efficient than a matrix-free alternative. In addition, a y+ adaptive turbulent wall boundary condition methodology is presented. This boundary condition methodology is deigned to automatically switch between a solve-to-the-wall and a wall-matching-function boundary condition based on the local y+ of the 1st cell center off the wall. The aforementioned methods and techniques are then applied to a series of hypersonic and supersonic turbulent flat plate unit tests to examine the efficiency, robustness and convergence behavior of the implicit scheme and to determine the ability of the solve-to-the-wall and y+ adaptive turbulent wall boundary conditions to reproduce the turbulent law-of-the-wall. Finally, the thermally perfect, chemically frozen, Mach 7.8 turbulent flow of air through a scramjet flow-path is computed and compared with experimental data to demonstrate the robustness, accuracy and convergence behavior of the unstructured-grid solver for a realistic 3-D geometry on a non-hex-dominant grid.

White, Jeffery A.↗
NASA NTRS2020

Reduction of structural impacts and distinction of photosynthetic pathways in a global estimation of GPP from space-borne solar-induced chlorophyll fluorescence

Quantifying global photosynthesis remains a challenge due to a lack of accurate remote sensing proxies. Solar-induced chlorophyll fluorescence (SIF) has been shown to be a good indicator of photosynthetic activity across various spatial scales. However, a global and spatially challenging estimate of terrestrial gross primary production (GPP) based on satellite SIF remains unresolved due to the confounding effects of species-specific physical and physiological traits and external factors, such as canopy structure or photosynthetic pathway (C3orC4). Here we analyze an ensemble of far-red SIF data from OCO-2 satellite and ground observations at multiple sites, using the spectral invariant theory to reduce the effects of canopy structure and to retrieve a structure-corrected total canopy SIF emission (SIF total). We find that the relationships between observed canopy-leavingSIF and ecosystem GPP vary significantly among biomes. In contrast, the relationships between SIF total and GPP converge around two unique models, one for C3and one for C4plants. We show that the two single empirical models can be used to globally scale satellite SIF observations to terrestrial GPP. We obtain an independent estimate of global terrestrial GPP of 129.56 ± 6.54 PgC/year for the 2015–2017 period, which is consistent with the state-of-the-art data- and process-oriented models. The new GPP product shows improved sensitivity to previously undetected ‘hotspots’ of productivity, being able to resolve the double-peak in GPP due to rotational cropping systems. We suggest that the direct scheme to estimate GPP presented here, which is based on satelliteSIF, may open up new possibilities to resolve the dynamics of global terrestrial GPP across space and time.

scale satellite SIF observations to terrestrial GP↗
DOE OSTI2022

Size-Dependent Characteristics of Surface-Rooted Three-Dimensional Convective Objects in Continental Shallow Cumulus Simulations

A segmentation algorithm is applied to high resolution simulations of shallow continental convection to identify individual convective 3D objects within the convective boundary layer, in order to investigate which properties of the objects vary with the object width. The study analyses the geometry of the objects, along with their profiles of vertical velocity and total water, to assess various assumptions often used in spectral mass-flux convection schemes. The methodology of this paper is unique in that we use (a) a novel application of the watershed algorithm to detect individual objects in the constantly evolving continental boundary layer efficiently, and (b) an unprecedentedly large number of simulations being analyzed. In total, 26 days of LASSO simulations at the Atmospheric Radiation Measurement Southern Great Plains site are analyzed, yielding roughly one million objects. Plume-like surface-rooted objects are found to be omnipresent, the vertical extent of which is strongly dependent on the object width. The vertical velocity and moisture anomaly profiles of the widest objects are roughly consistent with the classic buoyancy-driven rising plume model. The kinematic and thermodynamic properties of the objects vary with object width. This width dependence is strongest above cloud base, but much weaker below. Finally the impact of neglecting the contribution of covariances to the vertical moisture flux, which is commonly used in mass-flux parameterizations, is investigated. The average effect of neglecting covariances increases linearly with object width, leading to a 20% flux underestimation for 2 km wide objects. Implications of the results for spectral convection scheme development are briefly discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Advanced Laboratory and Field Arrays: Technical Report for Performance Enhancement for Marine Energy Converter (MEC) Arrays (Task 5)

Task 5 consisted of two parts, wave energy converter (WEC) arrays and current energy converter (CEC) arrays. The two parts of the project were executed independently, and this report covers only the work on WECs. In this part of the task, the team characterized the behavior of and developed control schemes for wave energy converter (WEC) arrays that improved performance (i.e., maximized aggregate power generation and reduced the levelized cost of energy) over baseline, non-coordinated control approaches. Specifically, for WEC arrays: the team developed optimal layouts of arrays that considered the effect of WEC placement within an array on coordinated array control with real-time estimation schemes. The numerical codes developed through this task were made available to users for array design. These tools utilized existing commercial software as well as software developed within the project team.

13 HYDRO ENERGY↗
NASA NTRS1973

An unsupervised classification technique for multispectral remote sensing data.

Description of a two-part clustering technique consisting of (a) a sequential statistical clustering, which is essentially a sequential variance analysis, and (b) a generalized K-means clustering. In this composite clustering technique, the output of (a) is a set of initial clusters which are input to (b) for further improvement by an iterative scheme. This unsupervised composite technique was employed for automatic classification of two sets of remote multispectral earth resource observations. The classification accuracy by the unsupervised technique is found to be comparable to that by traditional supervised maximum-likelihood classification techniques.

Su, M. Y.↗
NASA NTRS1982

Mature thunderstorm cloud top structure - Three-dimensional numerical simulation versus satellite observations

Preliminary results of four runs with a three-dimensional model of the effects of vertical wind shear on cloud top height/temperature structure and the internal properties of isolate midlatitude thunderstorms are reported. The model is being developed as an aid to analyses of GEO remote sensing satellite data. The grid is a 27 x 27 x 20 mesh with 2 km horizontal resolution and 0.9 vertical resolution. The total grid is 54 km on a side and 18 km deep. A second-order Crowley scheme for advecting momentum is extended with a third-order correction for spatial truncation error, and the earth-relative horizontal surface wind components are decreased to 50 percent of their values at 0.45 km. A temperature increase with height is included, together with an initial impulse consisting of a nonrotating cylindrical weak buoyant updraft 10 km in radius. The results of the runs are discussed in terms of the time variation of the vertical velocity extrema, the effects of strong and weak shear on a storm, the cloud top height, the Lagrangian dynamics of a thermal couplet, and data from a real storm.

Schlesinger, R. E.↗
NASA NTRS1984

A combined direct/inverse three-dimensional transonic wing design method for vector computers

A three-dimensional transonic-wing design algorithm for vector computers is developed, and the results of sample computations are presented graphically. The method incorporates the direct/inverse scheme of Carlson (1975), a Cartesian grid system with boundary conditions applied at a mean plane, and a potential-flow solver based on the conservative form of the full potential equation and using the ZEBRA II vectorizable solution algorithm of South et al. (1980). The accuracy and consistency of the method with regard to direct and inverse analysis and trailing-edge closure are verified in the test computations.

Weed, R. A.↗