Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transfer 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 181 records · Page 10

The DECADE cosmic shear project III: validation of analysis pipeline using spatially inhomogeneous data

We present the pipeline for the cosmic shear analysis of the Dark Energy Camera All Data Everywhere (DECADE) weak lensing dataset: a catalog consisting of 107 million galaxies observed by the Dark Energy Camera (DECam) in the northern Galactic cap. The catalog derives from a large number of disparate observing programs and is therefore more inhomogeneous across the sky compared to existing lensing surveys. First, we use simulated data-vectors to show the sensitivity of our constraints to different analysis choices in our inference pipeline, including sensitivity to residual systematics. Next we use simulations to validate our covariance modeling for inhomogeneous datasets. Finally, we show that our choices in the end-to-end cosmic shear pipeline are robust against inhomogeneities in the survey, by extracting relative shifts in the cosmology constraints across different subsets of the footprint/catalog and showing they are all consistent within 1σ to 2σ. This is done for forty-six subsets of the data and is carried out in a fully consistent manner: for each subset of the data, we re-derive the photometric redshift estimates, shear calibrations, survey transfer functions, the data vector, measurement covariance, and finally, the cosmological constraints. Our results show that existing analysis methods for weak lensing cosmology can be fairly resilient towards inhomogeneous datasets. This also motivates exploring a wider range of image data for pursuing such cosmological constraints.

79 ASTRONOMY AND ASTROPHYSICS↗

Impedance Analysis and PHIL Demonstration of Reactive Power Oscillations in a Wind Power Plant Using a 4-MW Wind Turbine

This paper presents impedance-based analysis, mitigation, and power-hardware-in-the-loop (PHIL) demonstration of reactive power oscillations in a wind power plant using a 4-MW Type III wind turbine drivetrain. Because such low-frequency oscillations result from interactions among slower control loops of wind turbines regulating phasor quantities—active and reactive power output of the wind turbine and the magnitude of voltages at the point of interconnection (POI)—a new type of admittance is defined in terms of phasor quantities for their analysis. The so-called power-domain admittance of a wind turbine is defined as the transfer function from the frequency and magnitude of voltages at the POI to the active and reactive power output of the turbine. The power-domain admittance responses of the 4-MW wind turbine are measured using a 7-MVA grid simulator to identify the source of the reactive power oscillations. Power-domain impedance analysis and PHIL experiments are performed to explain how a resonant mode manifests as turbine-to-turbine and plant-to-grid reactive power oscillations. It is discovered that weaker grids exhibiting high inductive impedance mitigate oscillations in the reactive power output of wind power plants; however, a higher grid impedance does not help in damping turbine-to-turbine reactive power oscillations. This paper presents a simple droop-based solution to eliminate both turbine-to-turbine and plant-to-grid reactive power oscillations.

17 WIND ENERGY↗

Real-Time Fault Location Using the Retardation Method

A new method for short-circuit fault location is proposed. The method is based on instantaneous signal measurement and its first and second derivatives, which are the novel elements of the current approach. The derivatives allow associating a precise time stamp to the occurrence of the fault. Due to retardation phenomena, the difference between the times in which a signal is registered in two detectors can be used to locate the fault. We offer several mathematical models to describe the fault. Although a description of faults in terms of a lumped circuit is useful for elucidating the methods for detecting the fault, this description will not suffice to describe the fault signal propagation; hence, a distributed models is needed, which is given in terms of the telegraph equations. Those equations were used to derive a transmission line transfer function, and an exact analytical description of the fault signal propagating in the transmission line was obtained. The analytical solution was verified both by numerical simulations and experimentally.

42 ENGINEERING↗

Bayesian Estimation of Oscillator Parameters: Toward Anomaly Detection and Cyber-Physical System Security

Cyber-physical system security presents unique challenges to conventional measurement science and technology. Anomaly detection in software-assisted physical systems, such as those employed in additive manufacturing or in DNA synthesis, is often hampered by the limited available parameter space of the underlying mechanism that is transducing the anomaly. As a result, the formulation of anomaly detection for such systems often leads to inverse or ill-posed problems, requiring statistical treatments. Here, we present Bayesian inference of unknown parameters associated with a generic actuator considered as a representative vital element of a cyber-physical system. Via a series of experimental input-output measurements, a transfer function for the actuator is obtained numerically, which serves as our model for the proposed method. Linear, nonlinear, and delayed dynamics may be assumed for the actuator response. By devising a code-based malicious signal, we study the efficacy of Bayesian inference for its potential to produce a detection, including uncertainty quantification, with a remarkably small number of input data points. Our approach should be adaptable to a variety of real-time cyber-physical anomaly detection scenarios.

