Search NASA⌕ Search

SEARCH · Search NASA

Results for “distributed predictive control”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Incremental planning to control a blackboard-based problem solver

To control problem solving activity, a planner must resolve uncertainty about which specific long-term goals (solutions) to pursue and about which sequences of actions will best achieve those goals. A planner is described that abstracts the problem solving state to recognize possible competing and compatible solutions and to roughly predict the importance and expense of developing these solutions. With this information, the planner plans sequences of problem solving activities that most efficiently resolve its uncertainty about which of the possible solutions to work toward. The planner only details actions for the near future because the results of these actions will influence how (and whether) a plan should be pursued. As problem solving proceeds, the planner adds new details to the plan incrementally, and monitors and repairs the plan to insure it achieves its goals whenever possible. Through experiments, researchers illustrate how these new mechanisms significantly improve problem solving decisions and reduce overall computation. They briefly discuss current research directions, including how these mechanisms can improve a problem solver's real-time response and can enhance cooperation in a distributed problem solving network.

Durfee, E. H.↗
DOE OSTI2026

Multiscale modeling of packed-bed microwave reactors and estimation of intrinsic materials' permittivity

Modeling of packed-bed microwave reactors relies on an accurate representation of particle size, shape, and distribution within the bed, as well as the particles' dielectric properties. The measured permittivity of microwave susceptors (powders or structured materials) depends on the geometric features of the particles and the porosity of the bed, as well as the specific form factor of a structured material. These are effective properties and cannot be used to analyze other reactor configurations unless the geometric effects are removed. Therefore, we introduce a methodology for extracting the intrinsic particle permittivity from experimentally measured effective permittivity by combining cavity-based measurements with multiscale simulations and machine learning. Further, we develop the first multiscale model of packed-bed microwave reactors that incorporate particle effects (geometric features, random packing, and particle contact). This approach bridges macroscopic observables with mesoscopic physics, enabling analysis of local hotspots, arcing, and contact effects that control reactor performance. Using polymer-based spherical activated carbon (PBSAC) and silicon carbide (SiC) as examples, we demonstrate that the inferred particle permittivity is consistent with independent experimental heating profiles we collect from microwave reactors without adjustable parameters. Finally, this methodology establishes a foundation for predictive, multiscale design of microwave packed-bed reactors that explicitly accounts for particle-scale effects, enabling the estimation of intrinsic permittivity for the first time.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A Generation-Storage Coordination Dispatch Strategy for Power System Based on Causal Reinforcement Learning

In the backdrop of global energy transformation, power systems integrating high proportions of renewable energy sources are facing unprecedented challenges in operational stability and dispatch efficiency. To address these challenges, this study introduces a generation-storage coordination real-time dispatch strategy based on Causal Power System Dynamic Reinforcement Learning (CPSDRL). Diverging from traditional reinforcement learning approaches, CPSDRL innovatively incorporates causal inference within the state prediction model - the crux of model-based reinforcement learning - thereby establishing the Power Causal Dynamic Model (PCDM). Assisted by the prior knowledge of power systems, the model significantly enhances prediction accuracy and reliability through a two-stage training process. Utilizing PCDM, this study further applies a direct policy search algorithm to optimize the real-time dispatch strategy. Experimental results indicate that the proposed method improves the stability of generation-storage coordination real-time dispatch and exhibits competitive advantages in sample efficiency and computational speed, compared to traditional model-based and model-free reinforcement learning algorithms. This method is expected to enhance the practicality and adaptability of causal reinforcement learning techniques in power system scheduling and control.

causal reinforcement learning↗
NASA NTRS2019

Bowles-Tatnall Wake Vortex Encounter Hazard Metric

