Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear Parameter Varying”

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

Progressive Aerodynamic Model Identification From Dynamic Water Tunnel Test of the F-16XL Aircraft

Development of a general aerodynamic model that is adequate for predicting the forces and moments in the nonlinear and unsteady portions of the flight envelope has not been accomplished to a satisfactory degree. Predicting aerodynamic response during arbitrary motion of an aircraft over the complete flight envelope requires further development of the mathematical model and the associated methods for ground-based testing in order to allow identification of the model. In this study, a general nonlinear unsteady aerodynamic model is presented, followed by a summary of a linear modeling methodology that includes test and identification methods, and then a progressive series of steps suggesting a roadmap to develop a general nonlinear methodology that defines modeling, testing, and identification methods. Initial steps of the general methodology were applied to static and oscillatory test data to identify rolling-moment coefficient. Static measurements uncovered complicated dependencies of the aerodynamic coefficient on angle of attack and sideslip in the stall region making it difficult to find a simple analytical expression for the measurement data. In order to assess the effect of sideslip on the damping and unsteady terms, oscillatory tests in roll were conducted at different values of an initial offset in sideslip. Candidate runs for analyses were selected where higher order harmonics were required for the model and where in-phase and out-of-phase components varied with frequency. From these results it was found that only data in the angle-of-attack range of 35 degrees to 37.5 degrees met these requirements. From the limited results it was observed that the identified models fit the data well and both the damping-in-roll and the unsteady term gain are decreasing with increasing sideslip and motion amplitude. Limited similarity between parameter values in the nonlinear model and the linear model suggest that identifiability of parameters in both terms may be a problem. However, the proposed methodology can still be used with careful experiment design and carefully selected values of angle of attack, sideslip, amplitude, and frequency of the oscillatory data.

Murphy, Patrick C.↗

Constrained model calibration of grain structure dependent spall dynamics in shock-loaded tantalum

Here, we perform a gas gun experiment by shock loading tantalum samples of varying grain structures to assess the suitability of a numerical model for simulating spall behavior. The observed differences in spall strength, as well spallation and re-compression history, are not captured in uncalibrated hydrodynamic simulations. An optimization is performed on the Johnson spall model to determine the best parameters that fit the observed trends. Linear stability analysis is employed to motivate bounds on those parameters. Herein, optimized simulations agree well with the experimental results, reproducing pullback depth and recompression timescales across the different samples tested. Further, the observed pullback time of the single crystal sample was found to imply, via the stability analysis, a percolation threshold in good agreement with the theoretical value for a body centered cubic lattice. Therefore, the combined linear stability and percolation analysis shows promise and may be applied to other materials with diverse microstructures. Collectively, the findings demonstrate that the model is suitable for reproducing spall-induced free surface behavior across various microstructures, but also points to caution in using model coefficients for uncalibrated microstructures.

58 GEOSCIENCES↗

Closed Loop Geothermal Working Group: GeoCLUSTER App, Subsurface Simulation Results, and Publications

