Search NASA⌕ Search

SEARCH · Search NASA

Results for “Optimal matching”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

High-Performance All-Polymer Solar Cells and Photodetectors Enabled by a High-Mobility n-Type Polymer and Optimized Bulk-Heterojunction Morphology

All-polymer solar cells (all-PSCs) and all-polymer photodetectors (all-PPDs) are promising for application in flexible electronics, but their performance is greatly limited by the lack of polymer acceptors and the difficulties in morphology control. Herein, we report the use of a high-mobility n-type polymer semiconductor, PNDI-DTBT, for fabricating high-performance all-PSCs and all-PPDs. By pairing PNDI-DTBT with a matched polymer donor, a prominent power conversion efficiency of 8.5% was obtained in all-PSCs, demonstrating the great potential of high-mobility n-type polymers with a conventional donor–acceptor skeleton for photovoltaic application. Moreover, all-PPDs with a low dark current density of 1.32 × 10 –8 A cm –2 at –0.1 V bias and a high specific detectivity of 4.77 × 10 12 Jones were achieved, which belong to the best results of organic photodetectors. Morphology investigations revealed that the formation of a bicontinuous interpenetrating network with optimal phase separation scale, high domain purity, and preferential vertical composition distribution in the bulk-heterojunction active layer contribute to the remarkable device performance. Furthermore, these results suggest that combining high-mobility n-type polymers and morphology optimization is fruitful to achieving high-performance all-PSCs and all-PPDs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗
NASA NTRS1975

The effects of stability augmentation on the gust response of a STOL aircraft during a curved manual approach

The task of investigating the effect of SAS was formulated as two optimal control problems for stochastic systems: (1) to compute SAS gains with a rate model in the performance index algorithm; and (2) to calculate the pilot gains and system gust response using an optimal pilot model. Both problems were solved to yield reasonable low gains for the pilot and SAS, and the lateral directional mode poles and the longitudinal short period poles could be placed accurately by the model matching algorithm. Of the longitudinal (SAS) poles achieved, the vertical rms path error was least for the unaugmented poles. The lateral rms path error was an order of magnitude larger than the vertical error and showed a plus or minus 50 percent variation with SAS. It increased with dutch roll frequency and damping, and it decreased most significantly with increased roll stability. The variation in lateral error with bank angle for curved flight was also a significant function of the augmented poles.

Porter, M. B., Jr.↗
DOE OSTI2020

Transportation analysis and related design optimization of the Fermilab high-beta 650 MHz cryomodule

The Proton Improvement Plan-II (PIP-II) at Fermi National Accelerator Laboratory (FNAL) will create a new and vastly improved accelerator, which will be the source of high-energy particles for the experiments taking place at FNAL. The new linear accelerator (LINAC) contains several types of cryomodules, which are individual particle accelerators. The last cryomodule in the LINAC will be the High-Beta 650 (HB650), which will operate at 650 MHz. Each module is approximately 15 meters in length and 1.5 meters wide, weighs 13 tonnes, and shares many design features with three of the other cryomodules. The HB650 consists of two primary sections, the outer vacuum vessel and the components that reside within it, which will be cooled to cryogenic temperatures – the cold mass. Once assembly is complete at FNAL, it will be transported to one of three places: another location at FNAL, another national laboratory within the continental United States, or to a scientific partner in Europe. Any excitation the module experiences, such as going over rough roads when on a semitrailer, can create high stresses in components and cause failure if severe enough. Additionally, any delicate components with low resonant frequencies that match the excitation spectrum could achieve resonance, potentially causing a fatigue failure by repeated flexing. To ensure the successful transportation of the HB650, a transport analysis utilizing ANSYS simulation software has been performed on major subassemblies, as well as analysis of the subassemblies as a combined system. The optimization of components, design of support structures, and overall increasing of the lowest resonant frequencies have resulted in a cryomodule that will be better equipped to handle transportation and any large impact loads that come with it.

43 PARTICLE ACCELERATORS↗
NASA NTRS1989

Predicted performance of InP solar cells in Cassegrainian and slats space concentrator arrays at 20 to 100 AM0, 80 to 100 C

