Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic parameter estimation”

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 523 records · Page 29

3-D Modeling of the CFR Engine for the Investigation of Knock on Natural Gas

Abstract: Engine knock and misfire are barriers to pathways leading to higher efficiency natural gas (NG) engines, where the general tendency to knock is highly dependent on engine operating conditions and the fuel reactivity. The problem is further complicated by low emissions limits and the wide range of chemical reactivity in pipeline quality NG, making the understanding of fuel- engine interactions essential to achieving high efficiencies in spark-ignited NG engines. The present work highlights a three-dimensional computational fluid dynamics (CFD) model of the Cooperative Fuel Research (CFR) engine to investigate knock behavior when running on NG. Using a Converge CFD model of the CFR engine fuel-engine interactions that lead to knock with NG, including effects of fuel reactivity, engine operating parameters, and exhaust gas recirculation, are being investigated. A GT-Power model is used to estimate appropriate initial conditions in the Converge CFD model, including cylinder wall, piston, valves, and port surface temperatures as well as residual gas fraction. In order to validate the model, in-cylinder pressure data taken with a CFR engine is compared to model predictions while varying compression ratio, CA50, and BMEP.

CFR Engine, Converge CFD, knock, Natural Gas↗

DEMONSTRATION OF A DATA-DRIVEN PHYSICS-BASED APPROACH FOR COMPUTATIONALLY EFFICIENT CFD PREDICTION OF TWO-PHASE BUBBLY FLOW

To realize efficient computational fluid dynamics (CFD) prediction of two-phase flow, a multi-scale physics-guided data-driven approach, Feature Similarity Measurement (FSM) technique was developed for error estimation in two-phase flow simulation using coarse-mesh CFD, to achieve a comparable accuracy as fine-mesh simulations with fast-running feature. By defining physics-guided parameters and variable gradients as physical features, FSM has the capability to capture the underlying local patterns in the coarse-mesh CFD simulation. Massive low-fidelity data and respective high-fidelity data are used to explore the underlying information relevant to the main simulation errors and the effects of phenomenological scaling. By learning from previous simulation data, a surrogate model using deep feedforward neural network (DFNN) can be developed and trained to estimate the simulation error of coarse-mesh CFD. In a demonstration case of two-phase bubbly flow, the DFNN model well captured and corrected the unphysical “peaks” in the velocity and void fraction profiles near the wall in the coarse-mesh configuration, even for extrapolative predictions. The research documented supports the feasibility of the physics-guided deep learning methods for coarse mesh CFD simulations which has a potential for the efficient industrial design.

42 ENGINEERING↗

Global ocean dimethyl sulfide climatology estimated from observations and an artificial neural network

Marine dimethyl sulfide (DMS) is important to climate due to the ability of DMS to alter Earth's radiation budget. Knowledge of the global-scale distribution, seasonal variability, and sea-to-air flux of DMS is needed in order to improve understanding of atmospheric sulfur, aerosol/cloud dynamics, and albedo. Here we examine the use of an artificial neural network (ANN) to extrapolate available DMS measurements to the global ocean and produce a global climatology with monthly temporal resolution. A global database of 82 996 ship-based DMS measurements in surface waters was used along with a suite of environmental parameters consisting of latitude–longitude coordinates, time of day, time of year, solar radiation, mixed layer depth, sea surface temperature, salinity, nitrate, phosphate, and silicate. Linear regressions of DMS against the environmental parameters show that on a global-scale mixed layer depth and solar radiation are the strongest predictors of DMS. These parameters capture ~9 % and ~7 % of the raw DMS data variance, respectively. Multilinear regression can capture more of the raw data variance (~39 %) but strongly underestimates DMS in high-concentration regions. In contrast, the artificial neural network captures ~66 % of the raw data variance in our database. Like prior climatologies our results show a strong seasonal cycle in surface ocean DMS with the highest concentrations and sea-to-air fluxes in the high-latitude summertime oceans. We estimate a lower global sea-to-air DMS flux (20.12±0.43 Tg S yr –1 ) than the prior estimate based on a map interpolation method when the same gas transfer velocity parameterization is used. Our sensitivity test results show that DMS concentration does not change unidirectionally with each of the environmental parameters, which emphasizes the interactions among these parameters. The ANN model suggests that the flux of DMS from the ocean to the atmosphere will increase with global warming. Given that larger DMS fluxes induce greater cloud albedo, this corresponds to a negative climate feedback.