To better understand the heat production, electricity generation performance, and economic viability of closed-loop geothermal systems in hot-dry rock, the Closed-Loop Geothermal Working Group -- a consortium of several national labs and academic institutions has tabulated time-dependent numerical solutions and levelized cost results of two popular closed-loop heat exchanger designs (u-tube and co-axial). The heat exchanger designs were evaluated for two working fluids (water and supercritical CO2) while varying seven continuous independent parameters of interest (mass flow rate, vertical depth, horizontal extent, borehole diameter, formation gradient, formation conductivity, and injection temperature). The corresponding numerical solutions (approximately 1.2 million per heat exchanger design) are stored as multi-dimensional HDF5 datasets and can be queried at off-grid points using multi-dimensional linear interpolation. A Python script was developed to query this database and estimate time-dependent electricity generation using an organic Rankine cycle (for water) or direct turbine expansion cycle (for CO2) and perform a cost assessment. This document aims to give an overview of the HDF5 database file and highlights how to read, visualize, and query quantities of interest (e.g., levelized cost of electricity, levelized cost of heat) using the accompanying Python scripts. Details regarding the capital, operation, and maintenance and levelized cost calculation using the techno-economic analysis script are provided. This data submission will contain results from the Closed Loop Geothermal Working Group study that are within the public domain, including publications, simulation results, databases, and computer codes. GeoCLUSTER is a Python-based web application created using Dash, an open-source framework built on top of Flask that streamlines the building of data dashboards. GeoCLUSTER provides users with a collection of interactive methods for streamlining the exploration and visualization of an HDF5 dataset. The GeoCluster app and database are contained in the compressed file geocluster_vx.zip, where the "x" refers to the version number. For example, geocluster_v1.zip is Version 1 of the app. This zip file also contains installation instructions. **To use the GeoCLUSTER app in the cloud, click the link to "GeoCLUSTER on AWS" in the Resources section below. To use the GeoCLUSTER app locally, download the geocluster_vx.zip to your computer and uncompress this file. When uncompressed this file comprises two directories and the geocluster_installation.pdf file. The geo-data app contains the HDF5 database in condensed format, and the GeoCLUSTER directory contains the GeoCLUSTER app in the subdirectory dash_app, as app.py. The geocluster_installation.pdf file provides instructions on installing Python, the needed Python modules, and then executing the app.

15 GEOTHERMAL ENERGY↗

Experimental results for the transpired turbulent boundary layer in an adverse pressure gradient

The fluid mechanics of transpired incompressible turbulent boundary layers under zero and adverse pressure gradient conditions is investigated using an open-ended wind tunnel with a porous floor in the test section and a secondary air system for supply and metering of the transpiration air. All velocity profiles and turbulence profiles are obtained by linearized constant-temperature hot-wire anemometry. The wall shear stress is determined by measuring the shear stress away from the wall and extrapolating to the wall by integrating the boundary layer equations for the shear-stress profile. Equilibrium boundary layers are obtained when the transpiration velocity is varied such that the blowing parameter and the Clauser pressure gradient parameter are held constant. The experimental results obtained are presented in tabular and graphical forms.

Andersen, P. S.↗

Effect of Operating Frequency on PDE Driven Ejector Thrust Performance

Results of an on-going study of pulse detonation engine driven ejectors are presented and discussed. The experiments were conducted using a pulse detonation engine (PDE) designed to operate at frequencies up to 50 Hz. The PDE used in these experiments utilizes an equi-molar mixture of oxygen and nitrogen as the oxidizer, and ethylene (C2H4) as the fuel, with the propellant mixture having an equivalence ratio of one. A line of sight laser absorption technique was used to determine the time needed for proper filling of the tube. Thrust measurements were made using an integrated spring damper system coupled with a linear variable displacement transducer. The baseline thrust of the PDE was first measured at each desired frequency and agrees with experimental and modeling results found in the literature. Thrust augmentation measurements were then made for constant diameter ejectors. The ejectors had varying lengths, and two different inlet geometries were tested for each ejector configuration. The parameter space for the study included PDE operation frequency, ejector length, overlap distance and the radius of curvature for the ejector inlets. For the studied experimental matrix, the results showed a maximum thrust augmentation of 106% at an operational frequency of 30 Hz.

Santoro, Robert J.↗

Utilizing the Vertical Variability of Precipitation to Improve Radar QPE

