Search NASA⌕ Search

SEARCH · Search NASA

Results for “distributed predictive control”

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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Peak load reduction and load shaping in HVAC and refrigeration systems in commercial buildings by using a novel lightweight dynamic priority-based control strategy

Reducing peak power demand in a building can reduce electricity expenses for the building owner and contribute to the efficiency and reliability of the electrical power grid. For the building owner, reduced expenses come from the reduction or elimination of peak power charges on electricity bills. For the power system operator, reducing peak power demand leads to a more predictable load profile and reduces stress on the electric grid system. Herein we present a computationally inexpensive, dynamic, and retrofit-deployable control strategy to effect peak load reduction and load shaping. The effectiveness of the control strategy is examined in a simulation with 80 air-conditioning units and 40 refrigeration units. The results show that a peak demand reduction of 60 kW can be achieved relative to peak demand in a typical set point–based approach. The proposed strategy was deployed in a gymnasium building with four rooftop HVAC units, where it showed over 15% peak demand (kW) reduction savings while maintaining or lowering energy consumption (in kilowatt-hours) relative to the set point–based thermostat controls.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2018

Study on 3D Internal Magnetic Field Distribution and Dynamic Mechanics of a Giant Magnetostrictive Actuator

The three-dimensional (3D) coupled magneto-elastic dynamic model of a giant magnetostrictive actuator (GMA) was built based on Maxwell’s equations, elasticity, and vibration theories. The general analytical expressions for the 3D internal magnetic field distribution of GMA were derived, and the corresponding analytical solutions were obtained. This analytical model is able to predict the internal magnetic field distribution and the output displacement of GMA system with different physical parameters under external excitations. The dynamic behaviors, such as the amplitude vs. frequency responses and the time domain waveforms of the GMA, were also studied. Analytical solutions showed that the internal magnetic field has a prominent end effect along the axial direction of a giant magnetostrictive material (GMM) rod, whereas a minimal skin effect along its radial direction, which are matched very well with the FEM results. The amplitude-frequency relationship exhibited prominent resonance and anti-resonance characteristic, which is a typical hysteresis feature of the GMA system. In order to testify the validity of the model prediction, a GMA testing system was also established. Excellent agreement was found between the experiment measurements and model predictions in the range of our research (frequency is lower than 300 Hz). The analytical results of the present study have potential applications in the fields of active vibration control and other precision machining of the GMA system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Validating corrosion models: Influence of physical properties

In chloride containing environments, two metals in physical contact can undergo galvanic corrosion limiting the lifetime of components. Being able to accurately predict galvanic corrosion damage distributions over time has not been widely presented in literature. Therefore, Finite Element Method (FEM) corrosion models were experimentally validated for two galvanic couples as a function of governing equations, environment, anode material, anode:cathode ratio, and time. Carbon steel/stainless steel (CS/SS) and zinc/stainless steel (Zn/SS) galvanic couples were exposed to NaCl solutions at room temperatures for up to 14 days. For the galvanic couples and environments, the Laplace equation with variable conductivity and reactions was sufficient to model the experimental corrosion damage. Increasing the chloride concentration for the CS/SS galvanic couple, regardless of the anode:cathode ratio, decreased the observed and modeled corrosion damage. Decreasing the anode:cathode ratio (i.e., increasing the cathode length), increased the experimentally observed and modeled corrosion damage. For small anode:cathode ratios, the governing equations deviate over long time periods with the Nernst-Plank equation being conservative. For the Zn/SS galvanic couple, the cathode length controls the dominant cathodic reduction reaction. For small cathode lengths, the hydrogen evolution reaction is dominant. For larger cathode lengths, the oxygen reduction reaction is dominant. The results are discussed with regard to the influence of solution chemistry, ohmic drop, and governing reactions. Overall, validated FEM models were presented, and the resultant models, physics, and mechanisms can be applied to other corrosion scenarios with confidence.

Carbon steel↗
NASA NTRS2004

Exhaust Simulation Testing of a Hypersonic Airbreathing Model at Transonic Speeds

