Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Global Model”

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,711 records · Page 95

Modeling the terrestrial hydrology for the global atmosphere - The future role of satellite data

A global terrestrial hydrology model has been developed for the transport and storage of moisture and heat in the ground surface layer where the hydrological parameters react to diurnal and seasonal changes in the atmosphere. The spatial and temporal variability of land surface features is considered in the model by means of large scale parameterizations. The model can be either forced by the atmosphere using conventional meteorological data or coupled to an atmospheric general circulation model (GCM) for interactive studies. The global surface is divided into 4 deg longitude by 5 deg latitude cells while the ground is represented by a thin surface layer, a bulk layer (the root zone), and a deep layer (the ground water zone). Results are presented from a seven-day global experiment which was conducted utilizing the GLAS GCM (NASA Goddard Laboratory for Atmospheric Sciences). The model has demonstrated its capability to predict, over a large region, the overall soil moisture storage and major flux exchanges with the atmosphere above and the ground water below.

Lin, J. D.↗

Classification of Clouds and Deep Convection from GEOS-5 Using Satellite Observations

With the increased resolution of global atmospheric models and the push toward global cloud resolving models, the resemblance of model output to satellite observations has become strikingly similar. As we progress with our adaptation of the Goddard Earth Observing System Model, Version 5 (GEOS-5) as a high resolution cloud system resolving model, evaluation of cloud properties and deep convection require in-depth analysis beyond a visual comparison. Outgoing long-wave radiation (OLR) provides a sufficient comparison with infrared (IR) satellite imagery to isolate areas of deep convection. We have adopted a binning technique to generate a series of histograms for OLR which classify the presence and fraction of clear sky versus deep convection in the tropics that can be compared with a similar analyses of IR imagery from composite Geostationary Operational Environmental Satellite (GOES) observations. We will present initial results that have been used to evaluate the amount of deep convective parameterization required within the model as we move toward cloud system resolving resolutions of 10- to 1-km globally.

Putman, William↗

Ocean Tide Loading Computation

September 15,2003 through May 15,2005 This grant funds the maintenance, updating, and distribution of programs for computing ocean tide loading, to enable the corrections for such loading to be more widely applied in space- geodetic and gravity measurements. These programs, developed under funding from the CDP and DOSE programs, incorporate the most recent global tidal models developed from Topex/Poscidon data, and also local tide models for regions around North America; the design of the algorithm and software makes it straightforward to combine local and global models.

Agnew, Duncan Carr↗

Variational data assimilation with a semi-Lagrangian semi-implicit global shallow-water equation model and its adjoint

An adjoint model is developed for variational data assimilation using the 2D semi-Lagrangian semi-implicit (SLSI) shallow-water equation global model of Bates et al. with special attention being paid to the linearization of the interpolation routines. It is demonstrated that with larger time steps the limit of the validity of the tangent linear model will be curtailed due to the interpolations, especially in regions where sharp gradients in the interpolated variables coupled with strong advective wind occur, a synoptic situation common in the high latitudes. This effect is particularly evident near the pole in the Northern Hemisphere during the winter season. Variational data assimilation experiments of 'identical twin' type with observations available only at the end of the assimilation period perform well with this adjoint model. It is confirmed that the computational efficiency of the semi-Lagrangian scheme is preserved during the minimization process, related to the variational data assimilation procedure.

Li, Y.↗

NASA Downscaling Project: Final Report

A team of researchers from NASA Ames Research Center, Goddard Space Flight Center, the Jet Propulsion Laboratory, and Marshall Space Flight Center, along with university partners at UCLA, conducted an investigation to explore whether downscaling coarse resolution global climate model (GCM) predictions might provide valid insights into the regional impacts sought by decision makers. Since the computational cost of running global models at high spatial resolution for any useful climate scale period is prohibitive, the hope for downscaling is that a coarse resolution GCM provides sufficiently accurate synoptic scale information for a regional climate model (RCM) to accurately develop fine scale features that represent the regional impacts of a changing climate. As a proxy for a prognostic climate forecast model, and so that ground truth in the form of satellite and in-situ observations could be used for evaluation, the MERRA and MERRA - 2 reanalyses were used to drive the NU - WRF regional climate model and a GEOS - 5 replay. This was performed at various resolutions that were at factors of 2 to 10 higher than the reanalysis forcing. A number of experiments were conducted that varied resolution, model parameterizations, and intermediate scale nudging, for simulations over the continental US during the period from 2000 - 2010. The results of these experiments were compared to observational datasets to evaluate the output.