54 ENVIRONMENTAL SCIENCES↗

Data-Driven Robust M-LS-SVR-Based NARX Modeling for Estimation and Control of Molten Iron Quality Indices in Blast Furnace Ironmaking

Optimal operation of an industrial blast furnace (BF) ironmaking process largely depends on a reliable measurement of molten iron quality (MIQ), which is not feasible using the conventional sensors. This paper proposes a novel data-driven robust modeling method for online estimation and control of MIQ indices. First, to completely capture the nonlinear dynamics of the BF process, the nonlinear autoregressive exogenous (NARX) model is constructed for the MIQ indices. Then, considering that the standard least-squares support vector regression (LS-SVR) cannot directly cope with the multi-output problem, the multi- task transfer learning is proposed to design a novel multi-output LS-SVR (M-LS-SVR) for learning of the NARX model. Furthermore, a novel M-estimator is proposed to reduce the interference of outliers and improve the robustness of the M-LS- SVR model. Since the weights of different outlier data are properly given by the weight function, their corresponding contributions on modeling can be properly distinguished, thus a robust modeling result can be achieved. Finally, a novel multi- objective evaluation index on modeling performance is developed by comprehensively considering the root mean square error of modeling and the correlation coefficient on trend fitting, based on which the NSGA-II algorithm is used to globally optimize the model parameters. Both experiments using industrial data and industrial applications illustrate that the proposed method can eliminate the adverse effect caused by the fluctuation of data in BF process efficiently, which indicates its stronger robustness and higher accuracy. Moreover, control testing shows that the developed model can be well applied to realize data-driven control of the BF process.

blast cavities, Control↗

Affinity propagation clustering of full-field, high-spatial-dimensional measurements for robust output-only modal identification: A proof-of-concept study

Determination of the model order is a challenging problem in system identification, especially in output-only or operational modal identification where some modes are weakly excited. While existing methods such as the stabilization diagram method (spectral information) are effective, they do not scale to high-dimensional data, which is usually needed for high-fidelity characterization of structural dynamics and has been made available in the emerging full-field measurement techniques using optical methods such as photogrammetry and laser vibrometers. In this proof-of-concept study we present a new non-parametric, data-driven approach for robust output-only identification of high-spatial-dimensional modal parameters of basic structures by efficiently processing and interactively exploiting the full-field measurement (i.e., very dense spatial measurement points). Specifically, we first over-estimate the system model once, producing a pool of candidate modes associated with their modal frequencies and full-field, high-spatial-dimensional mode shapes. This is accomplished by a data-driven method termed affinity propagation clustering (APC), where the active clusters, which are the active modes in our formulations, emerge from the "message-passing" procedure and does not require a pre-determination of the cluster number (mode or model order). Next, rather than using the spectral information to distinguish the physical and spurious modes in the stabilization diagram method, we exploit and visualize the spatial, full-field mode shape associated with each candidate mode to do so. We conduct extensive experiments on basic structural models with comparisons to a few existing methods. The results indicate that the new method is computationally efficient for identifying high-spatial-dimensional modal parameters, and robust to identify weak modes by exploiting the full-field measurement. We also discuss its applicability and limitations for structures with complex geometry (shapes).

42 ENGINEERING↗

The mass splittings of SU(3) baryons within a chiral soliton model

In the framework of a chiral soliton model ({sub {chi}S}M), the complete mass splittings of SU(3) baryons are investigated, the isospin and SU(3) flavor symmetry breaking being simultaneously considered with a 'model-independent approach'. The dynamical quantities appearing in the expressions for the masses are adjusted to the masses of the baryon octet, the decuplet baryon {Sigma}* and the masses of {Theta}{sup +} and N{sup *} as inputs rather than by extracting them from a calculated self-consistent soliton profile. We estimate the masses of the baryon decuplet and antidecuplet without any further adjustable free parameter and it turns out that calculated mass values of the baryon decuplet are in very good agreement with experimental data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fault-Tolerant Decentralized Control for Large-Scale Inverter-Based Resources for Active Power Tracking