47 OTHER INSTRUMENTATION↗

LiG Metrology, Correlated Error, and the Integrity of the Global Surface Air-Temperature Record

The published 95% uncertainty of the global surface air-temperature anomaly (GSATA) record through 1980 is impossibly less than the 2σ = ±0.25 °C lower limit of laboratory resolution of 1 °C/division liquid-in-glass (LiG) thermometers. The ~0.7 °C/century Joule-drift of lead- and soft-glass thermometer bulbs renders unreliable the entire historical air-temperature record through the 19th century. A circa 1900 Baudin meteorological spirit thermometer bulb exhibited intense Pb X-ray emission lines (10.55, 12.66, and 14.76 keV). Uncorrected LiG thermometer non-linearity leaves 1σ = ±0.27 °C uncertainty in land-surface air temperatures prior to 1981. The 2σ = ±0.43 °C from LiG resolution and non-linearity obscures most of the 20th century GSATA trend. Systematic sensor-measurement errors are highly pair-wise correlated, possibly across hundreds of km. Non-normal distributions of bucket and engine-intake difference SSTs disconfirm the assumption of random measurement error. Semivariogram analysis of ship SST measurements yields half the error difference mean, ±½Δε 1,2, not the error mean. Transfer-function adjustment following a change of land station air-temperature sensor eliminates measurement independence and forward-propagates the antecedent uncertainty. LiG resolution limits, non-linearity, and sensor field calibrations yield GSATA mean ±2σ RMS uncertainties of, 1900–1945, ±1.7 °C; 1946–1980, ±2.1 °C; 1981–2004, ±2.0 °C; and 2005–2010, ±1.6 °C. Finally, the 20th century (1900–1999) GSATA, 0.74 ± 1.94 °C, does not convey any information about rate or magnitude of temperature change.

54 ENVIRONMENTAL SCIENCES↗

The XFaster Power Spectrum and Likelihood Estimator for the Analysis of Cosmic Microwave Background Maps

We present the XFaster analysis package, XFaster is a fast, iterative angular power spectrum estimator based on a diagonal approximation to the quadratic Fisher matrix estimator. XFaster uses Monte Carlo simulations to compute noise biases and filter transfer functions and is thus a hybrid of both Monte Carlo and quadratic estimator methods. In contrast to conventional pseudo-C ℓ based methods, the algorithm described here requires a minimal number of simulations, and does not require them to be precisely representative of the data to estimate accurate covariance matrices for the bandpowers. The formalism works with polarization-sensitive observations and also data sets with identical, partially overlapping, or independent survey regions. The method was first implemented for the analysis of BOOMERanG data (Netterfield et al. 2002; Jones et al. 2006), and also used as part of the Planck analysis (Rocha et al. 2011). Here, we describe the full, publicly available analysis package, written in Python, as developed for the analysis of data from the 2015 flight of the SPIDER instrument (SPIDER Collaboration 2021). The package includes extensions for self-consistently estimating null spectra and for estimating fits for Galactic foreground contributions. We show results from the extensive validation of XFaster using simulations, and its application to the SPIDER data set.

79 ASTRONOMY AND ASTROPHYSICS↗

Efficacy of a portable, moderate-resolution, fast-scanning differential mobility analyzer for ambient aerosol size distribution measurements

Abstract. Ambient aerosol size distributions obtained with a compact scanning mobility analyzer, the “Spider” differential mobility analyzer (DMA), are compared to those obtained with a conventional mobility analyzer, with specific attention to the effect of mobility resolution on the measured size distribution parameters. The Spider is a 12 cm diameter radial differential mobility analyzer that spans the 10–500 nm size range with 30 s mobility scans. It achieves its compact size by operating at a nominal mobility resolution R=3 (sheath flow = 0.9 L min−1; aerosol flow = 0.3 L min−1) in place of the higher ratio of sheath flow to aerosol flow commonly used. The question addressed here is whether the lower resolution is sufficient to capture key characteristics of ambient aerosol size distributions. The Spider, operated at R=3 with 30 s up- and downscans, was co-located with a TSI 3081 long-column mobility analyzer, operated at R=10 with a 360 s sampling duty cycle. Ambient aerosol data were collected over 26 consecutive days of continuous operation, in Pasadena, CA. Over the 17–500 nm size range, the two instruments exhibit excellent correlation in the total particle number concentrations and geometric mean diameters, with regression slopes of 1.13 and 1.00, respectively. Our results suggest that particle sizing at a lower resolution than typically employed may be sufficient to obtain key properties of ambient size distributions, at least for these two moments of the size distribution. Moreover, it enables better counting statistics, as the wider transfer function for a given aerosol flow rate results in a higher counting rate.