dynamical downsizing↗

NASA Downscaling Project

A team of researchers from NASA Ames Research Center, Goddard Space Flight Center, the Jet Propulsion Laboratory, and Marshall Space Flight Center, along with university partners at UCLA, conducted an investigation to explore whether downscaling coarse resolution global climate model (GCM) predictions might provide valid insights into the regional impacts sought by decision makers. Since the computational cost of running global models at high spatial resolution for any useful climate scale period is prohibitive, the hope for downscaling is that a coarse resolution GCM provides sufficiently accurate synoptic scale information for a regional climate model (RCM) to accurately develop fine scale features that represent the regional impacts of a changing climate. As a proxy for a prognostic climate forecast model, and so that ground truth in the form of satellite and in-situ observations could be used for evaluation, the MERRA and MERRA-2 reanalyses were used to drive the NU-WRF regional climate model and a GEOS-5 replay. This was performed at various resolutions that were at factors of 2 to 10 higher than the reanalysis forcing. A number of experiments were conducted that varied resolution, model parameterizations, and intermediate scale nudging, for simulations over the continental US during the period from 2000-2010. The results of these experiments were compared to observational datasets to evaluate the output.

Ferraro, Robert↗

Electromagnetic Analysis of ITER Electron Cyclotron Emission Components and Support Structures for the Preliminary Design Review

The ITER electron cyclotron emission (ECE) diagnostic system is located at Design Shield Module (DSM) 2, Equatorial Port (EP) 9 to measure electron temperature profile and electron temperature fluctuations and also to assess nonthermal electron distributions via the oblique view. Therefore, ECE has both radial and oblique views with two couples of mirrors for two different optical views. This ECE diagnostic system shall be exposed to significant power due to unabsorbed electron cyclotron heating power in the plasma. It shall also receive large electromagnetic (EM) loads up to 100 MN/m 3 force density due to the eddy currents generated by the short 16-ms transient plasma disruption. The global EM models of EP12 and EP11 with the worst case of plasma disruption Major Disruption Downward Exponential 16 ms Category III (MD_DW_EXP16MS_CATIII) have been provided by ITER. However, the support structure for these critical components, including four mirrors, is still under the design and development phase, especially during the period of preliminary design reviews (PDRs). For those small in-vessel components that are not modeled explicitly for the PDRs, we can extract the magnetic field (B) data and flux variations (dB/dt) at those locations from the global model and then use the B and dB/dt method to calculate the eddy currents, EM forces, and torques on those components in the local submodels. In conclusion, this paper will present the detailed PDR EM analysis results of the ECE components and support structures.

Electromagnetic↗

The potential role of organics in new particle formation and initial growth in the remote tropical upper troposphere