Characteristics of the melting layer and raindrop size distribution can be exploited to further improve radar quantitative precipitation estimation (QPE). Using dual-polarimetric radar and disdrometers, we found that the characteristic size of raindrops reaching the ground in stratiform precipitation often varies linearly with the depth of the melting layer. As a result, a radar rainfall estimator was formulated using D(sub m) that can be employed by polarimetric as well as dual-frequency radars (e.g., space-based radars such as the GPM DPR), to lower the bias and uncertainty of conventional single radar parameter rainfall estimates by as much as 20%. Polarimetric radar also suffers from issues associated with sampling the vertical distribution of precipitation. Hence, we characterized the vertical profile of polarimetric parameters (VP3)-a radar manifestation of the evolving size and shape of hydrometeors as they fall to the ground-on dual-polarimetric rainfall estimation. The VP3 revealed that the profile of ZDR in stratiform rainfall can bias dual-polarimetric rainfall estimators by as much as 50%, even after correction for the vertical profile of reflectivity (VPR). The VP3 correction technique that we developed can improve operational dual-polarimetric rainfall estimates by 13% beyond that offered by a VPR correction alone.

Precipitation↗

Electrical Properties of Unidirectional Carbon Fiber Reinforced Epoxy-Graphene Composites

This paper aims at establishing connections between electrical resistivity and stresses in carbon fiber reinforced epoxy-graphene composites. To avoid the interconnection of multiple parameters, we examined composite specimens consisting of a single strand or tow (12,000 carbon fiber filaments) embedded in an epoxy matrix containing up to 1.0 wt% of graphene (TG679) particles. The specimens show resistance increase with the tensile stress. The characteristic of the dependence varies with the content of graphene particles, tending to be linear as the percentage of graphene content increases. We did not observe any significant changes in the ultimate strength of the composite with an increase of the graphene particles concentration.

Daniel Wentzel↗

Satellite Constellation Orbit Design to Enable a Space-Based Radio Interferometer

Two different design methods for a networked constellation of N small satellites are presented. In the first method, the (linear) Clohessy-Wilthsire equations are used as an initial design tool, followed by conversion to a two-body model. Discrepancies between the linear and nonlinear solutions are minimized in the conversion process. The second method utilizes invariant manifold theory, by perturbing a reference trajectory in different directions along the center eigenvectors. Both methods require 5N parameters to fully define a constellation. In a relative, rotating frame of the reference path, the spacecraft appear as periodic ellipses of varying sizes. Deployment, reconfiguration using propulsive maneuvers, and station keeping costs for an example mission scenario are addressed.

Chien, Steve A.↗

SLR Station Recovery, Center of Frame Motion, and Time Varying Gravity

Weekly station position estimates, beginning with 1993, are derived from the ITRF2008-based SLR processing of up to four satellites: Lageos 1, Lageos2, Starlette, and Stella. Helmert parameters obtained from c omparison of weekly SLR station positions and the a-priori SLRF2008 station complement are evaluated for geocenter motion and scale. Two me thods for modeling time varying gravity are employed in the SLR satel lite POD processing, with GGM03S serving as the static gravity field. Both methods forward model atmosphere gravity derived from 6-hour ECM WF pressure data. The standard approach applies an annual 20x20 field estimated from 4 years of GRACE data, and the IERS2003 recommended linear rates for C20, C30, C40, C21, and S21. The alternate approach us es a new set of low-order/degree 4x4 coefficients estimated weekly fr om SLR & DORIS processing to 10 satellites from 1993-2012. This exper imental tvg4x4 model has been shown to improve the TOPEX, Jason-1, and Jason-2 altimeter satellite orbits,. In this paper we apply the more detailed time-variable gravity modeling to the SLR satellite POD pro cessing and subsequent reference frame analyses. For this study we will evaluate the orbit differences (periodic and secular) for the satel lites concerned, characterize the impact on the station coordinate solutions, and the impact on reference frame parameters (geocenter and s cale).

Zelensky, Nikita P.↗

Gortler vortices and transition in wall boundary layers of two Mach 5 nozzles