Researchers have calculated the expected performance dependence of near-optimally designed shallow homojunction n+pp+InP solar cells on incident intensities up 200 AM0 and temperatures up to 100 C (373K). Both circular and rectangular cells have been considered, the former for use in a Cassegrainian concentrator array at 100 AM0, 80 to 100 C and the latter for use in a Slats type concentrator array at 20 AM0 80 to 100 C. Calculation of the temperature dependence of the performance parameters I sub sc, V sub oc, FF and eta was done by first verifying that the use of the measured temperature variation of I sub sc, of the best published value of the temperature dependence of the bandgap of InP, and of the temperature dependences of the lifetimes and mobilities of electrons and holes the same as in equivalently doped GaAs, gave calculated results that closely matched measured data on the temperature variation of I sub sc, V sub oc, and FF of four existing InP cells at 1 AM0. It was then assumed that the same temperature dependences of I sub sc, the bandgap and lifetimes and mobilities would hold in the near-optimally designed cells at the higher concentrations.

Goradia, Chandra↗
DOE OSTI2021

Fast Local Spatial Verification for Feature-Agnostic Large-Scale Image Retrieval

Images from social media can reflect diverse viewpoints, heated arguments, and expressions of creativity, adding new complexity to retrieval tasks. Researchers working on Content-Based Image Retrieval (CBIR) have traditionally tuned their algorithms to match filtered results with user search intent. However, we are now bombarded with composite images of unknown origin, authenticity, and even meaning. With such uncertainty, users may not have an initial idea of what the search query results should look like. For instance, hidden people, spliced objects, and subtly altered scenes can be difficult for a user to detect initially in a meme image, but may contribute significantly to its composition. It is pertinent to design systems that retrieve images with these nuanced relationships in addition to providing more traditional results, such as duplicates and near-duplicates — and to do so with enough efficiency at large scale. In this work, we propose a new approach for spatial verification that aims at modeling object-level regions using image keypoints retrieved from an image index, which is then used to accurately weight small contributing objects within the results, without the need for costly object detection steps. We call this method the Objects in Scene to Objects in Scene (OS2OS) score, and it is optimized for fast matrix operations, which can run quickly on either CPUs or GPUs. It performs comparably to state-of-the-art methods on classic CBIR problems (Oxford 5K, Paris 6K, and Google-Landmarks), and outperforms them in emerging retrieval tasks such as image composite matching in the NIST MFC2018 dataset and meme-style imagery from Reddit.

42 ENGINEERING↗
DOE OSTI2025

Nonadiabatic Force Matching for Alchemical Free-Energy Estimation

We propose a method to compute free-energy differences from nonadiabatic alchemical transformations by using flow-based generative models. The method, nonadiabatic force matching, hinges on estimating the dissipation along an alchemical switching process in terms of a nonadiabatic force field that can be learned through stochastic flow matching. The learned field can be used in conjunction with short-time trajectory data to evaluate upper and lower bounds on the alchemical free energy that variationally converge to the exact value if the field is optimal. Applying the method to evaluate the alchemical free energy of atomistic models shows that it can substantially reduce the simulation cost of a free-energy estimate at a negligible loss of accuracy when compared with thermodynamic integration.

Computational chemistry↗
DOE OSTI2022

An Accurate, Error-Tolerant, and Energy-Efficient Neural Network Inference Engine Based on SONOS Analog Memory

In this work, we demonstrate SONOS (silicon-oxide-nitrideoxide- silicon) analog memory arrays that are optimized for neural network inference. The devices are fabricated in a 40nm process and operated in the subthreshold regime for in-memory matrix multiplication. Subthreshold operation enables low conductances to be implemented with low error, which matches the typical weight distribution of neural networks, which is heavily skewed toward near-zero values. This leads to high accuracy in the presence of programming errors and process variations. We simulate the end-to-end neural network inference accuracy, accounting for the measured programming error, read noise, and retention loss in a fabricated SONOS array. Evaluated on the ImageNet dataset using ResNet50, the accuracy using a SONOS system is within 2.16% of floating-point accuracy without any retraining. The unique error properties and high On/Off ratio of the SONOS device allow scaling to large arrays without bit slicing, and enable an inference architecture that achieves 20 TOPS/W on ResNet50, a >10× gain in energy efficiency over state-of-the-art digital and analog inference accelerators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Evaluating pulse-shaping capabilities of next-generation pulsed power architectures