Abstract. Global observations and model studies indicate that new particle formation (NPF) in the upper troposphere (UT) and subsequent particles supply 40 %–60 % of cloud condensation nuclei (CCN) in the lower troposphere, thus affecting the Earth's radiative budget. There are several plausible nucleation mechanisms and precursor species in this atmospheric region, which, in the absence of observational constraints, lead to uncertainties in modeled aerosols. In particular, the type of nucleation mechanism and concentrations of nucleation precursors, in part, determine the spatial distribution of new particles and resulting spatial distribution of CCN from this source. Although substantial advances in understanding NPF have been made in recent years, NPF processes in the UT in pristine marine regions are still poorly understood and are inadequately represented in global models. Here, we evaluate commonly used and state-of-the-art NPF schemes in a Lagrangian box model to assess which schemes and precursor concentrations best reproduce detailed in situ observations. Using measurements of aerosol size distributions (0.003 < Dp < 4.8 µm) in the remote marine troposphere between ∼0.18 and 13 km altitude obtained during the NASA Atmospheric Tomography (ATom) mission, we show that high concentrations of newly formed particles in the tropical UT over both the Atlantic and Pacific oceans are associated with outflow regions of deep convective clouds. We focus analysis on observations over the remote Pacific Ocean, which is a region less perturbed by continental emissions than the Atlantic. Comparing aerosol size distribution measurements over the remote Pacific with box model simulations for 32 cases shows that none of the NPF schemes most commonly used in global models, including binary nucleation of sulfuric acid and water (neutral and ion-assisted) and ternary involving sulfuric acid, water, and ammonia, are consistent with observations, regardless of precursor concentrations. Through sensitivity studies, we find that the nucleation scheme among those tested that is able to explain most consistently (21 of 32 cases) the observed size distributions is that of Riccobono et al. (2014), which involves both organic species and sulfuric acid. The method of Dunne et al. (2016), involving charged sulfuric acid–water–ammonia nucleation, when coupled with organic growth of the nucleated particles, was most consistent with the observations for 5 of 32 cases. Similarly, the neutral sulfuric acid–water–ammonia method of Napari (2002), when scaled with a tuning factor and with organic growth added, was most consistent for 6 of 32 cases. We find that to best reproduce both nucleation and growth rates, the mixing ratios of gas-phase organic precursors generally need to be at least twice that of SO2, a proxy for dimethyl sulfide (DMS). Unfortunately, we have no information on the nature of oxidized organic species that participated in NPF in this region. Global models rarely include organic-driven nucleation and growth pathways in UT conditions where globally significant NPF takes place, which may result in poor estimates of NPF and CCN abundance and contribute to uncertainties in aerosol–cloud–radiation effects. Furthermore, our results indicate that the organic aerosol precursor vapors may be important in the tropical UT above marine regions, a finding that should guide future observational efforts.

54 ENVIRONMENTAL SCIENCES↗

A global three-dimensional model of the circulation and chemistry of CFCl3, CF2Cl2, CH3CCl3, CCl4, and N2O

The use of a three-dimensional spectral model to study the tropospheric and stratospheric circulation, chemistry, and photochemistry of the CFCl3, CF2Cl2, CH3CCl3, CCl4, and N2O atmospheric species is examined. The components of the model are described. Lifetime, regional, and global trends and budgets for the species are evaluated. Calculated horizontal, vertical, and temporal distributions of the atmospheric species are compared with observations; good correlation is detected. The differences observed between calculated and observed surface distributions of CH3CCl3 and the vertical distribution of CCl4 are analyzed. The calculated global atmospheric lifetimes of CFCl3, CF2Cl2, CCl4, and N2O are 73, 210, 48, and 182 years, respectively.

Golombek, A.↗

Dynamics Explorer guest investigator

The research has focused on using the SAI auroral images as a high resolution auroral precipitation input to the USU global scale ionospheric model. From the global scale modeling viewpoint, these images offer unique spatial and temporal resolution since all prior studies have used empirical auroral models. These latter models are devoid of storm, substorm, or discrete oval features. The research focused on the problems in converting images to energy flux; using LAPU data to calibrate these energy fluxes; using the USU Time Dependent Ionospheric Model (TDIM) to look at the ionospheric consequences of this structure; and then using DE-2 in-situ observations to compare with the TDIM ionospheric parameters. In carrying out these studies, several additional investigations cropped up which were pursued to help meet the overall goals. The foremost difficulty in carrying out the TDIM modeling in conjunction with the high resolution DE auroral model was that of defining an appropriate ionospheric convection pattern. Under northward conditions this pattern is very complex. In order to study Theta aurora or in general northward IMF conditions, a new model was required. Hence, a study was completed to supply this new model to drive the TDIM as a function of the IMF. With the DE auroral model having adequate resolution to show structure on the 100's of km and all model electric fields being devoid of such structure, an investigation was pursued to find out the effects of structures in the electric field on the F-region.

Sojka, J. J.↗

Rapid Tools for an AFP Manufacturing Defects Assessment Framework

This work formulates an automated fiber placement (AFP) defects assessment framework. Such framework assumes AFP manufacturing processes are able to identify manufacturing defects via automated inspection systems and intends to provide rapid analysis tools to create a defect assessment loop. The defect assessment loop defines an automated analysis process during AFP manufacturing to minimize the number of manual repairs on the part, thus accelerating AFP manufacturing of composite structures. This defects assessment loop advantages automated inspection data to evaluate the influence of encountered AFP defects on a ply-by-ply basis. Structural evaluation is based on strength criteria using local-global finite element models. Local models affect the global part model via material property reductions. Schemes to reduce material properties using the local models are the main enablers of this framework’s automated assessment capability. Finally, we discuss the main technical challenges to realize the feasibility of this framework.