An experimental study was performed to examine jet-effects for an airframe-integrated, scramjet-rocket combined-cycle vehicle configuration at transonic test conditions. This investigation was performed by testing an existing exhaust simulation wind tunnel model, known as Model 5B, in the NASA Langley 16-Ft. Transonic Tunnel. Tests were conducted at freestream Mach numbers from 0.7 to 1.2, at angles of attack from 2 to +14 degrees, and at up to seven nozzle static pressure ratio values for a set of horizontal-tail and body-flap deflections. The model aftbody, horizontal tails, and body flaps were extensively pressure instrumented to provide an understanding of jet-effects and control-surface/plume interactions, as well as for the development of analytical methodologies and calibration of computational fluid dynamic codes to predict this type of flow phenomenon. At all transonic test conditions examined, the exhaust flow at the exit of the internal nozzle was over-expanded, generating an exhaust plume that turned toward the aftbody. Pressure contour plots for the aftbody of Model 5B are presented for freestream transonic Mach numbers of 0.70, 0.95, and 1.20. These pressure data, along with shadowgraph images, indicated the impingement of an internal plume shock and at least one reflected shock onto the aftbody for all transonic conditions tested. These results also provided evidence of the highly three-dimensional nature of the aftbody exhaust flowfield. Parametric testing showed that angle-of-attack, static nozzle pressure ratio, and freestream Mach number all affected the exhaust-plume size, exhaust-flowfield shock structure, and the aftbody-pressure distribution, with Mach number having the largest effect. Integration of the aftbody pressure data showed large variations in the pitching moment throughout the transonic regime.

Huebner, Lawrence D.↗
DOE OSTI2024

Elucidating the Radiation-Induced Redox Chemistry of Plutonium Under Used Nuclear Fuel Reprocessing Conditions

Plutonium plays a critical role in the development of sustainable nuclear fuel cycles, and yet, our fundamental understanding of this element’s inherent radiation-induced redox chemistry and associated impacts on nuclear fuel cycle technologies is limited. Unanticipated changes in oxidation state distribution can influence the speciation and transport of plutonium in a given process. Control of these parameters is especially important for used nuclear fuel reprocessing technologies, wherein the separation and recovery of plutonium is typically achieved by the selective formation, maintenance, and complexation of specific oxidation states. Furthermore, plutonium’s inherent radiation-induced redox chemistry has the capacity to influence the radiolytic behavior of its complexes, the longevity of which are critical in the design of efficient and cost-effective advanced reprocessing technologies. These radiation-induced processes are unavoidable under fuel cycle conditions owing to the inherency of ionizing radiation fields to the decay of plutonium’s isotopes and to the various other radioisotopes generated by nuclear fission and neutron-capture process and the subsequent radioactive decay of their products. As such, mechanistically understanding the response of plutonium’s multiple oxidation states to multi-component ionizing radiation fields is essential for predicting the behavior of this critical element under used nuclear fuel reprocessing conditions. Here, through a combination of time-resolved (electron pulse) and steady-state (alpha and gamma) irradiation experiments complemented by quantitative, multiscale modeling calculations, we present advances in our understanding of radiation-induced plutonium redox chemistry!

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Tuning Optical Properties of Conjugated Molecules by Lewis Acids: Insights from Electronic Structure Modeling

Understanding and controlling the optoelectronic properties of organic semiconductors at the molecular level remains a challenge due to the complexity of chemical structures and intermolecular interactions. A common strategy to address this challenge is to utilize both experimental and computational approaches. In this contribution, we show that density functional theory (DFT) calculation is a useful tool to provide insights into the bonding, electron population distribution, and optical transitions of adducts between conjugated molecules and Lewis acids (CM-LA). Adduct formation leads to relevant modifications of key properties, including a red shift in optical transitions and an increase in charge carrier density and charge mobility, compared to the parent conjugated molecules. We show that electron density transfer from the CM to the LA, which was hypothesized to cause the experimental red shift in absorption spectra upon LA binding, can be quantified and interpreted by population analysis. Experimental red shifts in optical transitions for all molecular families can also be predicted by time dependent DFT calculations with different density functionals. Finally, these detailed insights help to optimize a priori design guidelines for future applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Improving Trustworthiness of Data-Driven Power Grid Contingency Analysis With Bayesian Residual Graph Neural Networks