Wake vortex spacing standards constrict the terminal area throughput and impose severe constraints on the overall capacity and efficiency of the National Airspace System. For more than two decades starting in the early 1990s, the National Aeronautics and Space Administration conducted extensive research on characterizing the formation and evolution of aircraft wakes. This multidisciplinary work included comprehensive field experiments (Pruis et al. 2016), flight tests (Vicroy et al. 1998), and wind tunnel tests (Rossow 1994; Chow et al. 1997). Parametric studies using large eddy simulations (Proctor 1998; Proctor et al. 2006) were conducted in order to develop fast-time models for the prediction of wake transport and decay (Ahmad et al. 2016). Substantial effort was spent on the formulation of acceptable vortex hazard metrics (Tatnall 1995; Hinton and Tatnall 1997). Several wake encounter severity metrics have been suggested in the past, which include the wake circulation strength, vortex-induced rolling moment coefficient (Clv), bank angle, and the roll control ratio (Tatnall 1995; Hinton and Tatnall 1997; Van der Geest 2012). The vortex-induced rolling moment coefficient introduced by Bowles and Tatnall (Tatnall 1995; Gloudemans et al. 2016) has been used extensively for risk and safety analysis of newly proposed air traffic management concepts and procedures. The original method of Bowles and Tatnall assumed a constant wing loading (the wing lift-curve slope, CL is constant), which resulted in an overestimation of the vortexinduced rolling moment coefficient. Bowles (2014) suggested a correction to the original method that provides more accurate values of Clv and which is also consistent with the underlying physics of the problem. The overestimation of Clv in the original method can be corrected by assuming an elliptical lift distribution. Figure 1.1 illustrates the correction in Clv achieved by the modified method.

Joel Malissa↗
DOE OSTI2022

Measurements of the groomed and ungroomed jet angularities in pp collisions at $ \sqrt{s}$ = 5.02 TeV

The jet angularities are a class of jet substructure observables which characterize the angular and momentum distribution of particles within jets. These observables are sensitive to momentum scales ranging from perturbative hard scatterings to nonperturbative fragmentation into final-state hadrons. We report measurements of several groomed and ungroomed jet angularities in pp collisions at $\sqrt{s}$ = 5.02 TeV with the ALICE detector. Jets are reconstructed using charged particle tracks at midrapidity (|η| < 0.9). The anti-$k_T$ algorithm is used with jet resolution parameters R = 0.2 and R = 0.4 for several transverse momentum $p^{ch jet}_{T}$ intervals in the 20–100 GeV/c range. Using the jet grooming algorithm Soft Drop, the sensitivity to softer, wide-angle processes, as well as the underlying event, can be reduced in a way which is well-controlled in theoretical calculations. We report the ungroomed jet angularities, λ α , and groomed jet angularities, λ α,g , to investigate the interplay between perturbative and nonperturbative effects at low jet momenta. Various angular exponent parameters α = 1, 1.5, 2, and 3 are used to systematically vary the sensitivity of the observable to collinear and soft radiation. Results are compared to analytical predictions at next-to-leading-logarithmic accuracy, which provide a generally good description of the data in the perturbative regime but exhibit discrepancies in the nonperturbative regime. Moreover, these measurements serve as a baseline for future ones in heavy-ion collisions by providing new insight into the interplay between perturbative and nonperturbative effects in the angular and momentum substructure of jets. They supply crucial guidance on the selection of jet resolution parameter, jet transverse momentum, and angular scaling variable for jet quenching studies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2007

NASA Tech Briefs, May 2007