Bahamonde Jacome, Luis G.↗

The AgMIP Coordinated Climate-Crop Modeling Project (C3MP): Methods and Protocols

Climate change is expected to alter a multitude of factors important to agricultural systems, including pests, diseases, weeds, extreme climate events, water resources, soil degradation, and socio-economic pressures. Changes to carbon dioxide concentration ([CO2]), temperature, and water (CTW) will be the primary drivers of change in crop growth and agricultural systems. Therefore, establishing the CTW-change sensitivity of crop yields is an urgent research need and warrants diverse methods of investigation. Crop models provide a biophysical, process-based tool to investigate crop responses across varying environmental conditions and farm management techniques, and have been applied in climate impact assessment by using a variety of methods (White et al., 2011, and references therein). However, there is a significant amount of divergence between various crop models' responses to CTW changes (Rotter et al., 2011). While the application of a site-based crop model is relatively simple, the coordination of such agricultural impact assessments on larger scales requires consistent and timely contributions from a large number of crop modelers, each time a new global climate model (GCM) scenario or downscaling technique is created. A coordinated, global effort to rapidly examine CTW sensitivity across multiple crops, crop models, and sites is needed to aid model development and enhance the assessment of climate impacts (Deser et al., 2012). To fulfill this need, the Coordinated Climate-Crop Modeling Project (C3MP) (Ruane et al., 2014) was initiated within the Agricultural Model Intercomparison and Improvement Project (AgMIP; Rosenzweig et al., 2013). The submitted results from C3MP Phase 1 (February 15, 2013-December 31, 2013) are currently being analyzed. This chapter serves to present and update the C3MP protocols, discuss the initial participation and general findings, comment on needed adjustments, and describe continued and future development. AgMIP aims to improve substantially the climate, crop, and economic simulation tools that are used to characterize the agricultural sector, to assess future world food security under changing climate conditions, and to enhance adaptation capacity both globally and regionally. To understand better and improve the modeled crop responses, AgMIP has conducted detailed crop model intercomparisons at closely observed field sites for wheat (Asseng et al., 2013), rice (Li et al., in review), maize (Bassu et al., 2014), and sugarcane (Singels et al., 2013). A coordinated modeling exercise was one of the original motivations for AgMIP, and C3MP provides rapid estimation of crop responses to CO2, water, and temperature (CTW) changes, adding dimension and insight into the crop model intercomparisons, while facilitating interactions within the global community of modelers. C3MP also contributes a fast-track, multi-model climate sensitivity assessment for the AgMIP climate and crop modeling teams on Research Track 2 (Fig. 1), which seeks to understand the impact of projected climatic changes on crop production and food security (Rosenzweig et al., 2013; Ruane et al., 2014).

agriculture↗

A protocol and analysis of year-long simulations of global storm-resolving models and beyond

We propose a protocol to evaluate and analyze year-long simulations of global storm-resolving models (GSRMs). The proposed protocol complements an earlier 40-day simulation protocol under the DYAMOND (DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains) project to allow the analysis of the seasonal cycle and associated climatic relevant phenomena. This intercomparison aims to reveal how GSRMs, which can simulate mesoscale convective systems (MCSs) in the global domain, reproduce atmospheric large-scale structures related to convection beyond month-long simulations. The intercomparison for one-year simulations is conducted by either atmosphere-only models or atmosphere–ocean coupled models with atmospheric horizontal mesh sizes less than 5 km. We recommend the continuous four seasons from March 2020 to February 2021 as a target period for the intercomparison but with options for many groups to join more flexibly. The output variables are collected at 0.25° resolution, and archives of a small set of native grid variables are encouraged to analyze tropical cyclones and MCSs. Through the proposed global storm-resolving simulation, we will evaluate the climatological distributions of the atmospheric large-scale circulations, such as the Intertropical Convergence Zone (ITCZ), monsoon, midlatitude jets, their time evolution, and the upscale impacts on them. We present sample analyses from a one-year simulation using the 3.5 km mesh Nonhydrostatic Icosahedral Atmospheric Model (NICAM), revealing the realistic zonal contrast of tropical precipitation, no double ITCZ structure, the reasonable midlatitude jet position and intensity but a weak bias of storm track activities, and a warm bias over the Eurasia during boreal winter. We also clarify the cross-scale interaction, such as the effects of cold pools on mean precipitation over the Maritime Continent through the precipitation diurnal cycle and the effects of resolved gravity waves on midlatitude mean flows. The proposed one-year simulation protocol is referred to as the “Sendai Protocol.” This protocol is not unique or definite for evaluating GSRMs; we prospect a hierarchical set of experiments from short-term to multi-year simulations as GSRM intercomparisons.