Integration of inverter-based resources (IBRs) which lack the intrinsic characteristics such as the inertial response of the traditional synchronous-generator (SG)-based sources presents a new challenge in the form of analyzing the grid stability under their presence. While the dynamic composition of IBRs differs from that of the SGs, the control objective remains similar in terms of tracking the desired active power. This letter presents a decentralized primal-dual-based fault-tolerant control framework for the power allocation in IBRs. Overall, a hierarchical control algorithm is developed with a lower level addressing the current control and the parameter estimation for the IBRs and the higher level acting as the reference power generator to the low level based on the desired active power profile. The decentralized network-based algorithm adaptively splits the desired power between the IBRs taking into consideration the health of the IBRs transmission lines. The proposed framework is tested through a simulation on the network of IBRs and the high-level controller performance is compared against the existing framework in the literature. The proposed algorithm shows significant performance improvement in the magnitude of power deviation and settling time to the nominal value under faulty conditions as compared to the algorithm in the literature.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Modeling High-Impact Weather and Climate: Lessons From a Tropical Cyclone Perspective

Although the societal impact of a weather event increases with the rarity of the event, our current ability to assess extreme events and their impacts is limited by not only rarity but also by current model fidelity and a lack of understanding of the underlying physical processes. This challenge is driving fresh approaches to assess high-impact weather and climate. Recent lessons learned in modeling high-impact weather and climate are presented using the case of tropical cyclones as an illustrative example. Through examples using the Nested Regional Climate Model to dynamically downscale large-scale climate data the need to treat bias in the driving data is illustrated. Domain size, location, and resolution are also shown to be critical and should be guided by the need to: include relevant regional climate physical processes; resolve key impact parameters; and to accurately simulate the response to changes in external forcing. The notion of sufficient model resolution is introduced together with the added value in combining dynamical and statistical assessments to fill out the parent distribution of high-impact parameters. Finally, through the example of a tropical cyclone damage index, direct impact assessments are presented as powerful tools that distill complex datasets into concise statements on likely impact, and as highly effective communication devices. Capsule: "Combining dynamical modeling of high-impact weather using traditional regional climate models with statistical techniques allows for comprehensive sampling of the full distribution, uncertainty estimation, direct assessment of impacts, and increased confidence in future changes."

climate↗

Development of a Dynamic Lidar Uncertainty Framework

As wind turbine sizes increase and wind energy expands to more complex and remote sites, remote-sensing devices such as lidars are expected to play a key role in wind resource assessment and power performance testing. The switch to remote-sensing devices represents a paradigm shift in the way the wind industry typically obtains and interprets measurement data for wind energy. For example, the measurement techniques and sources of uncertainty for a remote-sensing device are vastly different from those associated with a cup anemometer on a meteorological tower. Current IEC standards for quantifying remote sensing device uncertainty for power performance testing consider uncertainty due to mounting, calibration, and classification of the remote sensing device, among other parameters. Values of the uncertainty are typically given as a function of the mean wind speed measured by a reference device and are generally fixed, leading to climatic uncertainty values that apply to the entire measurement campaign. However, real-world experience and a consideration of the fundamentals of the measurement process have shown that lidar performance is highly dependent on atmospheric conditions, such as wind shear, turbulence, and aerosol content. At present, these conditions are not directly incorporated into the estimated uncertainty of a lidar device. In this presentation, we describe the development of a new dynamic lidar uncertainty framework that adapts to current flow conditions and more accurately represents the actual uncertainty inherent in lidar measurements under different conditions. In this new framework, sources of uncertainty are identified for estimation of the line-of-sight wind speed and reconstruction of the three-dimensional wind field. These sources are then related to physical processes caused by the atmosphere and lidar operating conditions. The framework is applied to lidar data from a field measurement site to assess the ability of the framework to predict errors in lidar-measured wind speed. The results show how uncertainty varies over time and can be used to help select data with different levels of uncertainty for different applications, for example, low uncertainty data for power performance testing versus all data for plant performance monitoring.