Topics include: Noise-Canceling Helmet Audio System; Program Analyzes Spacecraft/Ground Radio Links; Two-Way Communication Using RFID Equipment and Techniques; Six-Message Electromechanical Display System; Scanning Terahertz Heterodyne Imaging Systems; Master Clock and Time-Signal-Distribution System; Synchronous Phase-Resolving Flash Range Imaging; Integrated Radial Probe Transition From MMIC to Waveguide; Bar-Code System for a Microbiological Laboratory; MMIC Amplifier Produces Gain of 10 dB at 235 GHz; Mapping Nearby Terrain in 3D by Use of a Grid of Laser Spots; Digital Beam Deflectors Based Partly on Liquid Crystals; Narrow-Band WGM Optical Filters With Tunable FSRs; Better Finite-Element Analysis of Composite Shell Structures; Computing Spacecraft-Pointing Vectors for Limb Tracking; Enhanced Master Controller Unit Tester; Rover Graphical Simulator; Increasing Durability of Flame-Sprayed Strain Gauges; Multifunctional, High-Temperature Nanocomposites; Multilayer Impregnated Fibrous Thermal Insulation Tiles; Radiation-Shielding Polymer/Soil Composites; Film/Adhesive Processing Module for Fiber-Placement Processing of Composites; Fabrication of Submillimeter Axisymmetric Optical Components; Electrochemical Disposal of Hydrazines in Water; Statistical Model of Evaporating Multicomponent Fuel Drops; Resistively Heated SiC Nozzle for Generating Molecular Beams; Compact Packaging of Photonic Millimeter-Wave Receiver; Diffractive Combiner of Single-Mode Pump Laser-Diode Beams; Wide-Band, High-Quantum-Efficiency Photodetector; A Robustly Stabilizing Model Predictive Control Algorithm; Modeling Evaporation of Drops of Different Kerosenes; Development of Vapor-Phase Catalytic Ammonia Removal System; Several Developments in Space Tethers; Design Concept for a Nuclear Reactor-Powered Mars Rover; Formation-Initialization Algorithm for N Spacecraft; and DNSs of Multicomponent Gaseous and Drop-Laden Mixing Layers Achieving Transition to Turbulence.

Source record↗
NASA NTRS2016

Air Traffic Controller Acceptability of Unmanned Aircraft System Detect-and-Avoid Thresholds

A human-in-the-loop experiment was conducted with 15 retired air traffic controllers to investigate two research questions: (a) what procedures are appropriate for the use of unmanned aircraft system (UAS) detect-and-avoid systems, and (b) how long in advance of a predicted close encounter should pilots request or execute a separation maneuver. The controller participants managed a busy Oakland air route traffic control sector with mixed commercial/general aviation and manned/UAS traffic, providing separation services, miles-in-trail restrictions and issuing traffic advisories. Controllers filled out post-scenario and post-simulation questionnaires, and metrics were collected on the acceptability of procedural options and temporal thresholds. The states of aircraft were also recorded when controllers issued traffic advisories. Subjective feedback indicated a strong preference for pilots to request maneuvers to remain well clear from intruder aircraft rather than deviate from their IFR clearance. Controllers also reported that maneuvering at 120 seconds until closest point of approach (CPA) was too early; maneuvers executed with less than 90 seconds until CPA were more acceptable. The magnitudes of the requested maneuvers were frequently judged to be too large, indicating a possible discrepancy between the quantitative UAS well clear standard and the one employed subjectively by manned pilots. The ranges between pairs of aircraft and the times to CPA at which traffic advisories were issued were used to construct empirical probability distributions of those metrics. Given these distributions, we propose that UAS pilots wait until an intruder aircraft is approximately 80 seconds to CPA or 6 nmi away before requesting a maneuver, and maneuver immediately if the intruder is within 60 seconds and 4 nmi. These thresholds should make the use of UAS detect and avoid systems compatible with current airspace procedures and controller expectations.

unmanned aircraft systems↗
DOE OSTI2019

Effects of Fault Roughness on Coseismic Slip and Earthquake Locations

We report that fault zone structure is well known to exert strong controls on earthquake properties including coseismic slip distribution, rupture propagation direction, and hypocenter location. It has also been well established that the principal slip surface, which accommodates the majority of earthquake displacement, exhibits roughness at all scales following self-affine fractal distributions. Here we explore the relationship between fault roughness and specific earthquake properties including coseismic slip distribution and hypocenter location based on long-term simulations of earthquake catalogs on fractally rough faults. We begin by using the von Kármán autocorrelation function to procedurally generate single faults with different fractal roughness properties, which we place in a homogeneous elastic solid and apply pure right-lateral shear at a constant back slip rate with the earthquake simulator RSQSim. Running the simulations for 10,000 years each, we generate millions of earthquakes including thousands of events with Mw > 6.0, which rupture the surface. We show that the patterns of surface rupture in these large events follow self-affine fractal distributions with consistent fractal dimension related distinctly to the fractal dimension of the fault. In addition, the hypocenters of these large events occur in very specific predictable locations where the longest wavelength structure produces a stress asperity (i.e., restraining bend). The resulting patterns can explain many features observed on real fault systems, including clustered hypocenter locations, spatially variable coseismic slip distributions, and characteristic slip recurrence behavior. Lastly, these results demonstrate a quantitative link between a directly measurable fault property—roughness—and the properties of future earthquakes.