54 ENVIRONMENTAL SCIENCES↗

Baryon fraction from the BAO amplitude: a consistent approach to parameterizing perturbation growth

Galaxy clustering constrains the baryon fraction Omega_b/Omega_m through the amplitude of baryon acoustic oscillations and the suppression of perturbations entering the horizon before recombination. This produces a different pre-recombination distribution of baryons and dark matter. After recombination, the gravitational potential responds to both components in proportion to their mass, allowing robust measurement of the baryon fraction. This is independent of new-physics scenarios altering the recombination background (e.g. Early Dark Energy). The accuracy of such measurements does, however, depend on how baryons and CDM are modeled in the power spectrum. Previous template-based splitting relied on approximate transfer functions that neglected part of information. We present a new method that embeds an extra parameter controlling the balance between baryons and dark matter in the growth terms of the perturbation equations in the CAMB Boltzmann solver. This approach captures the baryonic suppression of CDM prior to recombination, avoids inconsistencies, and yields a clean parametrization of the baryon fraction in the linear power spectrum, separating out the simple physics of growth due to the combined matter potential. We implement this framework in an analysis pipeline using Effective Field Theory of Large-Scale Structure with HOD-informed priors and validate it against noiseless LCDM and EDE cosmologies with DESI-like errors. The new scheme achieves comparable precision to previous splitting while reducing systematic biases, providing a more robust way to baryon-fraction measurements. In combination with BBN constraints on the baryon density and Alcock-Paczynski estimates of the matter density, these results strengthen the use of baryon fraction measurements to derive a Hubble constant from energy densities, with future DESI and Euclid data expected to deliver competitive constraints.