17 WIND ENERGY↗

Estimating epistemic uncertainty in soil parameters for nonlinear site response analyses: Introducing the Latin Hypercube Sampling technique

This study quantifies the effects of epistemic uncertainty in soil parameters on nonlinear (NL) site response analysis (SRA) results, validated against the data recorded at a well-instrumented geotechnical downhole array located in Japan. To this end, a one-dimensional soil column model of the Service Hall Array (SHA) near the Kashiwazaki-Kariwa Nuclear Power Plant (KKNPP) is developed using the finite element (FE) program LS-DYNA. The dynamic stress–strain relationship is characterized by a modified two-stage hyperbolic (MTH) NL backbone curve formulation capable of capturing soil behavior at both small- and large-shear strains. The model is then validated against the ground motion recordings to capture the model bias. The uncertainties associated with the shear-wave velocity profile (a small-strain soil property) and soil shear strength (a large-strain soil property) are incorporated in NL SRA to quantify their separate and joint randomization effects on the results. This study proposes using the Latin Hypercube Sampling (LHS) method as an efficient alternative to commonly used methods, such as Standard Monte Carlo (SMC), to account for uncertainty propagation in such reliability analysis. Both low-intensity and design-level records from the recordings at the SHA are applied to study the contribution of the small- and large-strain NL dynamic soil properties. Results from 46,200 NL FE analyses (23,100 per input motion) are presented. Measured and predicted site response, using recorded ground motions at this downhole array, is compared to assess the significance of soil parameter uncertainty on the observed ground motion dispersions. It is demonstrated that increasing the number of soft realizations and implementing higher level earthquake intensity lead to higher ground motion dispersion. Unlike past studies in randomization of Vs profiles with the SMC method, the LHS method is shown to have no significant effect on the predicted median surface response spectra and amplification factors (AFs) for this case study.

Engineering↗

Wave like signatures in aerosol optical depth and associated radiative impacts over the central Himalayan region

Doppler Lidar and Multi-Filter Rotating Shadowband Radiometer (MFRSR) observations are utilized to show wave like signatures in aerosol optical depth (AOD) during daytime boundary layer evolution over the Himalayan region. Fourier analysis depicted 60–80 min periods dominant during afternoon hours, implying that observed modulations could be plausible reason for the AOD forenoon–afternoon asymmetry which was previously reported. Inclusion of wave amplitude in diurnal variation of aerosol radiative forcing estimates showed ~40% additional warming in the atmosphere relative to mean AOD. The present observations emphasize the importance of wave induced variations in AOD and radiation budget over the site.

Aerosol optical depth↗

Better practices for inferring ecosystem water use strategy from eddy covariance data

Eddy covariance data are critical for inferring ecosystem water use strategies. Yet, such inferences are sensitive to a range of assumptions applied across studies, hindering our understanding of water use strategies within and across eddy covariance sites. A recent analysis across 151 FLUXNET2015 and AmeriFlux-FLUXNET datasets found that poor model performance was the key driver of non-robust inferences of ecosystem water use strategies. Here, we leverage this previous analysis to (i) identify the specific assumptions that improve inference model performance across most sites, (ii) explain the mechanisms behind the performance improvements, and (iii) check whether better performance improves water use inference. We find that the common practice of fitting a model to canopy conductance (G c ) derived from the evapotranspiration (ET) observations, rather than to observed ET itself, artificially amplifies data errors and degrades the model performance. Next, accounting for vegetation dynamics by applying a growing season filter or incorporating satellite LAI data improves performance, but the former practice may remove soil water stress periods. Lastly, using the leaf-to-air vapor pressure deficit (VPD l ) derived from ET observations as a model input may artificially inflate performance. Based on these results, we recommend selecting observed ET (rather than derived G c ) as the response variable, carefully accounting for vegetation dynamics, and avoiding derived VPD l as a model input; these best practices improve model performance by c. 20% and robustness by c. 80% across all eddy covariance sites. Nevertheless, the performance improvements do not always correspond to more robust inference of water use strategies, as model parameter selection and surface energy budget closure corrections still strongly influence the ecosystem water use parameter estimation in a site-specific manner.