54 ENVIRONMENTAL SCIENCES↗

Improving tide-estuary representation in MPAS-Ocean (Final Report)

The overarching goal of this research is to produce a skill-assessed tide-estuary representation, including BGC (biogeochemistry) components of the ecosystem, into the MPAS-O global climate modeling system. Specifically, in the original proposal we proposed one-way nesting, into global MPAS-O, of a high-resolution model for TES (tidal and estuary systems) including the coupling of the BGC components. We have completed the proposed tasks during the past 4 years. Specifically, after extensive consultation with the MPAS-O group during the first year, we have decided to go beyond the original proposal and to two-way nest SCHISM into MPAS-O. As a result, SCHISM-CoSiNE (with CoSiNE being the BGC model) is now fully integrated into MPAS-O as a subroutine (more precisely, an analysis member) and can be invoked as MPAS-O is running with two-way exchange of information between the TES and global models (with proper spatial interpolation done on the fly). As originally proposed, we have also developed the new finite-volume (FV) model MPAS-OI (aka ELCIRC-sub) for TES. We have conducted several tests using the coupled MPAS-SCHISM- CoSiNE model and also MPAS-OI for Columbia River, San Francisco Bay and Chesapeake Bay. Four journal publications have been produced so far as the result of the project and 1 is under review (cf. Publication section). With these achievements, we will be ready to conduct fully coupled simulations of TES and global ocean anywhere, as soon as some missing capabilities in MPAS-O (e.g. tides, improved ability to handle the boundary conditions, and streamlining of the global BGC model) are worked out. The significance of the outcome of this project cannot be overstated. It marks the first time, as far as we know, that a global climate model is tightly coupled to a bona fide nearshore and terrestrial model (SCHISM system); in fact SCHISM is the first bona fide compound flooding model that seamlessly covers the coastal ocean, estuaries, rivers and watersheds (Ye et al. 2020; Zhang et al. 2020). The coupled model will reduce the errors and uncertainties of the MPAS-O global climate model system. It also complements and extends on-going E3SM project, as the latter seeks in the long-term to include components required to simulate impacts of sea-level change and storm surge on coastal regions, including wave models and focusing resolution in coastal and storm-track regions, as an important component of the Cryosphere-ocean system. Note that a spectral wave model (the unstructured-grid version of Wave Watch III, aka Wind Wave Model) is already incorporated inside the SCHISM system and is therefore directly accessible in the coupled system.

58 GEOSCIENCES↗

Ice-shelf freshwater triggers for the Filchner–Ronne Ice Shelf melt tipping point in a global ocean–sea-ice model

