Evaluation of flow mixing in an ARID-HV algal raceway using statistics of temporal and spatial distribution of fluid particles
Not Available
SEARCH · Search NASA
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.
Not Available
The effects of simulated microgravity on two bacterial isolates, Sphingobacterium thalpophilium and Ralstonia pickettii (formerly Burkholderia pickettii), originally recovered from water systems aboard the Mir space station were examined. These bacteria were inoculated into water, high and low concentrations of nutrient broth and subjected to simulated microgravity conditions. S. thalpophilium (which was motile and had flagella) showed no significant differences between simulated microgravity and the normal gravity control regardless of the method of enumeration and medium. In contrast, for R. pickettii (that was non-motile and lacked flagella), there were significantly higher numbers in high nutrient broth under simulated microgravity compared to normal gravity. Conversely, when R. pikkettii was inoculated into water (i.e., starvation conditions) significantly lower numbers were found under simulated microgravity compared to normal gravity. Responses to microgravity depended on the strain used (e.g., the motile strain exhibited no response to microgravity, while the non-motile strain did), the method of enumeration, and the nutrient concentration of the medium. Under oligotrophic conditions, non-motile cells may remain in geostationary orbit and deplete nutrients in their vicinity, while in high nutrient medium, resources surrounding the cell may be sufficient so that high growth is observed until nutrients becoming limiting.
Stress, radiation and microgravity cause astronauts to experience secondary immunosuppression. Spaceflight conditions enhance bacterial growth and alter antimicrobial susceptibility. Clinostats are used to model microgravity effects at lxg. In controls rotated on the vertical axis, the g-vector acts on cells as in static cultures. Salmonella enterica serovar Typhimurium virulence genes are up-regulated in modelled microgravity (MMG); a MMG regulon has been postulated. We hypothesize that the virulence of P. aeruginosa (PA) may be affected similarly by microgravity, which could be observed in MMG. This study focused on regulation of the ETA protein by PA during growth in MMG. PA103 was grown in an ETA production medium at 37 C. One series of media was inoculated with frozen cultures and grown using horizontal (MMG) or static incubation. Another series inoculated with refrigerated cultures included vertical rotating controls. Analyses included optical density (OD), agar plate counts (PC) on R2A, ETA ELISA, and protein expression by 2-D gel analyses. Growth and ETA results differed depending on inoculum, with minor effects of MMG. Proteomic analysis of 2-D gels indicate differences in protein expression with MMG. Growth and ETA results show that consistent methodology is critical when studying environmental effects. This study provides information on the relationships between environmental changes and virulence regulation, especially for flight experiments, when ground experiments are used to predict potential spaceflight effects.
High-spectral-resolution lidar (HSRL) systems provide vertical profiles of optical depth, backscatter cross-section, depolarization, and backscatter phase function. All HSRL measurements are calibrated by reference to molecular scattering, which is measured at each point in the lidar profile. Like the Raman lidar (RL), but unlike simple backscatter lidars such as the micropulse lidar (MPL), this enables the HSRL to measure backscatter cross-sections and optical depths without prior assumptions about the scattering properties of the atmosphere. The depolarization observations allow robust discrimination between ice and water clouds. Rigorous error estimates can be computed for all measurements. A very narrow angular field of view reduces multiple scattering contributions. The small field of view, coupled with a narrow optical bandwidth, nearly eliminates noise due to scattered sunlight. The laser transmitter is a diode-pumped, frequency-doubled Nd:YAG laser. Narrow-band, single-frequency operation is provided by injection seeding with a single-frequency, cw-diode-pumped diode laser. The main laser cavity is maintained in resonance with the seed laser by adjusting the cavity length to minimize the time between the Q-switch trigger and the emission of the laser pulse. The emission wavelength is tuned via temperature control of the seed laser crystal and is locked to line #1109 of the iodine absorption spectra. Locking is accomplished by minimizing the transmission through a 2-cm-long iodine absorption cell. Use of a high-repetition-rate laser and expansion of the transmitted beam through a 400-mm telescope reduces the transmitted energy density to eye-safe levels. It is possible to look directly into the output beam without hazard. The receiver and transmitter use the same afocal telescope, simplifying the maintenance of stable alignment of the transmitter and receiver although the angular FOV is only 100 μrad. The small FOV and the 4-kHz repetition rate also limit the near-field signal strength, making it possible to record continuous profiles that start at an altitude of ~100 m and extend to 30 km using photon counting detectors. The small FOV also suppresses multiple scattering contributions.
A compact fatigue testing machine that can be easily inserted into a hot cell for characterization of irradiated materials is beneficial to help determine relative fatigue performance differences between new and irradiated material. Hot cell use has been carefully considered by limiting the size and weight of the machine, simplifying sample loading and test setup for operation via master-slave manipulator, and utilizing an efficient design to minimize maintenance. Funded from a US-Japan collaborative effort, the machine has been specifically designed to help characterize titanium material specimens. These specimens are flat cantilevered beams for initial studies, possibly utilizing samples irradiated at other sources of beam. The option to test spherically shaped samples cut from the T2K vacuum window is also available. The machine is able to test a sample to $10^7$ cycles in under a week, with options to count cycles and sense material failure. The design of this machine will be presented along with current status.
PRD3000{sup TM} is a novel Personal Radiation Detector (PRD) with personnel radiation dose exposure monitoring. It is intended for First Responders, Law Enforcement, Customs Inspectors protecting critical infrastructures for detecting unexpected radioactive sources, who also need real time Hp(10) dose equivalent information. Traditional PRD devices use scintillator materials instrumented through either a photomultiplier tube or a photodiode photodetector. While the former is bulky and sensitive to magnetic fields, the latter has to compromise radiation sensitivity and energy threshold given its current noise per unit of photo-detection surface. Recently, solid state photodetectors (SiPM), based on arrays of Geiger operated diodes, have emerged as a scalable digital photodetector for photon counting. Their strong breakdown voltage temperature dependence (on the order of tens of milli-volts per K) has however limited their use for portable instruments where strong temperature gradients can be experienced, and limited power is available to temperature stabilize. The PRD3000 is based on the industry standard DMC3000 active dosimeter that complies with IEC 61526 Ed. 3 and ANSI 42.20 for direct reading personal dose equivalent meters and active personnel radiation monitors. An extension module is based on a CsI(Tl) scintillator readout by a temperature compensated SiPM. Preliminary nuclear tests combined with a measured continuous operation in excess of 240 hours from a single AAA battery cell indicate that the PRD3000 complies with the IEC 62401 Ed.2 and ANSI 42.32 without sacrificing battery life time. We present a summary of the device test results, starting with performance stability over a temperature range of - 20 deg. C to 50 deg. C, false alarm rates and dynamic response time. (authors)
Lithium (Li) metal anodes have long been counted on to meet the increasing demand for high energy, high–power rechargeable battery systems but they have been plagued by uncontrollable plating, unstable solid electrolyte interphase (SEI) formation, and the resulting low Coulombic efficiency. These problems are even aggravated under commercial levels of current density and areal capacity testing conditions. In this work, the channel–like structure of a carbonized eggplant (EP) as a stable “host” for Li metal melt infusion, is utilized. With further interphase modification of lithium fluoride (LiF), the as–formed EP–LiF composite anode maintains ≈90% Li metal theoretical capacity and can successfully suppress dendrite growth and volume fluctuation during cycling. EP–LiF offers much improved symmetric cell and full–cell cycling performance with lower and more stable overpotential under various areal capacity and elevated rate capability. Furthermore, carbonized EP serves as a light–weight high–performance current collector, achieving an average Coulombic efficiency ≈99.1% in ether–based electrolytes with 2.2 mAh cm –2 cycling areal capacity. Furthermore, the natural structure of carbonized EP will inspire further artificial designs of electrode frameworks for both Li anode and sulfur cathodes, enabling promising candidates for next–generation high–energy density batteries.
Wind applications require the ability to simulate rotating blades. To support this use-case, a novel design-order sliding mesh algorithm has been developed and deployed. The hybrid method combines the control volume finite element methodology (CVFEM) with concepts found within a discontinuous Galerkin (DG) finite element method (FEM) to manage a sliding mesh. The method has been demonstrated to be design-order for the tested polynomial basis (P=1 and P=2) and has been deployed to provide production simulation capability for a Vestas V27 (225 kW) wind turbine. Other stationary and canonical rotating ow simulations are also presented. As the majority of wind-energy applications are driving extensive usage of hybrid meshes, a foundational study that outlines near-wall numerical behavior for a variety of element topologies is presented. Results indicate that the proposed nonlinear stabilization operator (NSO) is an effective stabilization methodology to control Gibbs phenomena at large cell Peclet numbers. The study also provides practical mesh resolution guidelines for future analysis efforts. Application-driven performance and algorithmic improvements have been carried out to increase robustness of the scheme on hybrid production wind energy meshes. Specifically, the Kokkos-based Nalu Kernel construct outlined in the FY17/Q4 ExaWind milestone has been transitioned to the hybrid mesh regime. This code base is exercised within a full V27 production run. Simulation timings for parallel search and custom ghosting are presented. As the low-Mach application space requires implicit matrix solves, the cost of matrix reinitialization has been evaluated on a variety of production meshes. Results indicate that at low element counts, i.e., fewer than 100 million elements, matrix graph initialization and preconditioner setup times are small. However, as mesh sizes increase, e.g., 500 million elements, simulation time associated with \setup-up" costs can increase to nearly 50% of overall simulation time when using the full Tpetra solver stack and nearly 35% when using a mixed Tpetra- Hypre-based solver stack. The report also highlights the project achievement of surpassing the 1 billion element mesh scale for a production V27 hybrid mesh. A detailed timing breakdown is presented that again suggests work to be done in the setup events associated with the linear system. In order to mitigate these initialization costs, several application paths have been explored, all of which are designed to reduce the frequency of matrix reinitialization. Methods such as removing Jacobian entries on the dynamic matrix columns (in concert with increased inner equation iterations), and lagging of Jacobian entries have reduced setup times at the cost of numerical stability. Artificially increasing, or bloating, the matrix stencil to ensure that full Jacobians are included is developed with results suggesting that this methodology is useful in decreasing reinitialization events without loss of matrix contributions. With the above foundational advances in computational capability, the project is well positioned to begin scientific inquiry on a variety of wind-farm physics such as turbine/turbine wake interactions.
Wind applications require the ability to simulate rotating blades. To support this use-case, a novel design-order sliding mesh algorithm has been developed and deployed. The hybrid method combines the control volume finite element methodology (CVFEM) with concepts found within a discontinuous Galerkin (DG) finite element method (FEM) to manage a sliding mesh. The method has been demonstrated to be design-order for the tested polynomial basis (P=1 and P=2) and has been deployed to provide production simulation capability for a Vestas V27 (225 kW) wind turbine. Other stationary and canonical rotating flow simulations are also presented. As the majority of wind-energy applications are driving extensive usage of hybrid meshes, a foundational study that outlines near-wall numerical behavior for a variety of element topologies is presented. Results indicate that the proposed nonlinear stabilization operator (NSO) is an effective stabilization methodology to control Gibbs phenomena at large cell Peclet numbers. The study also provides practical mesh resolution guidelines for future analysis efforts. Application-driven performance and algorithmic improvements have been carried out to increase robustness of the scheme on hybrid production wind energy meshes. Specifically, the Kokkos-based Nalu Kernel construct outlined in the FY17/Q4 ExaWind milestone has been transitioned to the hybrid mesh regime. This code base is exercised within a full V27 production run. Simulation timings for parallel search and custom ghosting are presented. As the low-Mach application space requires implicit matrix solves, the cost of matrix reinitialization has been evaluated on a variety of production meshes. Results indicate that at low element counts, i.e., fewer than 100 million elements, matrix graph initialization and preconditioner setup times are small. However, as mesh sizes increase, e.g., 500 million elements, simulation time associated with "setup-up" costs can increase to nearly 50% of overall simulation time when using the full Tpetra solver stack and nearly 35% when using a mixed Tpetra- Hypre-based solver stack. The report also highlights the project achievement of surpassing the 1 billion element mesh scale for a production V27 hybrid mesh. A detailed timing breakdown is presented that again suggests work to be done in the setup events associated with the linear system. In order to mitigate these initialization costs, several application paths have been explored, all of which are designed to reduce the frequency of matrix reinitialization. Methods such as removing Jacobian entries on the dynamic matrix columns (in concert with increased inner equation iterations), and lagging of Jacobian entries have reduced setup times at the cost of numerical stability. Artificially increasing, or bloating, the matrix stencil to ensure that full Jacobians are included is developed with results suggesting that this methodology is useful in decreasing reinitialization events without loss of matrix contributions. With the above foundational advances in computational capability, the project is well positioned to begin scientific inquiry on a variety of wind-farm physics such as turbine/turbine wake interactions.
To identify an appropriate AI/ML approach for a specific problem, the best practice is to measure algorithm performance through the benchmarking process. A scientific benchmark consists of an AI-ready dataset and a reference implementation on a specific scientific question. The NASA Science Mission Directorate (SMD) has started the “Benchmark Initiative for AI/ML to create scientific benchmark datasets in three applications: 1) scientific benchmarking, which finds the best algorithm for a specific problem; 2) application benchmarking, which measures algorithm performance against a set of parameters; and 3) system benchmarking, which evaluates performance of hardware and software architecture. Currently, there are no standardized datasets available to benchmark AI/ML algorithms in the domain of space biology. In this work, we constructed two AI/ML-ready biological datasets from experiments in space-flown mice: cellular imaging and RNA-seq. First, radiation-exposed immune cells harbor DNA damage foci that can be fluorescently marked to visualize the amount of damage following exposure to ionizing radiation. However, such large datasets are difficult to analyze visually, due to imaging inconsistencies and human bias, and classical image processing approaches can fail on imaging artifacts. AI/ML are therefore exciting alternative, providing the speed of machines and the accuracy of humans. We have made this dataset available at https://registry.opendata.aws/bps_microscopy/. Second, high-throughput nucleic acid sequencing (DNA-seq, RNA-seq) has become widespread in biomedical research due to the growing availability and affordability of these assays. However, most sequencing datasets suffer from high dimensionality and low sample count. In this work, we used a generative adversarial network to synthesize a standardized, AI-ready, publicly available benchmark dataset for space biology RNA-seq data with sufficient space-flown and ground control mouse liver samples from NASA GeneLab. This dataset is available at https://registry.opendata.aws/bps_rnaseq/. These datasets are now fully open the Space Biology community to test their favorite AI/ML approaches.
RNA abundance quantification has become routine and affordable thanks to high-throughput “short-read” technologies that provide accurate molecule counts at the gene level. Similarly accurate and affordable quantification of definitive full-length, transcript isoforms has remained a stubborn challenge, despite its obvious biological significance across a wide range of problems. “Long-read” sequencing platforms now produce data-types that can, in principle, drive routine definitive isoform quantification. However some particulars of contemporary long-read datatypes, together with isoform complexity and genetic variation, present bioinformatic challenges. We show here, using ONT data, that fast and accurate quantification of long-read data is possible and that it is improved by exome capture. To perform quantifications we developed lr-kallisto, which adapts the kallisto bulk and single-cell RNA-seq quantification methods for long-read technologies.
The datasets include the capacity, energy, and ancillary service values of a marginal kilowatt-hour (kWh) of generic, daily, shiftable building flexibility as a presumed market entrant in the 2030 U.S. power systems. The results should be interpreted along with the caveats listed in "Valuation of Building Flexibility: Grid Service Value for Initial Market Entrants in Projected 2030 United States Power Systems", which also documents the methods used for the study. The factors examined in the study include: grid scenario, region, original usage hour in the day (local time), and building flexibility parameters - efficiency, dissipation, and shifting window include max pre-shift and max post-shift. Filenames in the datasets: The filenames contain information on the grid scenario, and the building flexibility efficiency and dissipation. For example: MidCase_2030_efficiency1.25_dissipation0.05_value.csv contains all results under Mid RE 2020, efficiency = 1.25, and dissipation = 0.05. * MidCase = Mid RE Each .csv file contains: region: The location of the building flexibility, aligned with U.S. Energy Information Administration (EIA) National Energy Modeling System (NEMS) Electricity Market Module regions. max_pre_shift: Part of building flexibility shifting window parameter, indicates the max number of hours the building flexibility can shift earlier. max_post_shift: Part of building flexibility shifting window parameter, indicates the max number of hours the building flexibility can shift later. local_datetime: Original datetime of 1kWh of building energy consumption. local_orig_h: Original usage hour (1-24) of building energy consumption, ignores daylight saving. local_shift_to: The datetime to when building energy consumption shifts, if shifting happens. If no shifting occurs, this cell is left blank. energy: Net energy value of energy shifting. capacity: Net capacity value of energy shifting. shifting_value: Net energy plus capacity value of energy shifting. spin: Spin reserve value at the original datetime of consumption. flex: Flexible reserve value at the original datetime of consumption. reg: Regulation reserve value at the original datetime of consumption. total_profit: Assuming the building flexibility is capable of providing energy, capacity, and ancillary services, total profit is the maximum value of energy shifting value, spin reserve value, flexible reserve value, and regulation reserve value - one of the four choices at any given hour - because we do not allow its value to be double-counted.
This paper focuses on the parallel performance of two high-performance aerodynamic simulation packages on the newly installed NASA Columbia supercomputer. These packages include both a high-fidelity, unstructured, Reynolds-averaged Navier-Stokes solver, and a fully-automated inviscid flow package for cut-cell Cartesian grids. The complementary combination of these two simulation codes enables high-fidelity characterization of aerospace vehicle design performance over the entire flight envelope through extensive parametric analysis and detailed simulation of critical regions of the flight envelope. Both packages. are industrial-level codes designed for complex geometry and incorpor.ats. CuStomized multigrid solution algorithms. The performance of these codes on Columbia is examined using both MPI and OpenMP and using both the NUMAlink and InfiniBand interconnect fabrics. Numerical results demonstrate good scalability on up to 2016 CPUs using the NUMAIink4 interconnect, with measured computational rates in the vicinity of 3 TFLOP/s, while InfiniBand showed some performance degradation at high CPU counts, particularly with multigrid. Nonetheless, the results are encouraging enough to indicate that larger test cases using combined MPI/OpenMP communication should scale well on even more processors.
Power subsystem cost/weight tradeoffs are discussed for the Venus probe spacecraft. The cost estimations of power subsystem units were based upon DSCS-2, DSP, and Pioneer 10 and 11 hardware design and development and manufacturing experience. Parts count and degree of modification of existing hardware were factored into the estimate of manufacturing and design and development costs. Cost data includes sufficient quantities of units to equip probe bus and orbiter versions. It was based on the orbiter complement of equipment, but the savings in fewer slices for the probe bus balance the cost of the different probe bus battery. The preferred systems for the Thor/Delta and for the Atlas/Centaur are discussed. The weights of the candidate designs were based upon slice or tray weights for functionally equivalent circuitry measured on existing hardware such as Pioneers 10 and 11, Intelsat 3, DSCS-2, or DSP programs. Battery weights were based on measured cell weight data adjusted for case weight or off-the-shelf battery weights. The solar array weight estimate was based upon recent hardware experience on DSCS-2 and DSP arrays.
Microbial colonization of metal surfaces or adhesion of bacteria on metal surfaces is a crucial step in initiating microbiologically influenced corrosion. As the substrate is essential in the development of a biofilm, an understanding of how substrate properties affect adhesion of bacterial cells, including the mechanism of bacterial adhesion, may assist in designing or modifying substrates to inhibit bacterial adhesion and prevent the corrosion of materials. A study was carried out to investigate the effect of nickel on the adhesion and corrosion ability of stainless steels with the same chromium content in nutrient-rich simulated seawater containing Pseudomonas aeruginosa. These stainless steels were 201, 316, 304 and 430 stainless steels. In addition, a carbon steel sample with fully nickel coated was included for comparison. The hydrophobicity of different materials was evaluated by measuring the contact angles made by the water, glycerol and formaldehyde. The number of bacteria adhering to the material surface after short-term immersion was counted by fluorescence microscopy. For corrosion study, the corrosion rates and pitting morphology were studied after 42 days of immersion. The corrosion behavior of the samples was also observed by scanning electron microscope (SEM) coupled with energy-dispersive x-ray spectroscopy (EDX) and inductively coupled plasma mass spectrometry (ICPMS). It was found that the nickel-coated sample had the highest level of bacteria attachment and offered poor corrosion resistance under artificial seawater with Pseudomonas aeruginosa strain. Higher nickel content in stainless steel resulted in a higher level of bacteria attachment. It was found that stainless steel with molybdenum exhibited better corrosion resistance in artificial seawater containing bacteria.
We report microfabricated Coulter counters are attractive for point of care (POC) applications since they are label free and compact. However, these approaches inherently suffer from a trade off between sample throughput and sensitivity. The counter measures a change in impedance due to displaced fluid volume by passing cells, and thus the counter's signal increases with the fraction of the sensing volume displaced. Reducing the size of the sensing region requires reductions in volumetric throughput in the absence of increased hydraulic pressure and sensor bandwidth. The risk of mechanical clog formation, rendering the counter inoperable, increases markedly with reductions in the size of the constriction aperture. We present here a microfluidic coplanar Coulter counter device design that overcomes the problem of constriction clogging while capable of operating in microfluidic channels filled entirely with highly conductive sample. The device utilizes microfabricated planar electrodes projecting into one side of the microfluidic channel and is easily integrated with upstream electronic, hydrodynamic, or other focusing units to produce efficient counting which could allow for dramatically increased volumetric and sample throughput. The design lends itself to simple, cost effective POC applications.
Radionuclide therapy combines the specificity of cancer cell targeting with the well-known cytotoxic effects of ionizing radiation. It has the potential to simultaneously eradicate both the primary tumor site as well as the metastatic extent of the disease throughout the body. Despite these advantages, a major limitation of radionuclide therapy is that conventional internal dosimetry methods have failed to accurately predict treatment outcomes in patients. It has been proposed that this shortcoming is due to the fact that therapeutic dosimetry is typically performed using methods developed for diagnostic imaging that do not consider patient specific anatomy, pharmacokinetics and dose heterogeneity. Even worse, dosimetry methods are often completely neglected during treatment planning and dose prescription is based on Body Surface Area/Weight (BSA/BW) and platelet counts. The most widely used internal dosimetry software is OLINDA/EXM (Organ Level Internal Dose Assessment for Exponential Modeling). OLINDA/EXM calculates absorbed doses using the formulism developed by the Medical Internal Radiation Dose (MIRD) committee of the Society of Nuclear Medicine. In its current implementation, OLINDA/EXM requires the use of S factors, which are derived in phantoms that do not accurately represent the actual patient. In addition, OLINDA/EXM assumes that the activity distribution in each organ and tumor is uniform, when it has been well established that activity distributions are often highly heterogeneous. Furthermore, the tumor is not modeled in the phantom, which introduces additional uncertainties in the calculation. Due to these limitations, there is a clear need for improved dosimetry methods for therapeutic nuclear medicine that utilize patient-specific imaging information. A variety of efforts are currently underway to develop such methods. There are two commercial software packages on the market for voxel level internal dosimetry: STRATOS (Philips Technologies, Aachen, Germany) and VoxelDose (DosiSoft, Cachan, France), but both packages are not widely used in clinical practice. In addition, several academic platforms have been developed. Similar to these offerings, RAPID is a robust, fast and multimodal platform. However, RAPID has been designed with tremendous flexibility to perform accurate dosimetry for any radiopharmaceutical in the pre-clinical and clinical setting in order to aid in drug development and to assist in radionuclide therapy treatment planning. This presentation will overview existing dosimetry platforms while focusing on the Monte Carlo-based platform developed at the University of Wisconsin known as RAPID (Radiopharmaceutical Assessment Platform for Internal Dosimetry). RAPID is a robust and universal platform that can provide accurate 3D internal dosimetry for any radiopharmaceutical during both the pre-clinical and clinical stages of drug development as well as clinical treatment planning in TRT. (authors)
A reliable linear system is presented for grid generation in 2-D and 3-D. The method is robust in the sense that convergence is guaranteed but is not as reliable as other nonlinear elliptic methods in generating nonfolding grids. The construction of nonfolding grids depends on having reasonable approximations of cell aspect ratios and an appropriate distribution of grid points on the boundary of the region. Some guidelines are included on approximating the aspect ratios, but little help is offered on setting up the boundary grid other than to say that in 2-D the boundary correspondence should be close to that generated by a conformal mapping. It is assumed that the functions which control the grid distribution depend only on the computational variables and not on the physical variables. Whether this is actually the case depends on how the grid is constructed. In a dynamic adaptive procedure where the grid is constructed in the process of solving a fluid flow problem, the grid is usually updated at fixed iteration counts using the current value of the control function. Since the control function is not being updated during the iteration of the grid equations, the grid construction is a linear procedure. However, in the case of a static adaptive procedure where a trial solution is computed and used to construct an adaptive grid, the control functions may be recomputed at every step of the grid iteration.