58 GEOSCIENCES↗
NASA NTRSDate not supplied

NASA Langley FUN3D Analyses in Support of the 1st AIAA Stability and Control Prediction Workshop

This work summarizes the results of FUN3D analyses conducted for the 1st AIAA Stability and Control Workshop on behalf of participants from the NASA Langley Research Center. The workshop was created to establish best practices for the prediction of stability and control derivatives using computational fluid dynamics and assess the limitations of these methods when those best practices are applied. The inaugural workshop considered the ONERA version of the NASA/Boeing Common Research Model, which includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. Wind tunnel data at small sideslip angles remain unpublished and served as ‘blind’ data for computational comparisons. The present research generated workshop test case data using the NASA FUN3D code, which is a parallelized, unstructured, node-based, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. Steady- state numerical simulations were conducted for workshop test cases investigating the following: grid convergence, Mach number effect on static stability, wind tunnel sting increments, static stability-derivative calculations, and a sideslip angle sweep. Results were generated for two series of unstructured, mixed-element grids, one set provided by the workshop and another set created using the HeldenMesh grid generation software. The results provided include total- and component-level breakdowns of the force and moment coefficients, in addition to sectional pressure distributions for the wing and tail components for comparisons to wind tunnel data.

CFD↗
DOE OSTI2022

Predictive Large-Eddy Simulation of Supercritical-Pressure Reactive Flows in the Cold Ignition Regime

This report describes a one-year study performed under DOE sponsorship, continuing the investigation of high-pressure turbulent reactive flows. The interest is in the effect of the chemical species distribution in high-pressure turbulent flows in the presence of strong temperature gradients as would occur during reactions in realistic flows. The prime example of such flows are boundary layers in which the wall is at a lower temperature than that of the fluid, as would be the case in Diesel engines. Previous DOE BES work in the program further highlighted the importance of the boundary layer as a configuration for fundamental studies: soot formation in boundary layers is still a problem poorly understood and depends on the availability of particular chemical species at that location, and when boundary layer Large Eddy Simulations results were compared with measurements, the agreement was unfavorable, showing that this important ‘unit’ problem is not well understood. To understand this unit problem, modeling and Direct Numerical Simulations of this unit problem were conducted for the simplest possible multispecies system, that is a binary-species system. The results discovered a new phenomenon, that is, Soret effect induced uphill diffusion. The far-reaching implication is that through an imposed wall/fluid temperature difference it is possible to control the distribution of the species in the boundary layer. Because the results have been documented in one paper published in the refereed literature, and also in conference papers, this final report is in the form of an Executive Summary succinctly describing the results and putting them in perspective with respect to existing information. The refereed and conference papers published are individually listed as Appendices and attached to this report. One manuscript is still in progress and is thus not listed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Numerical study of groundwater flow cycling controlled by seawater/freshwater interaction in Woodville Karst Plain