The onset of transition in the wall boundary layers of two axisymmetric Mach 5 wind-tunnel nozzles has been measured under conditions of extremely low incident disturbance levels. The range of test unit Reynolds numbers, based on conditions at the nozzle exit, was from 6 x 10 to the 6th power m to 2.5 x 10 to the 7th power m. When the nozzle walls were maintained in a polished and clean condition, transition moved gradually upstream as the test Reynolds number was increased. When transition occurred in the supersonic concave wall region, the values of the local Gortler parameter at transition varied from about 5 to 6, whereas the momentum thickness Reynolds number varied from about 750 to 1050. Oil flow patterns obtained near the exit of the nozzles indicated that Gortler vortices were always present when the wall boundary layers were laminar. Calculations for the growth of Gortler vortices based on new results from linear theory for supersonic flat-plate profiles gave amplification ratios to transition from e to the 4th power to e to the 15th power. Possible reasons for this wide range in amplification ratios are discussed, but no definite conclusions are yet possible regarding the values of n in a simple e to the nth power type theory for the assumed linear amplification of Gortler vortices to transition in supersonic nozzles.

Beckwith, I. E.↗

Tabulated Database of Closed-Loop Geothermal Systems Performance for Cloud-Based Technical and Economic Modeling of Heat Production and Electricity Generation: Preprint

To better understand the heat production, electricity generation performance and economic viability of closed loop geothermal systems in hot-dry-rock, the Closed Loop Geothermal Group, a consortium of several national labs and academic institutions has tabulated time-dependent numerical solutions and levelized cost results of two popular closed loop heat exchanger designs (u-tube and co-axial). The heat exchanger designs were evaluated for two working fluids (water and super-critical CO2) while varying seven continuous independent parameters of interest (i.e., mass flow rate, vertical depth, horizontal extent, borehole diameter, formation gradient, formation conductivity, and injection temperature). The corresponding numerical solutions (approximately 1.2 million per heat exchanger design) are stored as multi-dimensional HDF5 datasets and can be queried at off-grid points using multi-dimensional linear interpolation. A Python script was developed to query this database and estimate time-dependent electricity generation using an Organic Rankine cycle (for water) or direct turbine expansion cycle (for CO2) and perform a cost assessment. This document aims to give an overview of the HDF5 database file and highlights how to read, visualize, and query quantities of interest (e.g., levelized cost of electricity, levelized cost of heat) using the accompanying python scripts. Details regarding the capital, operation, and maintenance and levelized cost calculation using the TEA (techno-economic analysis) script are provided.

co-axial↗

Problems in fluid dynamics

A scheme was developed for the parametric differentiation and integration of gas dynamics equations. A numerical integration of the gas dynamics equations is necessarily performed for a specific set of parameter values. The linear variational equations are obtained by differentiating the exact equations with respect to each of the relevant parameters. The resulting matrix of flow quantities is referred to as the Jacobi matrix. The subsequent procedure is then straightforward. The method was tested for two dimensional supersonic flow past an airfoil, with airfoil thickness, camber, and angle of attack varied. This approach has great potential value for rapidly assessing the effect of design changes. The other focus of the work was on problems in fluid stability, bifurcations, and turbulence.

Sirovich, Lawrence↗

The Q 10 of in situ microbial soil respiration varies with mean annual temperature, precipitation, pH, and plant cover: a meta-analysis and spatial prediction of Q 10

