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 415 records · Page 23

Airborne and spaceborne radar images for geologic and environmental mapping in the Amazon rain forest, Brazil

Spaceborne and airborne radar image of portions of the Middle and Upper Amazon basin in the state of Amazonas and the Territory of Roraima are compared for purposes of geological and environmental mapping. The contrasted illumination geometries and imaging parameters are related to terrain slope and surface roughness characteristics for corresponding areas that were covered by each of the radar imaging systems. Landforms range from deeply dissected mountain and plateau with relief up to 500 m in Roraima, revealing ancient layered rocks through folded residual mountains to deeply beveled pediplain in Amazonas. Geomorphic features provide distinct textural signatures that are characteristic of different rock associations. The principle drainages in the areas covered are the Rio Negro, Rio Branco, and the Rio Japura. Shadowing effects and low radar sensitivity to subtle linear fractures that are aligned parallel or nearly parallel to the direction of radar illumination illustrate the need to obtain multiple coverage with viewing directions about 90 degrees. Perception of standing water and alluvial forest in floodplains varies with incident angle and with season. Multitemporal data sets acquired over periods of years provide an ideal method of monitoring environmental changes.

Ford, John P.↗

Sensitivity of Optimal Midcourse Correction Scheduling for Robust Cislunar Trajectory Design

A new approach to optimal trajectory design is the determination of optimal trajectories that are robust to initial trajectory dispersions, navigation errors, maneuver execution errors, and environment modeling errors. This paper investigates the sensitivity of cislunar robust optimal trajectory design to launch date, duration of navigation measurement passes, and navigation measurement frequency. For a given cislunar trajectory from translunar injection (TLI) to lunar orbit insertion (LOI), the optimal locations of midcourse corrections, also known as trajectory correction maneuvers (TCM) are determined by minimizing the final 3-σ ∆v subject to a final 3-σ position dispersion constraint for a given launch date, specified measurement pass duration prior to each maneuver, and measurement frequency. Optical navigation (OpNav) is assumed, and OpNav field-of-view (FOV) and lighting constraints are employed. These constraints turn out to be important elements of the problem. The sensitivity of the optimal TCM locations are then investigated by varying the launch date, duration of OpNav measurement passes, and the OpNav measurement frequency, and then re-optimizing the locations of the TCMs. Given the problem parameters provided herein, results show that while the optimal TCM locations with respect to TLI vary greatly from one launch date to another, their locations with respect to LOI are nearly invariant over a 2-month launch window. Results also show that in all cases the optimal location of the last TCM is found to be at the point where the OpNav lunar FOV constraint is first violated. For all other TCMs, OpNav measurement pass duration and measurement frequency can have a moderate to large affect on the optimal TCM locations.

Linear Covariance Analysis↗

Physics‐Informed Machine Learning Method for Large‐Scale Data Assimilation Problems

Abstract We develop a physics‐informed machine learning approach for large‐scale data assimilation and parameter estimation and apply it for estimating transmissivity and hydraulic head in the two‐dimensional steady‐state subsurface flow model of the Hanford Site given synthetic measurements of said variables. In our approach, we extend the physics‐informed conditional Karhunen‐Loéve expansion (PICKLE) method to modeling subsurface flow with unknown flux (Neumann) and varying head (time‐dependent Dirichlet) boundary conditions. We demonstrate that the PICKLE method is comparable in accuracy with the standard maximum a posteriori (MAP) method, but is significantly faster than MAP for large‐scale problems. Both methods use a mesh to discretize the computational domain. In MAP, the parameters and states are discretized on the mesh; therefore, the size of the MAP parameter estimation problem directly depends on the mesh size. In PICKLE, the mesh is used to evaluate the residuals of the governing equation, while the parameters and states are approximated by the truncated conditional Karhunen‐Loéve expansions with the number of parameters controlled by the smoothness of the parameter and state fields, and not by the mesh size. For a considered example, we demonstrate that the computational cost of PICKLE increases near linearly (as N 1.15 ) with the number of grid nodes N , while that of MAP increases much faster (as N 3.28 ). We also show that once trained for one set of Dirichlet boundary conditions (i.e., one river stage), the PICKLE method provides accurate estimates of the hydraulic head for any value of the Dirichlet boundary conditions (i.e., for any river stage).