AmeriFlux↗

Modeling gas–shell mixing in ICF with separated reactants

Mixing between fuel and shell materials in ICF implosions can affect implosion dynamics and even prevent ignition. We use data from a series of separated reactant experiments on the National Ignition Facility to calibrate and test the predictive power of gas–shell mix models. Two models are used to estimate fuel–shell mix: a Reynolds-averaged turbulence model and molecular diffusion. Minor uncertainties in capsule manufacture, experimental conditions, and values for mix model parameters produce significant variation in simulation results. Using input/output pairs from 1D simulations, we train Gaussian process surrogate models to predict experimental quantities of interest. The surrogates are used to construct posteriors for mix model parameters by marginalizing over uncertainties in capsule manufacture and experimental conditions. Mix models are calibrated with a subset of experimental data (neutron yields, ion temperature, and bang time) and tested using the remaining data. In general, both the diffusion and turbulence model correctly predict experimental DT and TT neutron yields. Despite having more free parameters, the turbulence model underpredicts ion temperature at high convergence ratio. Furthermore, the simpler diffusion model correctly predicts these temperatures, suggesting nonhydrodynamic gas–shell mix. The computational model consistently overpredicts DD neutron yield, indicating possible shortcomings outside of the mix model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Pulsed anodic arc discharge for the synthesis of carbon nanomaterials

Pulsed arc discharges can improve arc control and tailor the ablation process in the production of 1D and 2D nanostructures from carbon anodes. In this work, low-dimensional carbon nanoparticles have been generated by means of anodic arc discharge in helium atmosphere excited with a square-wave modulated signal (1–5 Hz, 10% duty cycle). The discharges were performed between two graphite electrodes with maximal peak current of 250 A and maximal voltage of 65 V. The erosion rates and conversion efficiency of the ablated anode are compared to reference samples grown in DC steady arc mode. Ablation rates in pulsed arcs are typically of the order of 1 mg/s. Combination of fast Langmuir probe diagnostics and optical emission spectroscopy provided plasma parameters of the discharges at the arc column. Ranges of 10 16 –10 17 m -3 for electron density and 0.5–2.0 eV for electron temperature are estimated. The obtained samples were characterized with Raman spectroscopy and scanning electron microscopy. The deposit on the cathode after pulsed arc consisted of carbon nanostructures such as graphene nano-platelets and carbon nanotubes. Erosion dynamics of pulsed arc discharge has been described in terms of a global model and compared to steady arc discharge. A correlation is identified among discharge regimes, optical emission patterns and ablation modes. In conclusion, pulsed anodic arc discharge is a very efficient source of carbon nanomaterials. The large control of the discharge characteristics will permit to tailor accurately the production and the properties of carbon nanotubes and graphene. Here, this deposition method is promising for the fabrication of semiconducting nanomaterials with tuneable electrical and optical properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A THREE-DIMENSIONAL MODEL OF ACTIVE REGION 7986: COMPARISON OF SIMULATIONS WITH OBSERVATIONS

In the present study, we use a forward modeling method to construct a 3D thermal structure encompassing active region 7986 of 1996 August. The extreme ultraviolet (EUV) emissions are then computed and compared with observations. The heating mechanism is inspired by a theory on Alfvén wave turbulence dissipation. The magnetic structure is built from a Solar and Heliospheric Observatory (SOHO)/MDI magnetogram and an estimated torsion parameter deduced from observations. We found that the solution to the equations in some locations is in a thermal nonequilibrium state. The time variation of the density and temperature profiles leads to time dependent emissions, which appear as thin, loop-like structures with uniform cross-section. Their timescale is consistent with the lifetime of observed coronal loops. The dynamic nature of the solution also leads to plasma flows that resemble observed coronal rain. The computed EUV emissions from the coronal part of the fan loops and the high loops compare favorably with SOHO/EIT observations in a quantitative comparison. However, the computed emission from the lower atmosphere is excessive compared to observations, a symptom common to many models. Some factors for this discrepancy are suggested, including the use of coronal abundances to compute the emissions and the neglect of atmospheric opacity effects.