The evolving energy landscape requires novel tools to efficiently perform contingency analysis and reliability assessment of power grids, potentially in real-time. The high computational cost of traditional power flow solvers limits their applicability in practice. Machine learning (ML) surrogates such as deep neural networks (NNs) accelerate power flow solvers computations, enabling high-order contingency analysis and real-time decision-making by learning highly nonlinear functions and integrating grid topology via graph architectures. However, (graph) NNs lack predictive power away from training data and do not provide predictive confidence estimates. Here, we present a Bayesian residual graph NN that integrates knowledge from low-fidelity data via residual training and embeds granular quantification of uncertainties, improving trustworthiness critical for high-consequence decision-making. Applying Bayesian concepts to NNs is challenging due to the high-dimensionality of both the parameter space, complicating derivation of a meaningful prior, and the output space in large grid systems, requiring enhanced techniques to assess the predicted high-dimensional uncertainties. Our contributions include: (1) Deriving a prior for fully connected and graph NNs that leverages low-fidelity data to guide mean predictions and appropriately control prior predictive uncertainty. (2) Integrating this prior within an ensembling with anchoring scheme for efficient approximate posterior inference. (3) Deriving enhanced metrics to assess accuracy of both the mean and uncertainty predictions in high dimensions, appropriately accounting for correlations propagated through graph layers. The resulting Bayesian residual graph NN is tested on a contingency analysis task for 14-bus and 118-bus grids.

24 - POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Low-Frequency Stability Analysis of Inverter-Based Islanded Multiple-Bus AC Microgrids Based on Terminal Characteristics

For system planning of three-phase inverter-based islanded ac microgrids, the low frequency instability issue caused by interactions of inverter droop controllers is a major concern. When internal control information of procured commercial inverters is unknown, impedance-based small-signal stability criteria facilitate prediction of resonances in medium and high frequency ranges, but they usually assume the grid fundamental frequency as constant and thus they are incapable of analyzing the low-frequency oscillation of the fundamental frequency in islanded microgrids. Aiming at solving this issue, this paper proposes two stability analysis methods based on terminal characteristics of inverters and passive connection network including the dynamics of the fundamental frequency for analysis of low-frequency stability in islanded multiple-bus microgrids. Based on the Component Connection Method (CCM) to systematically separate inverters from the passive connection network, a general approach is developed to model the microgrid as a multiple-input-multiple-output (MIMO) negative feedback system in the common system d-q reference frame. By applying the generalized Nyquist stability criterion (GNC) to the return-ratio and return-difference matrices of the MIMO system model, the low-frequency stability related to the fundamental frequency can be analyzed using the measured terminal characteristics of inverters. Finally, analysis and simulation of a 37-bus microgrid verify the effectiveness of the proposed stability analysis methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Basin-Size Mapping: Prediction of Metastable Polymorph Synthesizability Across TaC–TaN Alloys

The sizes of the basins of attraction on the potential energy surface are helpful indicators in determining the experimental synthesizability of metastable phases. In principle, these basins can be controlled with changes in thermodynamic conditions such as composition, pressure, and surface energy. Herein, we use random structure sampling to computationally study how alloying smoothly perturbs basin of attraction sizes. The TaC 1-x N x pseudobinary is an ideal test system given the structural and polymorphic contrast of its parent compounds and their technological relevance as epitaxial substrates for Al 1-x Ga x N. While we find limited thermodynamic stability across all computationally observed phases, random structure sampling shows a significant composition region where the rocksalt basin dominates. As such, we predict the potential for the nonequilibrium synthesis of metastable rocksalt TaC 1-x N x alloys as substrates for Al 1-x Ga x N. At higher nitrogen concentrations, other low-energy metastable polymorphs emerge that continue to retain the hexagonal close packing suitable for III-N growth. Confidence in these trends was established through uncertainty quantification of the basin sizes and energy distributions; such analysis utilized the Beta and Dirichlet distributions. In conclusion, we also find (a) polymorph basin sizes can be rationalized in terms of energetic preferences for different coordination environments; and (b) basin sizes universally shrink with increasing nitrogen content, making the system more prone to amorphous growth.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Inferring Atmospheric Release Characteristics in a Large Computer Experiment using Bayesian Adaptive Splines