The temperature sensitivity of soil microbial respiration, commonly quantified using the Q 10 coefficient, is a key parameter in carbon cycle models. Uncovering how environmental factors affect in situ Q 10 values can therefore provide critical insight into potential shifts in global carbon stocks under climate change. We collected data from previously published field experiments that measured soil microbial respiration across a range of temperatures. We hypothesized that the Q 10 coefficient of in situ soil microbial respiration would vary based on environmental factors including mean annual temperature (MAT), mean annual precipitation (MAP), plant cover type, pH, soil C:N, and latitude. Linear regression revealed that Q 10 correlates negatively with MAT and MAP and positively with pH and absolute latitude. Additionally, average Q 10 varied significantly across different plant cover types; it was highest in mountain grasslands and lowest in tropical moist forests. Variation in microbial Q 10 across environmental factors may arise from underlying mechanisms such as enzyme kinetics, substrate availability and complexity, and microbial adaptation. To capture patterns in Q 10 more comprehensively, we developed a multiple linear regression model of Q 10 based on the most individually significant environmental drivers and applied it to public datasets to generate a global map of predicted Q 10 . Q 10 was higher in high-latitude and high-altitude regions, where large permafrost carbon stores are vulnerable to thawing and decomposition. We also compared fits between the Q 10 equation and a model produced from macromolecular rate theory (MMRT). We found that the MMRT model had the superior fit and may be better suited to model temperature sensitivity of complex biological reactions. Overall, our results emphasize that relationships between microbial Q 10 and environmental variables should be accounted for in climate models. Incorporating these variations in the Q 10 parameter, rather than using a fixed value, will help predict whether CO 2 emissions will be buffered or exacerbated by soil microbial respiration under climate change.

54 ENVIRONMENTAL SCIENCES↗

Flow and Arching of Biomass Particles in Wedge-Shaped Hoppers

Poor understanding of the flow behavior of granular biomass material poses great challenges for the bioenergy industry, as the equipment functioning time is significantly reduced by handling issues like screw feeder clogging and hopper arching. In this work, the flow behavior of loblolly pine chips, including the mass flow rate and the critical outlet width, in a wedge-shaped hopper is investigated by combining physical experiments and numerical simulations. Comprehensive characterization of the flow response affected by the two material attributes (initial packing, particle density) and the three operational parameters (hopper outlet width, hopper inclination, and surcharge) is conducted. The results show that the hopper outlet width linearly controls the mass flow rate while the hopper inclination angle controls the critical outlet size. The packing determines whether the flow is smooth or surging, and the surcharge-induced compaction creates flow impedance. The magnitude of these influences varies from a slender hopper with a low inclination angle to a flat-bottom silo. These findings provide guidance for hopper operation in the material handling industry and shed light on the construction of a novel design method for material handling equipment in biorefineries.

09 BIOMASS FUELS↗

Multi-objective optimization of peel and shear strengths in ultrasonic metal welding using machine learning-based response surface methodology

Ultrasonic metal welding (UMW) is a solid-state joining technique with varied industrial applications. Despite of its numerous advantages, UMW has a relative narrow operating window and is sensitive to variations in process conditions. As such, it is imperative to quantitatively characterize the influence of welding parameters on the resulting joint quality. The quantification model can be subsequently used to optimize the parameters. Conventional response surface methodology (RSM) usually employs linear or polynomial models, which may not be able to capture the intricate, nonlinear input-output relationships in UMW. Furthermore, some UMW applications call for simultaneous optimization of multiple quality indices such as peel strength, shear strength, electrical conductivity, and thermal conductivity. To address these challenges, this paper develops a machine learning (ML)-based RSM to model the input-output relationships in UMW and jointly optimize two quality indices, namely, peel and shear strengths. The performance of various ML methods including spline regression, Gaussian process regression (GPR), support vector regression (SVR), and conventional polynomial regression models with different orders is compared. A case study using experimental data shows that GPR with radial basis function (RBF) kernel and SVR with RBF kernel achieve the best prediction accuracy. The obtained response surface models are then used to optimize a compound joint strength indicator that is defined as the average of normalized shear and peel strengths. In addition, the case study reveals different patterns in the response surfaces of shear and peel strengths, which has not been systematically studied in the literature. While developed for the UMW application, the method can be extended to other manufacturing processes.

42 ENGINEERING↗

Inductive Non-Contact Position Sensor