Crespi, Andrea [U. Waterloo (main); Waterloo U., I↗

Enhanced Dynamic Equivalent Identification Method of Large-Scale Power Systems Using Multiple Event

The increasing complexity of the interconnected power system makes high-fidelity dynamic simulation models computationally more intensive. To improve computation efficiency, model reduction techniques have been investigated to only preserve the dynamics in a limited area of interest (study area), while deriving equivalent representation for the external area. For this purpose, a measurement-based reduction approach using system identification techniques has been previously proposed. In that case, external areas are represented by dynamic equivalent loads using transfer function estimation. In this paper, in order to enhance the accuracy of the reduced model in preserving dynamics of the study area, multiple “grid” events of different types and at different locations are considered for identifying the parameters of the equivalent loads. Case studies are carried out in the NPCC 140-bus system. Comparison are made between multiple events training and single event training, highlighting the advantages of the proposed method in providing a better representation of the grid dynamics under different operating conditions.

Xue, Yaosuo↗

Impact of Wide-Area Oscillation Damping Control using Measurement-Driven Approach on System Separation - Saudi Grid Case Study

In any interconnected power grid, low-frequency oscillations is a major problem that can limit the power transfer capability and deteriorate power system security due to potential low-damped or even undamped oscillations. Synchronized measurements provided by PMUs enable the design and development of wide-area oscillation damping controllers (WADC) based on measurement-driven model to overcome the limitations of traditional controllers. This work focuses on the impact of a W ADC based on a measurement-driven approach on system separation for the Saudi Electricity Company (SEC) grid. A grid model is provided by SEC for this study. Modal analysis is performed to identify the oscillation mode of interest for which the controller is designed. Damping the dominant oscillation mode helps slightly to improve transient stability of the system. Considering bus frequency at different locations in SEC's system as the candidate observation signals, the optimal observation signal is chosen for the W ADC by using the FFT method. The system transfer function model is constructed by utilizing probing measurements. Then, the W ADC parameters are calculated based on the identified model. The performance of the designed W ADC is tested under different contingencies.

Altarjami, Ibrahim↗

Investigation of Nonlinear Difference-Frequency Wave Excitation on a Semisubmersible Offshore-Wind Platform with Bichromatic-Wave CFD Simulations: Preprint

The natural surge and pitch frequencies of semisubmersible offshore wind platforms are typically designed to be below the wave frequencies to avoid direct excitation. However, surge or pitch resonance can be excited by the nonlinear low-frequency loads generated by irregular incident waves. Second-order potential-flow models with added Morison drag have been found to underpredict this low-frequency excitation and response. As part of the OC6 project, Computational Fluid Dynamics (CFD) simulations were performed to enable a better understanding of the low-frequency loads and the limitations of lower-fidelity models. The focus of the current paper is to set up a computationally cost-effective CFD simulation of a fixed semisubmersible platform to investigate nonlinear difference-frequency loads and establish the corresponding uncertainty in the results. Due to the high computing cost, CFD simulations of irregular waves can be challenging. Instead, simulations were performed with bichromatic waves having a shorter repeat period. A preliminary comparison with Quadratic Transfer Functions from second-order potential-flow theory shows that CFD models consistently predict higher nonlinear wave loads at the difference frequency, likely due to flow separation and viscous drag not accounted for in potential-flow theory.

2nd order↗

Surface Charge Distribution Mapping Using an Electrostatic Probe

Undesired accumulation of charge on dielectric materials causing electrostatic discharges can be an issue in pulsed power systems and electrical systems in general. As such, an understanding of surface charge distribution, charge accumulation, and decay is required. An in-house no-contact electrostatic probe designed with a 2 mm resolution was fabricated to measure and map surface charge distribution. Designed as a contactless instrument, the probe consisting of a metal rod and low leakage amplifier circuit probes the electric potential near the surface. A scan of the surface then provides the raw data, a potential distribution. The actual probe response - i.e., the surface charge to voltage transfer function, is accounted for in post-processing. This is accomplished via an Inverse Wien Filter - a technique often applied in image processing - to deconvolve the probe response from the measured data. A commercially available electrostatic probe, the Trek 341B meter with a 3455ET probe, capable of measuring +/- 20 kV was compared to the in- house probe that is designed to cover a wider range from +/- 35 kV. A resolution better than the simple probe resolution is achieved through the distribution’s scanning voltage method and deconvolution. Applying repeated scans, surface charge decay was tracked on various polymer materials to determine the material and environment dependence; materials included were PA6, PTFE, and others. As an example of material dependence, samples charged to 20 kV at 65% humidity experience full charge decay in approximately 45 seconds for PA5 (152 mm dia.), while it took about 100 times longer for PTFE (51 mm dia.) to observe the same decay/redistribution of charge.

36 MATERIALS SCIENCE↗

A Testing Framework for Grid-Forming Resources: Preprint

This paper presents a testing framework for grid-forming (GFM) resources. It first analytically derives in frequency domain the active and reactive power response of an ideal voltage source with a reactor. These transfer functions are then used to quantify the voltage source behavior expected from a GFM inverter within subtransient time scales (5-10 cycles) following a grid disturbance. The paper also shows that the testing of a GFM inverter might require a reactor of an appropriate size between the inverter and the grid simulator used for the inverter testing. Finally, the paper presents a systematic approach for developing specifications for GFM resources using active and reactive power response frequency scans. The testing framework and results presented in this paper are demonstrated using EMT simulations of a 1-MW GFM inverter.

frequency scans↗

Dynamic Shaping of Grid Response of Multi-Machine Multi-Inverter Systems Through Grid-Forming IBRs: Preprint

We consider the problem of controlling the frequency response of weakly-coupled multi-machine multi-inverter low-inertia power systems via grid-forming inverter-based resources (IBRs). In contrast to existing methods, our approach relies on dividing the larger system into multiple strongly-coupled subsystems, without ignoring either the underlying network or approximating the subsystem response as an aggregate harmonic mean model. Rather, through a structured clustering and recursive dynamic shaping approach, the frequency response of the overall system to load perturbations is shaped appropriately. We demonstrate the proposed approach for a three-node triangular configuration and a small-scale radial network. Furthermore, for small-scale radial microgrids, we demonstrate the ability of IBRs to tune the effective transfer functions of synchronous machines. This enables us to relax the uniform turbine time-constant assumptions and widen the scope of existing synchronization results for proportionally heterogeneous machines.

frequency-shaping control↗

Impedance Scan of Inverter-Based Resources and Diesel Generator for Stability Analysis: Preprint

Impedance-based methods are widely used for power system stability analysis with inverter-based resources (IBRs), e.g., assessing dynamic interactions between the power grid and an IBR, control interactions between multiple IBRs, and the sub-synchronous oscillation and damping phenomenon. Since it is difficult to get a numerical model 100% matching with the hardware IBR, using the hardware inverter directly to obtain its output impedance has become a prominent approach nowadays. Therefore, this article presents the impedance scan using hardware IBRs, and also a hardware diesel generator as it still stays with the grid before the grid completely goes to renewable. The devices under test (DuTs) for the impedance scan includes two 3-..phi.., 480 V, 60 Hz commercial grid-forming IBRs (one of 250 kVA and another of 125 kVA rating) in series with ..delta..-Y transformers, one 3-..phi.., 480 V, 60 Hz commercial grid-following IBR (of 125 kVA rating), and a 3-..phi.., 480 V, 60 Hz commercial diesel generator (of 187.5 kVA rating). Using voltage signals perturbed with sub-, inter-, and higher harmonic components, and measuring the current response, the positive-sequence impedances are computed via an offline- based post-analysis. Moreover, best-fit transfer functions are estimated that closely resemble the measured data points of the positive-sequence impedances. Based on the observations from various outcomes of the hardware experiments, this article also provides some fundamental insights on the equivalent positive- sequence impedance of a combination of multiple hardware components by comparing the estimated and the empirically computed impedances. A comparative insight on the damping capability of the DuTs using the positive-sequence impedances of the hardware is also discussed.

grid following inverter↗

Multi-Physics Investigation of a Natural Circulation Molten Salt Micro-Reactor that Utilizes an Experimental In-Pile Device to Improve Core Physics and System Thermal-Hydraulic Performance.

INL employee PhD Dissertation - The Molten Salt Reactor (MSR) concept is a rapidly evolving Generation IV design that has recently attracted favorable attention due to the potential for reducing waste generation, realizing passive safety features, and seizing on the opportunity for cost effective economics. This thesis investigates the performance benefit of a new device invented by the doctoral candidate. The device is referred to as a Wrapped Helix around an Inclined Plane or WHIP. The WHIP is protected under a provisional patent filed with the USA Patent Office on 28 September 2021 under application number 63/261,776, BEA docket number BA-1254. The WHIP device can be located in-core or near-core to promoted enhanced thermal-hydraulics and neutronics performance. While the WHIP can be employed in a variety of solid or liquid fuel designs, this thesis investigates the device’s benefit in the application of a natural-circulation, micro-molten salt nuclear battery concept (MsNB). This thesis will specifically investigate the temperature coefficient of the MsNB fuel (FLiNaK) using novel temperature sensitivity techniques unique to Serpent particle transport code, evaluate the WHIP’s thermal-hydraulic and neutronic performance effects using both established (STAR-CCM+) and novel (Python code developed to estimate the circulation effective delayed neutron fraction, ßef f ) analytical and numerical methods, evaluate the neutron noise behavior of the MsNB and how the WHIP may alter the character of the MsNB’s transfer function, and how the WHIP affects the autonomous, load-following performance under transient power conditions using Python code developed by the candidate. The sum body of this work, in part, has been published in three journal articles as the timing of the provisional patent process has allowed. Results show that creative utilization of WHIP engineering design and function reduces reactor system volume, fuel loading, control/stability in the buoyant regime

buoyant flow↗

Advancing 1D thermo-hydraulic tools for large cryogenic facilities

The Cryogenic Division at Fermilab develops large-scale cryogenic systems for particle accelerators and superconducting magnet test facilities. To support design and diagnostics, a Python-based code was created to calculate pressure drops in components such as valves and pipes. This paper presents recent enhancements to the code, including new heat transfer functions that improve the accuracy of thermal and hydraulic predictions. The first application models pressure and temperature changes in Proton Improvement Plant PIP-II relief pipes, aiding pipe sizing and protecting relief valves. The second example analyses heat load evolution in a pipe carrying sub-atmospheric helium, helping interpret temperature sensor data and understand gas return behavior to cold compressors. These improvements significantly expand the tool's capabilities, offering a practical resource for designing and troubleshooting cryogenic systems under dynamic thermal and flow conditions.

Beckwith, R. [Fermilab]↗