97 MATHEMATICS AND COMPUTING↗

Transonic adaptive flutter suppression using approximate unsteady time domain aerodynamics

A digital adaptive controller is applied to the active flutter suppression problem of a wing under time varying flight conditions in subsonic and transonic flow. Linear quadratic controller gain at each time step is obtained using an iterative Riccati solver. The digital adaptive optimal controller is robust with respect to the unknown external loads. Flutter and divergence instabilities are simultaneously suppressed using a trailing-edge control surface and displacement sensing. A new transonic unsteady aerodynamic approximation methodology is developed which enables one to carry out the rapid calculation required for transonic aeroservoelastic applications. This approximation is based on a combination of unsteady subsonic aerodynamics combined with a transonic correction procedure. Aeroservoelastic transient time response is obtained using Roger's approximation, state transition matrices and an iterative time marching algorithm. The aeroservoelastic system in the time domain is modeled using a deterministic ARMA model together with a parameter estimator. Transonic flutter boundaries of a wing structure are computed, in the time domain, using an estimated aeroelastic system matrix and are in good agreement with experimental data for the low transonic Mach number range.

Pak, Chan-Gi↗

An Experimental Study at a Mach Number of 3 of the Effect of Turbulence Level and Sandpaper Type Roughness on Transition on a Flat Plate

An investigation has been conducted at a Mach number of 3 of the effect of turbulence level and sandpaper-type roughness on transition for a flat plate. The Reynolds number varied from 0.8 x 10(exp 6) to 1.8 x 10(exp 6) per inch; the settling-chamber turbulence level varied from 0.7 percent to 35 percent; and the heat transfer between the plate and the stream was negligible. Transition locations were determined by an optical method. This method was indicative of a permanent change in the boundary-layer density distribution rather than the onset of turbulent bursts. Results showed that, when transition was influenced by roughness, it moved in a way similar to its movement on a smooth plate. That is, it gradually approached the roughness location with either an increase in unit Reynolds number or an increase in turbulence level. For roughness submerged in the linear portion of the boundary-layer velocity profile, the square root of the roughness Reynolds number and the ratio of roughness height to boundary-layer displacement thickness gave similar results as parameters for predicting the effects of roughness. A range of each of these parameters which moved transition less than 10 percent was found and this range was a function of turbulence level.

Jones, Robert A.↗

PID-Regulated Heating System for PIP-II Reference Line

The Proton Improvement Project-2 centers on building a new superconducting linear particle accelerator (Linac) at Fermilab. At the heart of the accelerator is the reference line, a critical system that defines the ideal path for the particle beam as it passes through magnets, RF cavities, and other beamline elements. Temperature stability is crucial for the reliable operation of RF components, such as mixers and filters. Fluctuations affect key performance parameters like conversion loss, isolation, and linearity. To mitigate any drift caused by ambient temperature changes, a heating plate assembly is utilized to maintain key components at a controlled temperature of 40°C. The system utilizes an aluminum 36”x36”x0.5” heat plate powered by a MOSFET-based control circuit, delivering approximately 460 W of thermal energy through a resistor array. Real-time temperature feedback is provided by a PT100 Resistance Temperature Detector (RTD), which interfaces with a Proportional–Integral–Derivative (PID) control algorithm to maintain closed-loop temperature regulation. The control signal actively modulates the gate voltage of an N channel MOSFET, dynamically adjusting power delivery in response to deviations from the temperature setpoint. Simulations and LTspice models validate the functionality and responsiveness of the circuit under varying conditions. The prototype has successfully demonstrated stable thermal control, paving the way for integration into the PIP-II infrastructure. The final design will feature an expanded resistor array, as well as communication with a PLC for continuous data acquisition and diagnostics. This work directly supports Fermilab’s broader mission by contributing to the stability and reliability of core accelerator systems, enhancing the precision of particle beam delivery for future physics experiments.