An atmospheric release of hazardous material, whether accidental or intentional, can be catastrophic for those in the path of the plume. Predicting the path of a plume based on characteristics of the release (location, amount and duration) and meteorological conditions is an active research area highly relevant for emergency and long-term response to these releases. As a result, researchers have developed particle dispersion simulators to provide plume path predictions that incorporate release characteristics and meteorological conditions. However, since release characteristics and meteorological conditions are often unknown, the inverse problem is of great interest, that is, based on all the observations of the plume so far, what can be inferred about the release characteristics? This is the question we seek to answer using plume observations from a controlled release at the Diablo Canyon Nuclear Power Plant in Central California. With access to a large number of evaluations of a computationally expensive particle dispersion simulator that includes continuous and categorical inputs and spatio-temporal output, building a fast statistical surrogate model (or emulator) presents many statistical challenges, but is an essential tool for inverse modeling and sensitivity analysis. We achieve accurate emulation using Bayesian adaptive splines to model weights on empirical orthogonal functions. Here, we use this emulator as well as appropriately identifiable simulator discrepancy and observational error models to calibrate the simulator, thus finding a posterior distribution of characteristics of the release. Since the release was controlled, these characteristics are known, making it possible to compare our findings to the truth.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The impact of halo concentration on the Sunyaev Zel’dovich effect signal from massive galaxy clusters

ABSTRACT The Sunyaev Zel’dovich (SZ) effect is sensitive to the pressure of ionized gas inside galaxy clusters, which is in turn controlled largely by the gravitational potential of the cluster. Changing the concentration parameter describing the cluster mass distribution impacts the gravitational potential and thus the cluster SZ signal, with implications for cosmological and other analyses of SZ-selected clusters. We investigate the concentration–SZ relation in theory and simulations. We find that the impact of concentration on the inner SZ profile (R ≲ 0.75R200c) can be captured with standard polytropic gas models. However, we find that such models do a poor job of reproducing the outer SZ profiles (R ≳ 0.75R200c) and the relation between the integrated SZ signal, Y, and concentration. This disagreement results from a sharp truncation of the gas pressure profile near the splashback radius, likely caused by virial shocks. We develop a simple description of the truncation that leads to a good match with the simulated SZ profiles out to several R200c for clusters of varying mass and concentration, and that also accurately predicts the concentration–Y relationship. Finally, we determine how inference of the linear bias parameter and splashback radius for SZ-selected clusters can be biased by ignoring the concentration dependence of the SZ signal, finding that bias to the former is essentially negligible, while bias to the latter can be as much as 2 per cent.

Astronomy & Astrophysics↗
NASA NTRS1999

Studies on Normal and Microgravity Annular Two Phase Flows

Two-phase gas-liquid flows occur in a wide variety of situations. In addition to normal gravity applications, such flows may occur in space operations such as active thermal control systems, power cycles, and storage and transfer of cryogenic fluids. Various flow patterns exhibiting characteristic spatial and temporal distribution of the two phases are observed in two-phase flows. The magnitude and orientation of gravity with respect to the flow has a strong impact on the flow patterns observed and on their boundaries. The identification of the flow pattern of a flow is somewhat subjective. The same two-phase flow (especially near a flow pattern transition boundary) may be categorized differently by different researchers. Two-phase flow patterns are somewhat simplified in microgravity, where only three flow patterns (bubble, slug and annular) have been observed. Annular flow is obtained for a wide range of gas and liquid flow rates, and it is expected to occur in many situations under microgravity conditions. Slug flow needs to be avoided, because vibrations caused by slugs result in unwanted accelerations. Therefore, it is important to be able to accurately predict the flow pattern which exists under given operating conditions. It is known that the wavy liquid film in annular flow has a profound influence on the transfer of momentum and heat between the phases. Thus, an understanding of the characteristics of the wavy film is essential for developing accurate correlations. In this work, we review our recent results on flow pattern transitions and wavy films in microgravity.