This project evaluated the pulse shaping capabilities of next-generation pulsed power (NGPP) architectures. NGPP architectures share several common attributes including multiple independent pulse-generation lines, a radial water-insulated impedance transformer, and a central vacuum insulated load region. A multi-module circuit model was developed, incorporating independent pulse-generation lines and a 2-D transmission line mesh of the radial impedance transformer to assess the effects of azimuthal asymmetry in pulse-shaped experiments. Circuit model simulations demonstrated that NGPP architectures are able to produce the the desired current pulse shapes for exemplar NGPP experiments. Additionally, the project explored automated methods for experiment design, including derivative -ree optimization and machine learning. Pulse-shaped experiments require designers to determine machine parameters that reliably produce the desired current pulse at the load, a process that typically relies on expert knowledge and iterative adjustments using the Z circuit model. Given the increased complexity of NGPP systems, this manual approach may be impractical. While the evaluated methods do not eliminate the need for manual iteration, they can reduce the time required for experiment design. Derivative-free optimization automates much of the trial-and-error process, providing a close starting point for manual adjustments or making small modifications to near-final designs. Meanwhile, deep neural network methods can generate a good qualitative match to the desired current pulse in under one second without requiring circuit model simulations.

42 ENGINEERING↗
NASA NTRS1981

A computer controlled television detector for light, X-rays and particles

A versatile, high resolution, software configurable, two-dimensional intensified vidicon quantum detector system has been developed for multiple research applications. A thin phosphor convertor allows the detection of X-rays below 20 keV and non-relativistic particles in addition to visible light, and a thicker scintillator can be used to detect X-rays up to 100 keV and relativistic particles. Faceplates may be changed to allow any active area from 1 to 40 mm square, and active areas up to 200 mm square are possible. The image is integrated in a digital memory on any software specified array size up to 4000 x 4000. The array size is selected to match the spatial resolution, which ranges from 10 to 100 microns depending on the operating mode, the active area, and the photon or particle energy. All scan and data acquisition parameters are under software control to allow optimal data collection for each application.

Kalata, K.↗
NASA NTRS1995

The NUONCE engine for LEO networks

Typical LEO networks use constellations which provide a uniform coverage. However, the demand for telecom service is dynamic and unevenly distributed around the world. We examine a more efficient and cost effective design by matching the satellite coverage with the cyclical demand for service around the world. Our approach is to use a non-uniform satellite distribution for the network. We have named this constellation design NUONCE for Non Uniform Optimal Network Communications Engine.

Lo, Martin W.↗
NASA NTRS2022

Multi-Frequency Radiometer-Based Soil Moisture Retrieval and Algorithm Parameterization Using In Situ Sites

L-band brightness temperature (TB) has been shown to provide the best sensitivity to soil moisture (SM) although C- and X-band based products offer a longer time-series from satellite-based measurements. Currently, global coverage SM is routinely produced from spaceborne measurements using all three frequency bands, but despite continued validation efforts of the products, the relative characteristics and performance of these observations have not been fully established. Therefore, this study focused on the parametrization of SM retrieval algorithms at L-, C- and X-bands using TB observations from the L-band radiometer on NASA's SM Active Passive (SMAP) mission and the C- and X-band channels of JAXA's Advanced Microwave Scanning Radiometer 2 (AMSR2) onboard the GCOM-W satellite. These can be applied in global SM retrieval algorithms using either one of the frequencies or a combination of them. The reference in situ SM data was obtained from 12 core validation sites across various land cover types around the world. The investigation highlighted the known challenges of retrieving SM from C- and X-band data compared to the higher sensitivity of the L-band data. Even with a site-specific retrieval algorithm parameterization, the mean correlation of the C- and X-band retrievals for the core validation site SM measurement were much lower than that for L-band, being 0.52 (0.54) and 0.45 (0.47) for vertical (horizontal) polarization, respectively, while for the L-band retrieval the corresponding values were 0.81 (0.77). The parameterization exercise showed that matching the C- and X-band TB measurements with an emission model was not difficult; the problem was relating the observations to SM under the influence of large roughness and vegetation effects. As a result, parameter optimization produced values for some sites that were not realistic or did not allow any practical sensitivity to SM at C- and X-band. Considering the L-band observations, the parameter optimization resulted in superior bias performance as compared to the operational SMAP product parameterization, but the sensitivity to SM changes (R and unbiased root mean square difference) did not improve markedly, or in some cases degraded at the expense of a smaller bias.