Mosher, Alexander [Fermilab]↗

One-Parameter Analytical Derivation in Modified Budyko Framework for Unsteady-State Streamflow Elasticity in Humid Catchments

The changes in climate and catchment properties have altered the hydrological processes significantly at different spatiotemporal scales around the world. In particular for finer time scales, changes in water storage, which has been commonly neglected for long-term temporal scales, may play an important role on hydrological processes. Nevertheless, few studies addressed this question in modifying the Budyko framework, with most of them remaining a steady-state assumption on catchment characteristics. Here, in this work, we derive new analytical formulas of unsteady-state streamflow elasticity in a modified Budyko framework, incorporating both storage change and one specific parameter for catchment properties. We study 78 humid catchments in the USA with simulation data of daily hydrological processes from a probability-distribution-based hydrological model (the modified HyMOD). As indicated by results, the annual storage change ratio is linearly correlated with the annual aridity index, and this relationship can be used to estimate elasticity coefficients with our formulas. The estimated elasticity coefficients perform well in simulating the annual streamflow with the power-law model. For different catchments, variability of the unsteady-state elasticity is higher than that of the steady-state elasticity. Unsteady-state streamflow coefficients show significant linear correlation with catchment properties, such as the average slope, average elevation, and catchment area. This study provides a new analytical approach to investigate the interannual stability of catchments with varying climate and catchment properties.

54 ENVIRONMENTAL SCIENCES↗

Structure optimization with stochastic density functional theory

Linear-scaling techniques for Kohn–Sham density functional theory are essential to describe the ground state properties of extended systems. Still, these techniques often rely on the localization of the density matrix or accurate embedding approaches, limiting their applicability. In contrast, stochastic density functional theory (sDFT) achieves linear- and sub-linear scaling by statistically sampling the ground state density without relying on embedding or imposing localization. In return, ground state observables, such as the forces on the nuclei, fluctuate in sDFT, making optimizing the nuclear structure a highly non-trivial problem. In this work, we combine the most recent noise-reduction schemes for sDFT with stochastic optimization algorithms to perform structure optimization within sDFT. We compare the performance of the stochastic gradient descent approach and its variations (stochastic gradient descent with momentum) with stochastic optimization techniques that rely on the Hessian, such as the stochastic Broyden–Fletcher–Goldfarb–Shanno algorithm. In conclusion, we further provide a detailed assessment of the computational efficiency and its dependence on the optimization parameters of each method for determining the ground state structure of bulk silicon with varying supercell dimensions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An experimental analysis of critical factors involved in the breakdown process of leading edge vortex flows

Experimental crosswire measurements of the flowfield above a 70 and 75 degree flat plate delta wing were performed at a Reynolds number of 250,000. Survey grids were taken normal to the platform at a series of chordwise locations for angles of attack of 20 and 30 degrees. Axial and azimuthal vorticity distributions were derived from the velocity fields. The dependence of circulation on distance from the vortex core as well as on chordwise location was examined. The effects of nondimensionalization in comparison with other experimental data was made. The circulation distribution scales with the local semispan and grows approximately linearly in the chordwise direction. For regions of the flow outside of the vortex subcore, the circulation at any chordwise station was observed to vary logarithmically with distance from the vortex axis. The circulation was also found to increase linearly with angle of incidence at a given chordwise station. A reduction in the local circulation about the vortex axis occurred at breakdown. The spanwise distribution of axial vorticity was severely altered through the breakdown region and the spanwise distribution of axial vorticity present appeared to reach a maximum immediately preceding breakdown. The local concentration of axial vorticity about the vortex axis was reduced while the magnitude of the azimuthal vorticity decreased throughout the breakdown zone. The axial vorticity components with a negative sense, found in the secondary vortex, remained unaffected by changes in wing sweep or angle of attack, in direct contrast to the positive components. The inclusion of the local wing geometry into a previously derived correlation parameter indicated that the circulation of growing leading edge vortex flows were similar at corresponding radii from the vortex axis. It was concluded that the flow over a delta wing, upstream of the breakdown regions and away from the apex and trailing edge regions, is conical. In addition, the dominating factors leading to the onset of breakdown are felt to be the local circulation of the vortex and the accompanying pressure field.