79 ASTRONOMY AND ASTROPHYSICS↗

Determination of the Alfvén Speed and Plasma-beta Using the Seismology of Sunspot Umbra

For 478 centrally located sunspots observed in the optical continuum with Solar Dynamics Observatory /Helioseismic Magnetic Imager, we perform seismological diagnostics of the physical parameters of umbral photospheres. The new technique is based on the theory of slow magnetoacoustic waves in a non-isothermally stratified photosphere with a uniform vertical magnetic field. We construct a map of the weighted frequency of three-minute oscillations inside the umbra and use it for the estimation of the Alfvén speed, plasma-beta, and mass density within the umbra. We find the umbral mean Alfvén speed ranges between 10.5 and 7.5 km s{sup −1} and is negatively correlated with magnetic field strength. The umbral mean plasma-beta is found to range approximately between 0.65 and 1.15 and does not vary significantly from pores to mature sunspots. The mean density ranges between (1–6) × 10{sup −4} kg m{sup −3} and shows a strong positive correlation with magnetic field strength.

79 ASTRONOMY AND ASTROPHYSICS↗

Thermal Constraint Violations as a Result of Ambient Adjusted Transmission Line Ratings

Ambient Adjusted Ratings allow transmission line current ratings to be increased above static ratings while requiring minimal weather monitoring. However, Ambient Adjusted Ratings are still contingent on static parameters that may not match real world conditions. This poster examines how simulation methods can estimate the extent to which using Ambient Adjusted Ratings could result in the violation of a transmission line's thermal constraints.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Vortex Structure of Head Bubble in Convective Cloud Starting Plume

Abstract A developing cumulus cloud (Cu) was modeled, and dynamic, thermodynamic, and microphysical properties of an ascending head bubble reproducing the upper part of a developing Cu were investigated. The data for analysis are taken from 10-m-resolution LES of trade wind Cu under BOMEX conditions. The detection of a rising bubble is carried out using wavelet filtering of the velocity fields and microphysical fields, while a low-frequency signal of the filtering is associated with the convective-scale structure of cloud. We substantiate and discuss the representation of the bubble as a vortex ring, and estimate the parameters of this vortex ring. The simplest Hill’s vortex was chosen as a model of a vortex ring inside cloud. Analytical approximations of the radial profiles of the vertical velocity and of conservative quantities (such as total water mixing ratio and liquid water potential temperature inside and outside the bubble) are obtained. The spatial structure of these quantities is investigated using analytical expressions. Analytical models for spatial distributions of liquid water content (LWC) and adiabatic fraction (AF) are also designed and analyzed. The results demonstrate the existence of a cloud core with high values of LWC and AF up to the height of 1800 m. The horizontally averaged value of the adiabatic fraction, calculated analytically using the Hill’s vortex concept, is evaluated as 0.39, which is the typical AF value in the upper parts of such Cu. The vertical profiles of different important quantities characterizing cloud structure are presented. The analysis performed in this study allows us to conclude that a rising vortex ring plays the dominating role in formation of the thermodynamic and microphysical structure of developing Cu. Significance Statement 1) Dynamic and thermodynamic fields of a developing cumulus cloud simulated by high-resolution LES with spectral bin microphysics are separated into convective and turbulent components by means of the wavelet technique. 2) The analysis of convective component of the cloud revealed the existence of a vortex ring at the developing stage of the cloud and evaluate its parameters. 3) The analysis performed in this study allows us to conclude that a rising vortex ring plays the important role in formation of the thermodynamic and microphysical structure of developing Cu. 4) The study provides a novel insight into the cloud–environment interaction. 5) The approximating equations describing the vortex ring can be usefully applied for developing new schemes of convective parameterization.

Meteorology & Atmospheric Sciences↗