Balakotaiah, V.↗
NASA NTRS2004

Adaptive Multi-Agent Systems for Constrained Optimization

Product Distribution (PD) theory is a new framework for analyzing and controlling distributed systems. Here we demonstrate its use for distributed stochastic optimization. First we review one motivation of PD theory, as the information-theoretic extension of conventional full-rationality game theory to the case of bounded rational agents. In this extension the equilibrium of the game is the optimizer of a Lagrangian of the (probability distribution of) the joint state of the agents. When the game in question is a team game with constraints, that equilibrium optimizes the expected value of the team game utility, subject to those constraints. The updating of the Lagrange parameters in the Lagrangian can be viewed as a form of automated annealing, that focuses the MAS more and more on the optimal pure strategy. This provides a simple way to map the solution of any constrained optimization problem onto the equilibrium of a Multi-Agent System (MAS). We present computer experiments involving both the Queen s problem and K-SAT validating the predictions of PD theory and its use for off-the-shelf distributed adaptive optimization.

Macready, William↗
DOE OSTI2025

Insights Into Thermal Runaway Mechanisms: Fast Tomography Analysis of Metal Agglomerates in Lithium-Ion Batteries

Thermal Runaway (TR) in lithium-ion batteries (LIB) is a critical technological and social concern. Whilst such events are rare, TR is characterized by uncontrollable heating leading to catastrophic failures. To deepen the understanding of the failure process and subsequently develop more accurate TR prediction models and as a result safer battery systems, we present in this work high-speed X-ray tomography for in-depth investigations of the copper current collector melting and agglomeration during TR. The melting process presents valuable real-time internal information about heat evolution during TR, previously challenging to access but crucially important for validating TR models. In this work, controlled failure studies combined with high-speed X-ray tomography were performed on two different commercial LIB models, subjecting them to both external heating and nail penetration to induce TR. Through real-time observation via high-speed tomography, followed by segmentation, rendering, and analysis, the formation of copper agglomerates was qualitatively and quantitatively characterized and visualized for the first time. Agglomerates tended to form either from the battery's outermost layers or centrally, depending on the method of TR initiation, and gives an indirect insight into the internal temperature evolution and distribution. Moreover, an initial comparative analysis between the battery models also revealed differences in agglomerate size, which has been linked to the thicker copper current collectors of one of the cell models. We further discuss the impact of larger copper agglomerates on heat distribution and safety. This study not only sheds light on the intricate dynamics of TR in LIBs but also underscores the pivotal role of 'gold-standard' imaging techniques in advancing battery safety, crucial for the robust modeling of TR and the future design of electric vehicle safety systems.

25 ENERGY STORAGE↗
DOE OSTI2022

Understanding the Hydrogeochemical Response of a Mountainous Watershed Using Integrated Surface‐Subsurface Flow and Reactive Transport Modeling

Abstract Climate change and other disturbances significantly impact hydrogeochemical exports from mountainous headwater catchments such as the Upper Colorado River Basin. Developing a mechanistic understanding of how the physical and chemical processes interact in time and space in an integrated manner is key to quantifying the future impacts of such disturbances. The hydrogeochemical response of a mountainous catchment in the 2010–2019 period is evaluated quantitatively using a high‐resolution model that simulates integrated hydrology, and transport and reactions for selected solutes and minerals. The model assumes that pyrite is present only at depth while calcite is distributed uniformly, and captures the observed C‐Q reasonably well. Distinct C‐Q dynamics are observed in an average (WY16), a wet (WY17), and a dry (WY18) water year. The model also quantifies the water fraction from surface, shallow and deep groundwater compartments using tracers, and suggests greater groundwater contributions to peak stream discharge in the dry WY18. Results demonstrate that calcium concentrations do not change significantly from year to year, while sulfate shows significant temporal variability. Pyrite dissolution is affected by the changing hydrological drivers where it is enhanced in the dry WY18; calcite dissolution supplements calcium dilution under high flow conditions. The model simulates the reaction hotspots controlled by hydrological conditions, and the spatially‐resolved results show that higher soil saturation and less snowpack occur earlier on the south‐facing side than on the north‐facing side. This is a first‐of‐its‐kind demonstration of a model that integrates hydrologic processes, including evapotranspiration, and reactive transport to enable a predictive understanding of hydrogeochemical exports.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1993