passive microwave↗
NASA NTRS2014

Optimization of a Radiative Transfer Forward Operator for Simulating SMOS Brightness Temperatures over the Upper Mississippi Basin, USA

The Soil Moisture and Ocean Salinity (SMOS) satellite mission is routinely providing global multi-angular observations of brightness temperature (TB) at both horizontal and vertical polarization with a 3-day repeat period. The assimilation of such data into a land surface model (LSM) may improve the skill of operational flood forecasts through an improved estimation of soil moisture (SM). To accommodate for the direct assimilation of the SMOS TB data, the LSM needs to be coupled with a radiative transfer model (RTM), serving as a forward operator for the simulation of multi-angular and multi-polarization top of atmosphere TBs. This study investigates the use of the Variable Infiltration Capacity (VIC) LSM coupled with the Community Microwave Emission Modelling platform (CMEM) for simulating SMOS TB observations over the Upper Mississippi basin, USA. For a period of 2 years (2010-2011), a comparison between SMOS TBs and simulations with literature-based RTM parameters reveals a basin averaged bias of 30K. Therefore, time series of SMOS TB observations are used to investigate ways for mitigating these large biases. Specifically, the study demonstrates the impact of the LSM soil moisture climatology in the magnitude of TB biases. After CDF matching the SM climatology of the LSM to SMOS retrievals, the average bias decreases from 30K to less than 5K. Further improvements can be made through calibration of RTM parameters related to the modeling of surface roughness and vegetation. Consequently, it can be concluded that SM rescaling and RTM optimization are efficient means for mitigating biases and form a necessary preparatory step for data assimilation.

SMOS↗
DOE OSTI2023

Effective sensor location for detection of change in structural dynamic response

This paper investigates the sensitivity of structural system response to the sensor location by investigating consequences of small changes in the location to the structural system response. The paper discusses how maximum observability (based on mode shape and the participation of that mode in the input provided) drives optimal location. The structural responses were investigated in terms of the g-rms response for various low-frequency inputs (pure sinusoids and real-life inputs such as an earthquake and trains). Results were then analyzed in the context of Modal Contributions Factors (MCF) and changes to the Force-to-response Transfer Functions (TRFs). A modal-matching process is first presented using a Matlab TM -based Finite Element Method (FEM) model of a cantilever beam and instrumentation to determine the location of a small mass based on three different criteria. Subsequently, the structural response is investigated using experiments and the FEM model. The accelerometer of small mass (at 1/3 height) was moved up or down to obtain changes in the structural response (TRF) to various realistic low-frequency inputs. Modal Contribution Factor (MCF) and derivative (slope) of the associated mode-shapes were correlated to the observed changes in TRFs. Results show how optimal sensor locations for detecting change in structural response can be based on the MCFs and the associated mode-shapes.

Woodall, J.↗
NASA NTRSDate not supplied

Comparisons of Mixing Efficiency for the Strut Fuel Injector Obtained from Large-Eddy and Reynolds-Averaged Simulations, and Experiments

Mixing efficiency is obtained for a strut fuel injector at hypervelocity flow conditions by using large-eddy simulations (LES), Reynolds-averaged simulations (RAS), and experiments. The injector and flow conditions have been previously investigated by using RAS and experiments as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). Because the fidelity of LES is a strong function of the grid, the mixing efficiency is obtained on two grids, the coarser of which is a factor of two coarser in each of the three dimensions with respect to the fine grid. The RAS uses the two-equation linear eddy viscosity and diffusivity modeling of Menter. In RAS, the species diffusivity model exhibits a strong dependence on the turbulent Schmidt number, which is often adjusted until some metric of engineering interest, such as the mixing efficiency, matches the experimental data. In the absence of experimental data, scale-resolving simulations, such as LES, have been proposed as surrogates for experiments that could provide the data needed to “calibrate” the turbulent Schmidt number in the RAS models. This approach is followed because LES requires significantly more computational resources (CPU, data storage, and time) than RAS, making it prohibitive for use in many engineering applications and specifically for parameter exploration or optimization. Here we examine the mixing efficiency obtained from several RAS with different values of the turbulent Schmidt number, and compare the results with those obtained from the LES and experiments. In addition, the least squares fitting approach was used to demonstrate how to obtain an estimate for the turbulent Schmidt number from LES analytically. These estimates were then used together with prior knowledge about RAS model sensitivity to select a turbulence model that was expected to best match the LES data.