Visser, Kenneth D.↗

Sub-mm Scale Fiber Guided Deep/Vacuum Ultra-Violet Optical Source for Trapped Mercury Ion Clocks

We demonstrate the functionality of a mercury capillary lamp with a diameter in the sub-mm range and deep ultraviolet (DUV)/ vacuum ultraviolet (VUV) radiation delivery via an optical fiber integrated with the capillary. DUV spectrum control is observed by varying the fabrication parameters such as buffer gas type and pressure, capillary diameter, electrical resonator design, and temperature. We also show spectroscopic data of the 199Hg+ hyper-fine transition at 40.5GHz when applying the above fiber optical design. We present efforts toward micro-plasma generation in hollow-core photonic crystal fiber with related optical design and theoretical estimations. This new approach towards a more practical DUV optical interface could benefit trapped ion clock developments for future ultra-stable frequency reference and time-keeping applications.

lamps↗

Nonlinear Local Bending Response and Bulging Factors for Longitudinal and Circumferential Cracks in Pressurized Cylindrical Shells

Results of a geometrically nonlinear finite element parametric study to determine curvature correction factors or bulging factors that account for increased stresses due to curvature for longitudinal and circumferential cracks in unstiffened pressurized cylindrical shells are presented. Geometric parameters varied in the study include the shell radius, the shell wall thickness, and the crack length. The major results are presented in the form of contour plots of the bulging factor as a function of two nondimensional parameters: the shell curvature parameter, lambda, which is a function of the shell geometry, Poisson's ratio, and the crack length; and a loading parameter, eta, which is a function of the shell geometry, material properties, and the applied internal pressure. These plots identify the ranges of the shell curvature and loading parameters for which the effects of geometric nonlinearity are significant. Simple empirical expressions for the bulging factor are then derived from the numerical results and shown to predict accurately the nonlinear response of shells with longitudinal and circumferential cracks. The numerical results are also compared with analytical solutions based on linear shallow shell theory for thin shells, and with some other semi-empirical solutions from the literature, and limitations on the use of these other expressions are suggested.

Young, Richard D.↗

Understanding the Effect of Electron Beam Melting Scanning Strategies on the Aluminum Content and Materials State of Single Ti-6Al-4V Feedstock

Research on the additive manufacturing of metals often neglects any characterization of the composition of final parts, erroneously assuming a compositional homogeneity that matches the feedstock material. Here, the composition of electron-beam-melted Ti-6Al-4V produced through three distinct scanning strategies (linear raster and two point melting strategies, random fill and Dehoff fill) is characterized both locally and globally through energy-dispersive spectroscopy and quantitative chemical analysis. As a result of the different scanning strategies used, differing levels of preferential vaporization occur across the various parts, leading to distinct final compositions, with extremes of ~5.8 wt.% Al and ~4.8 wt.% Al. In addition, energy-dispersive spectroscopy composition maps reveal specific features in both the XY and XZ planes (with Z being the build direction) as a result of local inhomogeneous preferential vaporization. The subsequent change in composition significantly modifies the materials’ state of parts, wherein parts and local regions with higher aluminum contents lead to higher hardness levels (with a ~50 HV difference) and elastic property values and vice versa. While varying scan strategies and scan parameters are known to modify the microstructure and properties of a part, the effect on composition cannot, and should not, be neglected.

36 MATERIALS SCIENCE↗

Modeling Streamflow and Sediment Transport Under Current and Future Climates at the West Valley Site - 20231