Transient performance of fan engine with water ingestion

In a continuing investigation on developing and applying codes for prediction of performance of a turbine jet engine and its components with water ingestion during flight operation, including power settings, and flight altitudes and speed changes, an attempt was made to establish the effects of water ingestion through simulation of a generic high bypass ratio engine with a generic control. In view of the large effects arising in the air compression system and the prediffuser-combustor unit during water ingestion, attention was focused on those effects and the resulting changes in engine performance. Under all conditions of operation, whether ingestion is steady or not, it became evident that water ingestion causes a fan-compressor unit to operate in a time-dependent fashion with periodic features, particularly with respect to the state of water in the span and the film in the casing clearance space, at the exit of the machine. On the other hand, the aerodynamic performance of the unit may be considered as quasi-steady once the distribution of water has attained an equilibrium state with respect to its distribution and motion. For purposes of engine simulation, the performance maps for the generic fan-compressor unit were generated based on the attainment of a quasi-steady state (meaning steady except for long-period variations in performance) during ingestion and operation over a wide enough range of rotational speeds.

Murthy, S. N. B.↗
NASA NTRS1993

A test of the applicability of independent scattering to high albedo planetary regoliths

We show that 'independent scattering' is a useful approximation for high albedo particles whose size and packing density are similar to typical particles in the lunar regolith. Laboratory measurements of the intensity and linear polarization of light scattered from a laboratory sample of glass spheres of know size and composition are compared to radiative transfer calculations of the same observable quantities. Mie scattering is integrated over the size distribution of the particles to determine the mean phase and polarization functions, or phase matrix, of the particle. Assuming that the particles scatter independently, the 'doubling method' is used to rigorously calculate multiple scattering in an optically thick layer of these anisotropically scattering particles. All of the major features 'predicted' in the calculated intensity (double peaks at small phase angles) and polarization (negative branch at small phase angles, large positive peak near 20 degrees phase, and small polarization elsewhere) are observed in the laboratory measurements, with good quantitative agreement indicated at phase angles less than 90 degrees. Even though the particles are supported by physical contact with each other in the sample, as are the particles in planetary regoliths, the independent scattering calculation yields a good approximation to both the intensity and polarization. The physical parameters input to the calculation are only the size distribution of the particles and their complex index of refraction (composition). Significant advantages of this approach are that the phase matrix is calculated from basic physical principles and both the intensity and polarization are determined simultaneously. This model may have broad applications to the interpretation of photometry, spectroscopy, and polarimetry of the icy regoliths of high albedo satellites. The intent of this effort is to perform a controlled experiment that tests the utility of the independent scattering assumption

Goguen, J. D.↗
NASA NTRS2001

Dynamic Analysis of a Building Under Rocket Engine Plume Acoustic Load

Studies have been performed to develop finite-element modeling and simulation techniques to predict the dynamic structural response of Building 4010 to the acoustic load from the plume of high-thrust rocket motors. The building is the Test Control Center and general office space for the E-complex at Stennis Space Center. It is a large single span; light-structured building located approximately 1,000 feet from the E-1 test stand. A three-dimensional shell/beam combined model of the building was built using Pro/Engineer platform and imported into Pro/Mechanica for analysis. An Equivalent Shell technique was developed to simplify the highly complex building structure so that the calculation is more efficient and accurate. A deterministic approach was used for the dynamic analysis. A pre-stressed modal analysis was performed to simulate the weight stiffening of the structure, through which about 200 modes ranging from 0 to 35 Hz were identified. In an initial dynamic frequency analysis, the maximum response over the model was found. Then the complete 3-D distributions of the displacement, as well as the stresses, were calculated through a final frequency analysis. The results were compared to a strain gage and accelerometer recordings from rocket engine tests and showed reasonable agreement.

Qian, Z.↗