Optical hardware has been developed to measure the depth of defects in the Space Shuttle Orbiter's windows. In this hardware, a mirror is translated such that its position corresponds to the defect's depth, so the depth measurement problem is transferred to a mirror-position measurement problem. This is preferable because the mirror is internal to the optical system and thus accessible. Based on requirements supplied by the window inspectors, the depth of the defects needs to be measured over a range of 200 microns with a resolution of about 100 nm and an accuracy of about 400 nm. These same requirements then apply to measuring the position of the mirror, and in addition, since this is a scanning system, a response time of about 10 ms is needed. A market search was conducted and no sensor that met these requirements that also fit into the available housing volume (less than one cubic inch) was found, so a novel sensor configuration was constructed to meet the requirements. This new sensor generates a nearly linearly varying magnetic field over a small region of space, which can easily be sampled, resulting in a voltage proportional to position. Experiments were done with a range of inductor values, drive voltages, drive frequencies, and inductor shapes. A rough mathematical model was developed for the device that, in most aspects, describes how it operates and what electrical parameters should be chosen for best performance. The final configuration met all the requirements, yielding a small rugged sensor that was easy to use and had nanometer resolution over more than the 200-micron range required. The inductive position sensor is a compact device (potentially as small as 2 cubic centimeters), which offers nanometer-position resolution over a demonstrated range of nearly 1 mm. One of its advantages is the simplicity of its electrical design. Also, the sensor resolution is nearly uniform across its operational range, which is in contrast to eddy current and capacitive sensors whose sensitivity is dependent upon position.

Youngquist, Robert↗

Comparing the influence of cation order and composition in simulated Zn(Sn, Ge)N 2 on structure, elastic moduli, and polarization for solid state lighting

Alloying and site ordering play complementary roles in dictating a material’s properties. However, deconvolving the impacts of these separate phenomena can be challenging. In this work, we simulate structures of Zn(Sn,Ge)N 2 with varied Sn content and site ordering to determine the impacts of order and composition on structural and electronic properties. We assess the formation enthalpy, lattice parameters, elastic constants, spontaneous polarization, and piezoelectric coefficients. In mostly disordered structures (order parameters ranging from 0.2 to 0.4), the formation enthalpy exhibits local extrema as a function of the order parameter, deviating from the more linear trends seen in both fully disordered and fully ordered systems. This anomalous deviation from the otherwise linear trend in formation enthalpy with order manifests in each of the other properties calculated. This range of order parameters of interest may be caused by a transition in the ordering of the quaternary material similar to phase changes seen in ternary compounds but stretched over a region inclduing 20% of the order parameter range. Most parameters calculated are more sensitive to order than to composition in the limited composition range tested; however, the lattice parameter c, piezoelectric coefficient e 33 , and elastic moduli C 12 , C 13 , and C 23 are more sensitive to composition. Of the properties compared, the piezoelectric coefficients are influenced most significantly by changes in both the composition and order parameter. Lattice parameters undergo the smallest changes with order and composition, but these small differences appear to impart large trends in the other properties. Better understanding the effects of disorder and group IV alloying in Zn(Sn,Ge)N 2 allows for more accurate modeling of characteristics of this material system for solid state lighting and other applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dark and bright modes, and their coherent control in dipolar metasurface bilayers

Several plasmonic nanoparticles supporting dipolar resonances can couple to form normal modes. Here, we develop an analytical model to explain the formation of nonradiative “dark” and radiative “bright” modes through radiative coupling in bilayers consisting of dipolar nanoantenna arrays that are separated by a subwavelength distance. We also include near-field contributions in our model and show that the absorption and reflectance spectra obtained from our model agree reasonably well with the respective finite-difference time-domain simulation results for both perfectly aligned and misaligned bilayers. Here, the ability to vary the reflection and absorption spectra of these bilayers by changing the material and geometrical parameters has potential applications in the design of efficient spectral filters. We also show that we can selectively excite these modes by adjusting the phase between two counterpropagating normally incident fields, which has applications in all-optical modulators and switches based on purely linear interferometric effects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