Performance Assessment models often require inputs representing specific hydrologic processes such as streamflow and sediment flux. This paper details the development of a process-level hydrologic model and its use in a Probabilistic Performance Assessment (PPA) system-level model for the West Valley Site in western New York. For this PPA model, a watershed model was developed using SWAT (Soil and Water Assessment Tool) and used to evaluate water and sediment fluxes in nearby streams draining the Site. A particular challenge of probabilistic modeling over long periods of time is scaling parameters appropriately so that outcomes are not grossly over- or underestimated. This challenge is especially important when the spatial and temporal scales of process and system level models differ. In this instance, the SWAT watershed model and GoldSim PPA model have similar spatial scales of stream reaches and sub-catchments within a single watershed (approximately hundreds to tens of thousands of square meters, but the temporal scales differ considerably. The SWAT model calculates streamflow and sediment transport rates at the daily scale, while the PPA model requires long-term average values that are applied over hundreds to thousands of year-long periods. As such, distributions of daily climatic parameters such as average precipitation rate, temperature, wind speed, and relative humidity were developed to force a transient SWAT model that captures the dynamic behavior of the watershed in response to daily changes in the weather. Output from the model, which includes stream flow (m{sup 3}/y) and sediment transport rate (Mg/y) is then averaged over a 100-year period as an approximation for the long-term average values which can be applied in GoldSim. A single SWAT model was run 5000 times with varying climatic inputs and deterministic soil, land use, and elevation data, and the range of outputs were compiled into a distribution which can be implemented stochastically into the PPA Model. Additionally, this paper includes a high-level discussion of how changes in future climate affect weather input distributions, and a summary of how those changes in weather affect the outcomes from the SWAT model and the inputs into the PPA model. In this analysis, future climate distributions were developed based on climate research documenting the expected changes to precipitation and temperature in western New York. The 5,000 SWAT simulations were then repeated with these future climate distributions, and the resulting output was implemented into the PPA Model as hydrologic conditions under a future climate. In the current version of the PPA Model, hydrologic parameters change linearly from the 'current climate' conditions to the 'future climate' conditions over a period of 100 years, after which they remain constant. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Phase Behavior of Charged Star Block Copolymers at Fluids Interface

The phase behavior of block copolymers (BCPs) at the water-oil interface is influenced by the segmental interaction parameter (χ) and chain architecture. We synthesized a series of star block copolymers (s-BCPs) having polystyrene (PS) as core and poly(2-vinylpyridine) (P2VP) as corona. The interaction parameters of block-block (χ PS-P2VP ) and block-solvent (χ P2VP-solvent ) were varied by adjusting the pH of the aqueous solution. Lowering pH increased the fraction of quaternized-P2VP (Q-P2VP) with enhanced hydrophilicity. By transferring the equilibrated interfacial assemblies, morphologies ranging from bicontinuous films at pH of 7 and 3.1 to nanoporous and nanotubular structure at pH of 0.65 were observed. The nanoporous films formed hexagonally packed pores in s-BCP matrix, while nanotubes comprised Q-P2VP as corona and PS as core. Control over pore size, d-spacing between pores, and nanotube diameters was achieved by varying polymer concentration, molecular weight, volume fraction and arm number of s-BCPs. Large-scale nanoporous films were obtained by freeze-drying emulsions. Remarkably, the morphologies of linear BCPs were inverted, forming hexagonal-packed rigid spherical micelles with Q-P2VP as core and PS as corona in multilayer. Here, this work provides insights of phase behaviors of BCP at fluids interface and offer a facile approach to prepare nanoporous film with well-controlled pore structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Increasing Functional Diversity in a Global Land Surface Model Illustrates Uncertainties Related to Parameter Simplification

Abstract Simulations of the land surface carbon cycle typically compress functional diversity into a small set of plant functional types (PFT), with parameters defined by the average value of measurements of functional traits. In most earth system models, all wild plant life is represented by between five and 14 PFTs and a typical grid cell (≈100 × 100 km) may contain a single PFT. Model logic applied to this coarse representation of ecological functional diversity provides a reasonable proxy for the carbon cycle, but does not capture the non‐linear influence of functional traits on productivity. Here we show through simulations using the Energy Exascale Land Surface Model in 15 diverse terrestrial landscapes, that better accounting for functional diversity markedly alters predicted total carbon uptake. The shift in carbon uptake is as great as 30% and 10% in boreal and tropical regions, respectively, when compared to a single PFT parameterized with the trait means. The traits that best predict gross primary production vary based on vegetation phenology, which broadly determines where traits fall within the global distribution. Carbon uptake is more closely associated with specific leaf area for evergreen PFTs and the leaf carbon to nitrogen ratio in deciduous PFTs.

58 GEOSCIENCES↗

Elevation distributed micro-climatology data in a coastal glaciated watershed

Air temperature, ground temperature and relative humidity data were collected in a longitudinal transect of the Nooksack watershed at varying elevations from 500 to 1800 m above sea level. Data were collected by anchoring sensors from trees above winter snow levels and shaded from direct solar radiation. Paired sensors were also buried 3 cm under ground near each air temperature sensor to determine snow absence or presence. Select sites included relative humidity sensors to indicate whether precipitation was occurring. Data were collected every 3-4 h from December 2015 to Sept 2018 (with ongoing collection). Code for analysis of daily mean, minimum, maximum, and temperature change with elevation (lapse rates) are available on Github (https://doi.org/10.5281/zenodo.3239539). The sensor download and intermediate data products are available on HydroShare at (http://www.hydroshare.org/resource/222e832d3df24dea9bae9bbeb6f4219d) with publicly accessible visualization available from the Nooksack Observatory at data.cuahsi.org. Hydrologic models are generally structured with a single annual average lapse rate parameter which assumes a linear temperature gradient with elevation. The daily data (2016-2018) is used as part of ongoing studies on the non-linear dynamics and temporal variability of temperature with elevation to improve assessments of watershed function and salmon habitat.

54 ENVIRONMENTAL SCIENCES↗

A Bayesian Approach to Estimating Coupling Between Neural Components: Evaluation of the Multiple Component, Event-Related Potential (mcERP) Algorithm

Accurate measurement of single-trial responses is key to a definitive use of complex electromagnetic and hemodynamic measurements in the investigation of brain dynamics. We developed the multiple component, Event-Related Potential (mcERP) approach to single-trial response estimation. To improve our resolution of dynamic interactions between neuronal ensembles located in different layers within a cortical region and/or in different cortical regions. The mcERP model assets that multiple components defined as stereotypic waveforms comprise the stimulus-evoked response and that these components may vary in amplitude and latency from trial to trial. Maximum a posteriori (MAP) solutions for the model are obtained by iterating a set of equations derived from the posterior probability. Our first goal was to use the ANTWERP algorithm to analyze interactions (specifically latency and amplitude correlation) between responses in different layers within a cortical region. Thus, we evaluated the model by applying the algorithm to synthetic data containing two correlated local components and one independent far-field component. Three cases were considered: the local components were correlated by an interaction in their single-trial amplitudes, by an interaction in their single-trial latencies, or by an interaction in both amplitude and latency. We then analyzed the accuracy with which the algorithm estimated the component waveshapes and the single-trial parameters as a function of the linearity of each of these relationships. Extensions of these analyses to real data are discussed as well as ongoing work to incorporate more detailed prior information.

Shah, Ankoor S.↗

Parametric analysis on optimized design of hybrid solar power plants

There is increasing interest in utility-scale solar power plants with storage which can flexibly dispatch renewable energy to the grid. However, plant design possesses many degrees of freedom and non-obvious trade-offs in performance. Software tools can estimate or optimize the performance of a specific plant configuration under market and weather conditions of interest; the associated cost parameters and operating assumptions strongly influence estimates of plant performance and decisions regarding optimal sizing. We employ the National Renewable Energy Laboratory's Hybrid Optimization and Performance Platform, which incorporates optimal dispatch when evaluating plant performance, and investigate the sensitivity to weather and market conditions, operating limitations, and the presence of a capacity-based incentive. We demonstrate changes in plant performance and optimal sizing with respect to these inputs and discuss implications. Here results show that PV-with-battery designs are more profitable under our assumptions, but that designs including a concentrated solar power (CSP) system produce significantly greater annual energy; and that CSP-with-thermal energy storage designs maximizing the benefit-to-cost ratio have an input-dependent linear relationship between the CSP field solar multiple and the hours of storage as the project budget varies.

14 SOLAR ENERGY↗