Seawater intrusion due to sea level rise and climate change could significantly contaminate coastal groundwater resources, particularly in Florida, the flat low-land state in the United States. Based on the field investigation and hydrological measurements, a three-dimensional SEAWAT model is developed to evaluate the groundwater flow cycling and seawater intrusion to freshwater system in the Woodville Karst Plain (WKP), a typical karst groundwater system in the Floridan aquifer. The karst conduit network in the aquifer acts as fast flow pathway for groundwater flow and solute transport, so seawater could deeply intrude into the aquifer. Wakulla Spring, an inland spring 17 km from the coast and a coastal submarine spring, Spring Creek Spring Complex are connected through the conduit network. The flow direction between the two springs switches under various rainfall conditions in this region, thus the discharges at two karst springs are used to estimate the location of seawater/freshwater mixing interface. The SEAWAT modeling results indicate that the mixing interface, defined as 2 PSU (Practical Salinity Unit), intrudes 3 to 5 km through the subsurface karst conduit during the dry season and severely contaminates nearly 1 km width of groundwater around the conduit. The salinity distribution and the distance of seawater intrusion through the conduit system are very sensitive to precipitation variation and the sea level boundary condition. Furthermore, predictions are made for seawater intrusion to the aquifer under various sea level rise, precipitation scenarios and water pumping. The results show that the seawater intrusion could reach and contaminate inland freshwater systems if sea level rises 1.0 m or during a long-term no-precipitation season. Finally, this study provides insights for modeling and predicting the vulnerability of a coastal karst aquifer through the simulation of variable-density flow.

58 GEOSCIENCES↗
DOE OSTI2023

Modelling galaxy cluster triaxiality in stacked cluster weak lensing analyses

Counts of galaxy clusters offer a high-precision probe of cosmology, but control of systematic errors will determine the accuracy of this measurement. Using Buzzard simulations, we quantify one such systematic, the triaxiality distribution of clusters identified with the redMaPPer optical cluster finding algorithm, which was used in the Dark Energy Survey Year-1 (DES Y1) cluster cosmology analysis. We test whether redMaPPer selection biases the clusters’ shape and orientation and find that it only biases orientation, preferentially selecting clusters with their major axes oriented along the line of sight. Modelling the richness–mass relation as log-linear, we find that the log-richness amplitude ln (A) is boosted from the lowest to highest orientation bin with a significance of 14σ, while the orientation dependence of the richness-mass slope and intrinsic scatter is minimal. We also find that the weak lensing shear-profile ratios of cluster-associated dark haloes in different orientation bins resemble a ‘bottleneck’ shape that can be quantified with a Cauchy function. We test the correlation of orientation with two other leading systematics in cluster cosmology – miscentering and projection – and find a null correlation. The resulting mass bias predicted from our templates confirms the DES Y1 finding that triaxiality is a leading source of bias in cluster cosmology. However, the richness-dependence of the bias confirms that triaxiality does not fully resolve the tension at low-richness between DES Y1 cluster cosmology and other probes. Our model can be used for quantifying the impact of triaxiality bias on cosmological constraints for upcoming weak lensing surveys of galaxy clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

High-Resolution Comonomer Sequencing of Blocky Brominated Syndiotactic Polystyrene Copolymers Using 13C NMR Spectroscopy and Computer Simulations

This work demonstrates the first high-resolution comonomer sequencing of Blocky brominated syndiotactic polystyrene (sPS-co-sPS-Br) copolymers based on pentad assign- ments of the quaternary carbon region of the nuclear magnetic resonance spectrum. Copolymers containing p-bromostyrene (Br-Sty) units were prepared in matched sets using postpolymerization bromination methods carried out in the heterogeneous gel state (Blocky) and homogeneous solution state (Random). Quantitative information from the quaternary carbon spectra, heteronuclear multiple bond correlation spectroscopy, electronic structure calculations, and simulated statistically random copolymers was correlated to confirm the carbon resonance assignments for all 20 possible pentad comonomer sequences. Using the experimental pentad sequence prevalences, a computer code was developed to simulate chains with microstructures typical of each sample as a means to visually represent the copolymer blockiness with quantitative precision. Based on the microstructure and distribution of run lengths in these chains, the simulations revealed that the Blocky copolymers contain a high degree of blockiness. By comparing the run lengths in the simulated chains to the average number of styrene units in a crystalline segment of sPS (found by small-angle X-ray scattering), copolymer crystallizability was predicted. For the simulated Blocky B-21% (21 mol % Br-Sty) chain, the probability of randomly selecting a styrene unit in a crystallizable block was 25.8%, while that in the simulated Random R-18% was zero, in excellent agreement with the experimental crystallization behavior measured by differential scanning calorimetry. Additionally, these predictions confirmed that the simulated chains accurately represent the ensemble of chains in their respective copolymer samples. Furthermore, each simulated Blocky chain contained one or more long sPS blocks that paralleled the measured 38-40 styrene units spanning a crystalline segment within the sPS/CCl4 gel. This finding affirmed that the long sPS segments originated from the precise lamellar structure within the heterogeneous gel morphology (i.e., block length is correlated with lamellar thickness). Overall, the ability to tailor the copolymer microstructure through control of the semicrystalline gel morphology opens the door to synthesizing ordered copolymers by postpolymerization functionalization processes with unprecedented levels of compositional control.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Assimilation of Reconstructed Radiances from IASI Principal Component Scores into the GEOS-ADAS