Some ocean modeling studies have identified a potential tipping point from a low to a high basal melt regime beneath the Filchner–Ronne Ice Shelf (FRIS), Antarctica, with significant implications for subsequent Antarctic ice sheet mass loss. To date, investigation of the climate drivers and impacts of this possible event have been limited because ice-shelf cavities and ice-shelf melting are only now starting to be included in global climate models. Using a global ocean–sea-ice configuration of the Energy Exascale Earth System Model (E3SM) that represents both ocean circulations and melting within ice-shelf cavities, we explore freshwater triggers (iceberg melt and ice-shelf basal melt) of a transition to a high-melt regime at FRIS in a low-resolution (30 km in the Southern Ocean) global ocean–sea-ice model. We find that a realistic spatial distribution of iceberg melt fluxes is necessary to prevent the FRIS melt regime change from unrealistically occurring under historical-reanalysis-based atmospheric forcing. Further, improvement of the default parameterization for mesoscale eddy mixing significantly reduces a large regional fresh bias and weak Antarctic Slope Front structure, both of which precondition the model to melt regime change. Using two different stable model con figurations, we explore the sensitivity of FRIS melt regime change to regional ice-sheet freshwater fluxes. Through a series of sensitivity experiments prescribing incrementally increasing melt rates from the smaller, neighboring ice shelves in the eastern Weddell Sea, we demonstrate the potential for an ice-shelf melt “domino effect” should the upstream ice shelves experience increased melt rates. The experiments also reveal that modest ice-shelf melt biases in a model, especially at coarse ocean resolution where narrow continental shelf dynamics are not well resolved, can lead to an unrealistic melt regime change at downstream ice shelves. Thus, we find that remote connections between melt fluxes at different ice shelves are sensitive to baseline model conditions. Our results highlight both the potential and the peril of simulating prognostic Antarctic ice-shelf melt rates in a low-resolution global model.

58 GEOSCIENCES↗

Ice-shelf freshwater triggers for the Filchner–Ronne Ice Shelf melt tipping point in a global ocean–sea-ice model

Abstract. Some ocean modeling studies have identified a potential tipping point from a low to a high basal melt regime beneath the Filchner–Ronne Ice Shelf (FRIS), Antarctica, with significant implications for subsequent Antarctic ice sheet mass loss. To date, investigation of the climate drivers and impacts of this possible event have been limited because ice-shelf cavities and ice-shelf melting are only now starting to be included in global climate models. Using a global ocean–sea-ice configuration of the Energy Exascale Earth System Model (E3SM) that represents both ocean circulations and melting within ice-shelf cavities, we explore freshwater triggers (iceberg melt and ice-shelf basal melt) of a transition to a high-melt regime at FRIS in a low-resolution (30 km in the Southern Ocean) global ocean–sea-ice model. We find that a realistic spatial distribution of iceberg melt fluxes is necessary to prevent the FRIS melt regime change from unrealistically occurring under historical-reanalysis-based atmospheric forcing. Further, improvement of the default parameterization for mesoscale eddy mixing significantly reduces a large regional fresh bias and weak Antarctic Slope Front structure, both of which precondition the model to melt regime change. Using two different stable model configurations, we explore the sensitivity of FRIS melt regime change to regional ice-sheet freshwater fluxes. Through a series of sensitivity experiments prescribing incrementally increasing melt rates from the smaller, neighboring ice shelves in the eastern Weddell Sea, we demonstrate the potential for an ice-shelf melt “domino effect” should the upstream ice shelves experience increased melt rates. The experiments also reveal that modest ice-shelf melt biases in a model, especially at coarse ocean resolution where narrow continental shelf dynamics are not well resolved, can lead to an unrealistic melt regime change at downstream ice shelves. Thus, we find that remote connections between melt fluxes at different ice shelves are sensitive to baseline model conditions. Our results highlight both the potential and the peril of simulating prognostic Antarctic ice-shelf melt rates in a low-resolution global model.

54 ENVIRONMENTAL SCIENCES↗

A global/local analysis method for treating details in structural design

A method for analyzing global/local behavior of plate and shell structures is described. In this approach, a detailed finite element model of the local region is incorporated within a coarser global finite element model. The local model need not be nodally compatible (i.e., need not have a one-to-one nodal correspondence) with the global model at their common boundary; therefore, the two models may be constructed independently. The nodal incompatibility of the models is accounted for by introducing appropriate constraint conditions into the potential energy in a hybrid variational formulation. The primary advantage of this method is that the need for transition modeling between global and local models is eliminated. Eliminating transition modeling has two benefits. First, modeling efforts are reduced since tedious and complex transitioning need not be performed. Second, errors due to the mesh distortion, often unavoidable in mesh transitioning, are minimized by avoiding distorted elements beyond what is needed to represent the geometry of the component. The method is applied reduced to a plate loaded in tension and transverse bending. The plate has a central hole, and various hole sixes and shapes are studied. The method is also applied to a composite laminated fuselage panel with a crack emanating from a window in the panel. While this method is applied herein to global/local problems, it is also applicable to the coupled analysis of independently modeled components as well as adaptive refinement.

Aminpour, Mohammad A.↗