LES↗
NASA NTRSDate not supplied

Comparisons of Mixing Efficiency for the Strut Fuel Injector Obtained from Large-Eddy and Reynolds-Averaged Simulations, and Experiments

Mixing efficiency is obtained for a strut fuel injector at hypervelocity flow conditions by using large-eddy simulations (LES), Reynolds-averaged simulations (RAS), and experiments. The injector and flow conditions have been previously investigated by using RAS and experiments as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). Because the fidelity of LES is a strong function of the grid, the mixing efficiency is obtained on two grids, the coarser of which is a factor of two coarser in each of the three dimensions with respect to the fine grid. The RAS uses the two-equation linear eddy viscosity and diffusivity modeling of Menter. In RAS, the species diffusivity model exhibits a strong dependence on the turbulent Schmidt number, which is often adjusted until some metric of engineering interest, such as the mixing efficiency, matches the experimental data. In the absence of experimental data, scale-resolving simulations, such as LES, have been proposed as surrogates for experiments that could provide the data needed to “calibrate” the turbulent Schmidt number in the RAS models. This approach is followed because LES requires significantly more computational resources (CPU, data storage, and time) than RAS, making it prohibitive for use in many engineering applications and specifically for parameter exploration or optimization. Here we examine the mixing efficiency obtained from several RAS with different values of the turbulent Schmidt number, and compare the results with those obtained from the LES and experiments. In addition, the least squares fitting approach was used to demonstrate how to obtain an estimate for the turbulent Schmidt number from LES analytically. These estimates were then used together with prior knowledge about RAS model sensitivity to select a turbulence model that was expected to best match the LES data.

LES↗
DOE OSTI2021

Stratigraphic Identification with Airborne Electromagnetic Methods at the Hanford Site, Washington

Stratigraphic units can influence the fate and transport of subsurface contaminants within groundwater. Units having coarse-grained sediments act as preferential flow pathways, and therefore can accelerate the transport of contaminants to reach human and ecological receptors. At legacy waste sites, detailed knowledge of subsurface stratigraphy can be used for effective monitoring and remediation planning to help minimize risk to human health and the environment. Airborne electromagnetic (AEM) methods can non-invasively provide information on kilometer-scale or larger subsurface stratigraphic features and fill informational gaps in directly sampled data from sparsely located boreholes. In this paper, we present inversion results of a 412 line-km frequency-domain AEM survey to delineate subsurface stratigraphic features at the Hanford Site, located in southeastern Washington State. The inversion was performed using a massively parallel 3D electromagnetic modeling and inversion code, where the modeling is based on solving frequency-domain Maxwell’s equations using an unstructured-mesh finite-element method and the inversion employs a Gauss-Newton optimization scheme. The results are compared to an underlying geologic framework model (GFM), built by interpolating contact depths of stratigraphic units interpreted from site borehole datasets. In areas with good borehole coverage, the inversion results show a good match with the GFM to a depth of about 60 m. Outside of these areas, the inversion results exhibit inconsistencies from the assumptions made to create the GFM, demonstrating that the AEM survey results can be used to improve the understanding of the geological conceptual model.

47 OTHER INSTRUMENTATION↗
NASA NTRS1972

Singular perturbations in the state regulator problem

Most of the results of singular perturbation theory have been concerned with initial value problems whereas optimal control problems are of two-point boundary value type. The portions of this theory applicable to the open loop state regulator problem are reviewed. For obtaining approximate solutions to the state regulator problem the method of matched asymptotic expansions is employed. This method has been developed in connection with certain fluid mechanics problems and is applicable to nonlinear as well as linear problems. It has been found in the past to be advantageous not to formulate this method generally but to apply it to each individual problem and this approach is adopted here. A general recipe for the method is given and its application is illustrated by using the method to obtain an approximate solution to a simple, specific state regulator problem.

Ardema, M. D.↗