Hyperspectral Infrared sounders such as IASI, AIRS, and CrIS have long been an integral part of radiance assimilation in numerical weather prediction (NWP), providing vertical profiles of water vapor and temperature information. Principal Component Scores (PCS) are a lossy form of compression that retains most information, such as temperature and moisture, by using a large training set of atmospheric profiles. However, PCS may not well represent profiles which are rare events, such as volcanic eruptions, and drops some sources of random noise. There has been an increased interest in the use of PCS as EUMETSAT plans to distribute future geostationary sounder radiances from MTG-IRS via PCS only. NWP centers use two approaches to deal with PCS: direct assimilation of the PCS by modifying the radiative transfer model to produce PCS and the associated Jacobians, or a simpler approach of decompressing the PCS and reconstructing the radiances back into channel space to allow assimilating radiances without modifications to the data assimilation system. EUMETSAT has developed a PCS product for IASI that has been operational since 2011. We utilize this product opting for the simpler approach, decompressing IASI PCS into channel space, and assimilating those radiances using the GEOS-ADAS. We then compare this with a control using the standard IASI radiance product. Resulting differences in global forecast statistics, differences in Forecast Sensitivity to Observation Impact, along with implications for implementation and quality control are discussed.

Bryan M. Karpowicz↗
DOE OSTI2019

Chance Constraints for Improving the Security of AC Optimal Power Flow

This paper introduces a scalable method for improving the solutions of ac optimal power flow (AC OPF) with respect to deviations in predicted power injections from wind and other uncertain generation resources. The aim of this paper is on providing solutions that are more robust to short-term deviations, and that optimize both the initial operating point and a parametrized response policy for control during fluctuations. We formulate this as a chance-constrained optimization problem. To obtain a tractable representation of the chance constraints, we introduce a number of modeling assumptions and leverage recent theoretical results to reformulate the problem as a convex, second-order cone program, which is efficiently solvable even for large instances. Our experiments demonstrate that the proposed procedure improves the feasibility and cost performance of the OPF solution, while the additional computation time is on the same magnitude as a single deterministic AC OPF calculation.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Mercury Speciation of Hanford 241-AP-107 Tank Waste Samples - 20384

The concentration of key mercury species in the Hanford tank waste has not been previously measured. Based on process knowledge, an estimated inventory of 2,100 kg of mercury is assumed and was distributed in varied amounts across 149 single-shell tanks (SSTs) and 28 double-shell tanks (DSTs). Limited analyses of individual SSTs and DSTs for total mercury were all below detection, yielding upper limits in concentration estimates based solely on instrument detection limits, resulting in an uncertain mercury inventory for Hanford tanks. Accurate and reliable chemical speciation and quantification is desired to anticipate, predict, and abate the mass movements of various mercury species through the Hanford waste processing flowsheet. Highly sensitive separation and quantification methods for total, elemental, and monomethyl mercury species, in tank waste, have been developed at the Pacific Northwest National Laboratory's (PNNL) Radiochemical Processing Laboratory (RPL). A non-radiological environmental method to quantify mercury species has been adapted for application to Hanford tank 241-AP-107 raw supernatant feed, filtered feed, and cesium-decontaminated samples. These samples represent feed inventory and pretreatment process effluents prior to immobilization by vitrification. This method separates and preconcentrates total, elemental, and monomethyl mercury onto selective solid substrates with little-to-no radioactive background. The solid samples can then be released from strict radiological control, allowing for analytical work to be conducted in low-level radiological facilities. Separation of mercury species was demonstrated at PNNL's RPL in Richland, Washington before transport and analytical quantification using cold vapor atomic fluorescence spectrometry (CVAFS) at PNNL's Marine Science Laboratory (MSL) in Sequim, Washington. Additional studies were conducted to assess the stability of monomethyl mercury in distilled matrices and preconcentrated onto the Carbotrap{sup R} solid substrate to assess sample degradation over time during transport. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2002

High-Efficiency High-Resolution Global Model Developments at the NASA Goddard Data Assimilation Office

The Data Assimilation Office (DAO) has been developing a new generation of ultra-high resolution General Circulation Model (GCM) that is suitable for 4-D data assimilation, numerical weather predictions, and climate simulations. These three applications have conflicting requirements. For 4-D data assimilation and weather predictions, it is highly desirable to run the model at the highest possible spatial resolution (e.g., 55 km or finer) so as to be able to resolve and predict socially and economically important weather phenomena such as tropical cyclones, hurricanes, and severe winter storms. For climate change applications, the model simulations need to be carried out for decades, if not centuries. To reduce uncertainty in climate change assessments, the next generation model would also need to be run at a fine enough spatial resolution that can at least marginally simulate the effects of intense tropical cyclones. Scientific problems (e.g., parameterization of subgrid scale moist processes) aside, all three areas of application require the model's computational performance to be dramatically improved as compared to the previous generation. In this talk, I will present the current and future developments of the "finite-volume dynamical core" at the Data Assimilation Office. This dynamical core applies modem monotonicity preserving algorithms and is genuinely conservative by construction, not by an ad hoc fixer. The "discretization" of the conservation laws is purely local, which is clearly advantageous for resolving sharp gradient flow features. In addition, the local nature of the finite-volume discretization also has a significant advantage on distributed memory parallel computers. Together with a unique vertically Lagrangian control volume discretization that essentially reduces the dimension of the computational problem from three to two, the finite-volume dynamical core is very efficient, particularly at high resolutions. I will also present the computational design of the dynamical core using a hybrid distributed-shared memory programming paradigm that is portable to virtually any of today's high-end parallel super-computing clusters.

Lin, Shian-Jiann↗
NASA NTRS2002

High-Efficiency High-Resolution Global Model Developments at the NASA Goddard Data Assimilation Office

The Data Assimilation Office (DAO) has been developing a new generation of ultra-high resolution General Circulation Model (GCM) that is suitable for 4-D data assimilation, numerical weather predictions, and climate simulations. These three applications have conflicting requirements. For 4-D data assimilation and weather predictions, it is highly desirable to run the model at the highest possible spatial resolution (e.g., 55 kin or finer) so as to be able to resolve and predict socially and economically important weather phenomena such as tropical cyclones, hurricanes, and severe winter storms. For climate change applications, the model simulations need to be carried out for decades, if not centuries. To reduce uncertainty in climate change assessments, the next generation model would also need to be run at a fine enough spatial resolution that can at least marginally simulate the effects of intense tropical cyclones. Scientific problems (e.g., parameterization of subgrid scale moist processes) aside, all three areas of application require the model's computational performance to be dramatically improved as compared to the previous generation. In this talk, I will present the current and future developments of the "finite-volume dynamical core" at the Data Assimilation Office. This dynamical core applies modem monotonicity preserving algorithms and is genuinely conservative by construction, not by an ad hoc fixer. The "discretization" of the conservation laws is purely local, which is clearly advantageous for resolving sharp gradient flow features. In addition, the local nature of the finite-volume discretization also has a significant advantage on distributed memory parallel computers. Together with a unique vertically Lagrangian control volume discretization that essentially reduces the dimension of the computational problem from three to two, the finite-volume dynamical core is very efficient, particularly at high resolutions. I will also present the computational design of the dynamical core using a hybrid distributed- shared memory programming paradigm that is portable to virtually any of today's high-end parallel super-computing clusters.

Lin, Shian-Jiann↗