Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reversible computation”

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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Solubility and Complexation of the Bismuthate Ion in Nitric Acid Systems

Here, the dissolution rate and solubility of NaBiO 3 have been investigated in nitric acid systems ranging from 4–6 M HNO 3 , and were found to be 58–76 µg/cm 2 ·d and 490–830 mM, respectively. The presence of 50 mM U(VI) drastically increased the solubility to 540–1200 mM, while rates of dissolu-tion were relatively unchanged. The solubility of NaBiO 3 increased with an increase in U(VI) con-centrations at 4 M HNO 3 , with log-log analysis indicating a one-to-one complex between Bi and U and IR spectroscopic evidence monitoring uranyl stretching suggesting complex formation. Absorb-ance spectra were obtained experimentally and computationally with an absorbance band in the range of 450–600 nm that has been attributed to Bi(V). The ingrowth and decay of Bi(V) in solution was also studied as a function of mass of solid NaBiO 3 present, acidity, and temperature. The ac-tivation energies of dissolution and decomposition were calculated to be 39 ± 4 and 61 ± 6 kJ/mol, re-spectively. Furthermore, these results indicate that dissolution of NaBiO 3 into the respective Na + and BiO 3 – is occur-ring prior to undergoing reduction, a process which conventionally has been believed to occur in the reverse order.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1995

The Heat of Formation of HNO

The HNO molecule is of interest in both combustion and atmospheric chemistry. For example, Guadagnini et al. have recently presented ab initio potential energy surfaces for the three lowest lying electronic states of HNO and then used these in examining several chemical reactions that take place in the combustion of nitrogen containing fuels and in the oxidation of atmospheric nitrogen. We have previously studied the ground state potential energy surface (i.e., stationary points along the HNO reversible reaction HON path), vibrational spectrum (using an accurate quartic force field), zero-point energy, and bonding of HNO using coupled-cluster ab initio methods. HNO is also very interesting because of the unique nature of its bonding characteristics. That is, the potential energy surface is very flat along the H-N bonding coordinate thereby giving unusual harmonic and fundamental vibrational frequencies, and the H-N bond energy is rather weak in comparison to other H-N bond energies. In fact, using experimental heats of formation for HNO, H, and NO, the H- bond energy is computed to be only 49.9 kcal/ mol (298 K). However, ab initio calculations of isodesmic reaction energies involving HNO, FNO, ClNO, and several other molecules have shown that there is an inconsistency in the experimental heats of formation of the XNO (X double bond H, F, and Cl) species. Hence the motivation for this study was to determine a very accurate (Delta)H(sub f, sup o) value for HNO using state-of-the-art ab initio methods. Based on many recent studies it is evident that the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), in conjunction with large one-particle basis sets should be reliable to better than +/- 0.8 kcal/mol for this quantity. The computational methodology is described in the next section followed by our results and discussion. Conclusions are presented in the final section.

Lee, Timothy J.↗
DOE OSTI2023

Denoising diffusion models with geometry adaptation for high fidelity calorimeter simulation

Simulation is crucial for all aspects of collider data analysis, but the available computing budget in the High Luminosity LHC era will be severely constrained. Generative machine learning models may act as surrogates to replace physics-based full simulation of particle detectors, and diffusion models have recently emerged as the state of the art for other generative tasks. Here, we introduce CaloDiffusion, a denoising diffusion model trained on the public CaloChallenge datasets to generate calorimeter showers. Our algorithm employs 3D cylindrical convolutions, which take advantage of symmetries of the underlying data representation. To handle irregular detector geometries, we augment the diffusion model with a new geometry latent mapping (GLaM) layer to learn forward and reverse transformations to a regular geometry that is suitable for cylindrical convolutions. The showers generated by our approach are nearly indistinguishable from the full simulation, as measured by several different metrics.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2003

Afferent innervation patterns of the saccule in pigeons

The innervation patterns of vestibular saccular afferents were quantitatively investigated in pigeons using biotinylated dextran amine as a neural tracer and three-dimensional computer reconstruction. Type I hair cells were found throughout a large portion of the macula, with the highest density observed in the striola. Type II hair cells were located throughout the macula, with the highest density in the extrastriola. Three classes of afferent innervation patterns were observed, including calyx, dimorph, and bouton units, with 137 afferents being anatomically reconstructed and used for quantitative comparisons. Calyx afferents were located primarily in the striola, innervated a number of type I hair cells, and had small innervation areas. Most calyx afferent terminal fields were oriented parallel to the anterior-posterior axis and the morphological polarization reversal line. Dimorph afferents were located throughout the macula, contained fewer type I hair cells in a calyceal terminal than calyx afferents and had medium sized innervation areas. Bouton afferents were restricted to the extrastriola, with multi-branching fibers and large innervation areas. Most of the dimorph and bouton afferents had innervation fields that were oriented dorso-ventrally but were parallel to the neighboring reversal line. The organizational morphology of the saccule was found to be distinctly different from that of the avian utricle or lagena otolith organs and appears to represent a receptor organ undergoing evolutionary adaptation toward sensing linear motion in terrestrial and aerial species.

NASA Discipline Neuroscience↗
DOE OSTI2016

A comparative study of coarse-graining methods for polymeric fluids: Mori-Zwanzig vs. iterative Boltzmann inversion vs. stochastic parametric optimization

We construct effective coarse-grained (CG) models for polymeric fluids by employing two coarse-graining strategies. The first one is a forward-coarse-graining procedure by the Mori-Zwanzig (MZ) projection while the other one applies a reverse-coarse-graining procedure, such as the iterative Boltzmann inversion (IBI) and the stochastic parametric optimization (SPO). More specifically, we perform molecular dynamics (MD) simulations of star polymer melts to provide the atomistic fields to be coarse-grained. Each molecule of a star polymer with internal degrees of freedom is coarsened into a single CG particle and the effective interactions between CG particles can be either evaluated directly from microscopic dynamics based on the MZ formalism, or obtained by the reverse methods, i.e., IBI and SPO. The forward procedure has no free parameters to tune and recovers the MD system faithfully. For the reverse procedure, we find that the parameters in CG models cannot be selected arbitrarily. If the free parameters are properly defined, the reverse CG procedure also yields an accurate effective potential. Moreover, we explain how an aggressive coarse-graining procedure introduces the many-body effect, which makes the pairwise potential invalid for the same system at densities away from the training point. From this work, general guidelines for coarse-graining of polymeric fluids can be drawn.

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

First-principles effective Hamiltonian for finite-temperature modeling of nonperovskite ferroelectrics

First-principles-based effective Hamiltonian techniques have been widely employed for over three decades to investigate ferroelectricity and related phenomena in perovskite materials. These techniques offer high accuracy, transferability, compatibility with various finite-temperature algorithms, computational efficiency, and ease in incorporating interactions with external fields. They have been adapted to study diverse phenomena, ranging from topological dipole patterns in ferroelectric nanostructures to multicaloric effects. In this work, we develop an effective Hamiltonian for the nonperovskite ferroelectric HfO 2 (hafnia). Applying this methodology to explore the finite-temperature and finite-electric-field properties of ferroelectric hafnia revealed (1) exceptionally large intrinsic coercive fields, an order of magnitude higher than those observed in perovskite ferroelectrics; (2) their atomistic origin; and (3) the existence of a regime where the relationship between the coercive field and the energy barrier for polarization reversal is counterintuitive. Here, these developments could accelerate progress both in methodological advancements for simulating ferroics and in the atomistic understanding of a broad range of ferroelectrics.

Electric polarization↗
NASA NTRS1995

Capillary Pressure Balance with Repeated Head-Up Tilt in Humans

We measured blood density (BD). plasma density (PD) and hematocrit (Ht) during repeated 70 deg head-up tilt (HUT) consisting of a pre-drink period with two supine (P and P3) and two HUT (P2 and P4) phases of 45 min each. At the end of P4 test subjects (N=8) drank 10 ml/kg body weight of isotonic (290 mOsm/kg) sodium chloride (ISO) or hypotonic (<10 mOsm/kg) unsweetened us (HYP), or nothing (control: CON). The following periods continued the supine (P5, P7) / upright (P6) sequence. Mass density of plasma and whole blood using a "mechanical oscillation" technique were measured from ear lobe blood. Tilting altered (p<0.001) BD, PD, and Ht in all consecutive Phases. BD and PD were different (p<0.05) between CON. HYP. and ISO Phases 6 and 7. The computed mass density of fluid (FD). which moved between intra. and extravascular compartments. averaged 1008 g/l and was the same in any conditions. In the non-drink treatment (P3, P5, P7), supine PV steadly decreased when compared to P1. Tilt-induced PV shifts ranged from 9.7 to 16.7% when compared to PV du. ring the respective previous Phases. After drinking, PV increased above CON values at the end of Phase 7 by 12.9% with ISO, and by 6.6% with HYP. Progressive hemoconcentration. occurred during non-drink supine; isotonic saline ingestion increased supine PV a control level but did not stop or reverse the decrease of upright hemoconcentration. Upright density values were not slowed with repeated upright tilting; decreased fluid loss occurred in consecutive upright periods. Upright capillary pressure balance was achieved with diminished plasma volume loss during repeated tilting.

Hinghofer-Szalkay, H. G.↗
NASA NTRS2005

Matlab GUI for a Fluid Mixer

The Test and Engineering Directorate at NASA John C. Stennis Space Center developed an interest to study the modeling, evaluation, and control of a liquid hydrogen (LH2) and gas hydrogen (GH2) mixer subsystem of a ground test facility. This facility carries out comprehensive ground-based testing and certification of liquid rocket engines including the Space Shuttle Main engine. A software simulation environment developed in MATLAB/SIMULINK (M/S) will allow NASA engineers to test rocket engine systems at relatively no cost. In the progress report submitted in February 2004, we described the development of two foundation programs, a reverse look-up application using various interpolation algorithms, a variety of search and return methods, and self-checking methods to reduce the error in returned search results to increase the functionality of the program. The results showed that these efforts were successful. To transfer this technology to engineers who are not familiar with the M/S environment, a four-module GUI was implemented allowing the user to evaluate the mixer model under open-loop and closed-loop conditions. The progress report was based on an udergraduate Honors Thesis by Ms. Jamie Granger Austin in the Department of Electrical Engineering and Computer Science at Tulane University, during January-May 2003, and her continued efforts during August-December 2003. In collaboration with Dr. Hanz Richter and Dr. Fernando Figueroa we published these results in a NASA Tech Brief due to appear this year. Although the original proposal in 2003 did not address other components of the test facility, we decided in the last few months to extend our research and consider a related pressurization tank component as well. This report summarizes the results obtained towards a Graphical User Interface (GUI) for the evaluation and control of the hydrogen mixer subsystem model and for the pressurization tank each taken individually. Further research would combine the two components - mixer and tank, for a more realistic simulation tool.

Barbieri, Enrique↗
DOE OSTI2012

Subsolutions That Are Close in the Uniform Norm Are Close in the Sobolev Norm as Well

A new type weighted reverse Poincare inequality is established for a difference of two continuous weak subsolutions of a linear second order uniformly elliptic partial differential equation in the ball.This result is the key to deriving the error estimate for the gradient of the analytically unknown value function of the optimal stochastic control problem from the uniform error of the value function itself in the related numerical approximation problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

A Comparative Study of the Hypoxia PET Tracers [{sup 18}F]HX4, [{sup 18}F]FAZA, and [{sup 18}F]FMISO in a Preclinical Tumor Model

Purpose: Several individual clinical and preclinical studies have shown the possibility of evaluating tumor hypoxia by using noninvasive positron emission tomography (PET). The current study compared 3 hypoxia PET tracers frequently used in the clinic, [{sup 18}F]FMISO, [{sup 18}F]FAZA, and [{sup 18}F]HX4, in a preclinical tumor model. Tracer uptake was evaluated for the optimal time point for imaging, tumor-to-blood ratios (TBR), spatial reproducibility, and sensitivity to oxygen modification. Methods and Materials: PET/computed tomography (CT) images of rhabdomyosarcoma R1-bearing WAG/Rij rats were acquired at multiple time points post injection (p.i.) with one of the hypoxia tracers. TBR values were calculated, and reproducibility was investigated by voxel-to-voxel analysis, represented as correlation coefficients (R) or Dice similarity coefficient of the high-uptake volume. Tumor oxygen modifications were induced by exposure to either carbogen/nicotinamide treatment or 7% oxygen breathing. Results: TBR was stabilized and maximal at 2 hours p.i. for [{sup 18}F]FAZA (4.0 ± 0.5) and at 3 hours p.i. for [{sup 18}F]HX4 (7.2 ± 0.7), whereas [{sup 18}F]FMISO showed a constant increasing TBR (9.0 ± 0.8 at 6 hours p.i.). High spatial reproducibility was observed by voxel-to-voxel comparisons and Dice similarity coefficient calculations on the 30% highest uptake volume for both [{sup 18}F]FMISO (R = 0.86; Dice coefficient = 0.76) and [{sup 18}F]HX4 (R = 0.76; Dice coefficient = 0.70), whereas [{sup 18}F]FAZA was less reproducible (R = 0.52; Dice coefficient = 0.49). Modifying the hypoxic fraction resulted in enhanced mean standardized uptake values for both [{sup 18}F]HX4 and [{sup 18}F]FAZA upon 7% oxygen breathing. Only [{sup 18}F]FMISO uptake was found to be reversible upon exposure to nicotinamide and carbogen. Conclusions: This study indicates that each tracer has its own strengths and, depending on the question to be answered, a different tracer can be put forward.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2011

Assessing the toxic effects of ethylene glycol ethers using Quantitative Structure Toxicity Relationship models

Experimental determination of toxicity profiles consumes a great deal of time, money, and other resources. Consequently, businesses, societies, and regulators strive for reliable alternatives such as Quantitative Structure Toxicity Relationship (QSTR) models to fill gaps in toxicity profiles of compounds of concern to human health. The use of glycol ethers and their health effects have recently attracted the attention of international organizations such as the World Health Organization (WHO). The board members of Concise International Chemical Assessment Documents (CICAD) recently identified inadequate testing as well as gaps in toxicity profiles of ethylene glycol mono-n-alkyl ethers (EGEs). The CICAD board requested the ATSDR Computational Toxicology and Methods Development Laboratory to conduct QSTR assessments of certain specific toxicity endpoints for these chemicals. In order to evaluate the potential health effects of EGEs, CICAD proposed a critical QSTR analysis of the mutagenicity, carcinogenicity, and developmental effects of EGEs and other selected chemicals. We report here results of the application of QSTRs to assess rodent carcinogenicity, mutagenicity, and developmental toxicity of four EGEs: 2-methoxyethanol, 2-ethoxyethanol, 2-propoxyethanol, and 2-butoxyethanol and their metabolites. Neither mutagenicity nor carcinogenicity is indicated for the parent compounds, but these compounds are predicted to be developmental toxicants. The predicted toxicity effects were subjected to reverse QSTR (rQSTR) analysis to identify structural attributes that may be the main drivers of the developmental toxicity potential of these compounds.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Validation of formability of laminated sheet metal for deep drawing process using GTN damage model

In this study, we studied formability of PET/PVC laminated sheet metal which named VCM (Vinyl Coated Metal). VCM offers various patterns and good-looking metal steel used for appliances such as refrigerator and washing machine. But, this sheet has problems which are crack and peeling of film when the material is formed by deep drawing process. To predict the problems, we used finite element method and GTN (Gurson-Tvergaard-Needleman) damage model to represent damage of material. We divided the VCM into 3 layers (PET film, adhesive and steel added PVC) in finite element analysis model to express the crack and peeling phenomenon. The material properties of each layer are determined by reverse engineering based on tensile test result. Furthermore, we performed the simple rectangular deep drawing and simulated it. The simulation result shows good agreement with drawing experiment result in position, punch stroke of crack occurrence. Also, we studied the fracture mechanism of PET film on VCM by comparing the width direction strain of metal and PET film.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Inhibiting and Remodeling Toxic Amyloid-Beta Oligomer Formation Using a Computationally Designed Drug Molecule That Targets Alzheimer’s Disease

Alzheimer’s disease (AD) is rapidly reaching epidemic status among a burgeoning aging population. Much evidence suggests the toxicity of this amyloid disease is most influenced by the formation of soluble oligomeric forms of amyloid β-protein, particularly the 42-residue alloform (Aβ42). Developing potential therapeutics in a directed, streamlined approach to treating this disease is necessary. Here we utilize the joint pharmacophore space (JPS) model to design a new molecule [AC0107] incorporating structural characteristics of known Aβ inhibitors, blood-brain barrier permeability, and limited toxicity. To test the molecule’s efficacy experimentally, we employed ion mobility mass spectrometry (IM-MS) to discover [AC0107] inhibits the formation of the toxic Aβ42 dodecamer at both high (1:10) and equimolar concentrations of inhibitor. Atomic force microscopy (AFM) experiments reveal that [AC0107] prevents further aggregation of Aβ42, destabilizes preformed fibrils, and reverses Aβ42 aggregation. This trend continues for long-term interaction times of 2 days until only small aggregates remain with virtually no fibrils or higher order oligomers surviving. Pairing JPS with IM-MS and AFM presents a powerful and effective first step for AD drug development. .

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2018

Leonurine hydrochloride promotes osteogenic differentiation and increases osteoblastic bone formation in ovariectomized mice by Wnt/β-catenin pathway

Highlights: • Leonurine hydrochloride promotes osteogenic differentiation of osteoblasts in vitro. • Wnt/β-catenin pathway acts as an important role in the osteogenic differentiation affected by Leonurine hydrochloride. • Leonurine hydrochloride alleviate the effect of osteoporosis by acclerating the bone formation. • Bone mineral density and micro-structures are improved by the treatment of Leonurine hydrochloride. Leonurine hydrochloride (LH) is a synthetic chemical compound derived from leonurine that can be extracted from Leonurus sibiricus and possesses antioxidant, anti-apoptosis, and neuroprotective activities. In previous studies, LH has been demonstrated to attenuate osteoclast activity and prevent bone loss. However, it is unknown whether LH accelerates bone formation and promotes osteogenic differentiation. We systematically examined the effects of LH on ovariectomized-induced osteoporotic mice and the MC3T3-E1 osteoblastic cell line. The results revealed that LH enhanced differentiation of MC3T3-E1 cells, with a dose-dependent increase in alkaline phosphatase (ALP) activity. Moreover, LH upregulated osteogenesis-related gene expression, including osterix, alpha 1 type 1 collagen, runt-related transcription factor 2 (Runx2) and ALP, as shown by quantitative reverse transcription-polymerase chain reaction analysis. At the same time, elevated expression of low-density lipoprotein receptor-related protein 5 and β-catenin mRNA was detected in the Wnt/β-catenin pathway. A western blot analysis revealed that LH dose-dependently increased the expression of Runx2 and β-catenin, and promoted phosphorylation of glycogen synthase kinase-3β in vitro. The in vivo results showed that administering LH (15 mg/kg/d) for 8 weeks alleviated destruction of the trabecular microstructure caused by osteoporosis. LH increased the bone mineral density and trabecular number, decreased trabecular separation according to a micro-computed tomography scan. In addition, LH enhanced the expression of β-catenin and Runx2 in vivo. In conclusion, LH promoted osteogenic differentiation and bone formation in vivo and in vitro, which alleviated osteoporosis through activation of the Wnt/β-catenin pathway.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Modeling low cycle fatigue (LCF) of additively manufactured Hastelloy X using An accelerated crystal plasticity fatigue damage model

This paper presents a microstructure-based model for low cycle fatigue (LCF) behavior and life of Nickel-based alloy Hastelloy X manufactured using laser-powder bed fusion (L-PBF) additive manufacturing (AM). AM Hastelloy X, a solution-strengthened alloy, is tested at elevated temperature under fully reversed LCF conditions at different strain levels. A generalized plane strain finite element model is generated from electron backscatter diffraction (EBSD) characterization. The constitutive behavior of the material under fatigue is modeled using crystal plasticity and calibrated with both monotonic tensile and cyclic stress–strain data. The fatigue micro-crack initiation and propagation in the microstructure is modeled using a modified Chaboche fatigue damage model. An embedded boundary condition with a homogenous medium is used to apply the cyclic deformation and prevent numerically introduced over-constraints during fatigue simulation. A ‘cycle-jump’ method is used to accelerate the fatigue simulation and reduce the computational cost. The simulation results are compared to LCF experiments, showing satisfactory matches in cyclic stress behavior and number of cycles to macro-crack initiation for all applied strain ranges. In addition, the model illustrates the potential for quantifying microscale fatigue life impacting factors such as microstructure and surface roughness, which is needed to accurately quantify the reliability of AM components in service.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Neutron-Friendly Li-Ion Battery Coin Cell for In Situ 3D Visualization of Li Plating

Advanced battery characterization using in situ/operando neutron imaging is critical for uncovering degradation modes such as lithium (Li) plating in Li-ion batteries (LIBs). However, conventional LIBs hinder operando neutron radiography (NR) and in situ neutron micro-computed tomography (N-μCT) for visualizing Li plating near the graphite-separator interface due to strong attenuation from hydrogen-rich components like PP–PE–PP separators, electrolyte, and Fe-based spacers. In this work, we designed and tested a neutron-friendly battery (NFB) optimized for in situ Li detection during extreme fast charging (XFC). Guided by neutron attenuation cross-sections and material transmission, the NFB enables clear visualization at the graphite–separator interface, which is typically opaque in standard LIBs. Electrochemical tests show the NFB exhibits voltage/current responses like standard cells for up to 50 XFC cycles. However, its lower reversibility and capacity are likely due to Cu-coated Al spacer degradation from delamination or corrosion. We propose titanium spacers as a more stable alternative, albeit requiring custom machining. Using this optimized cell, we achieved simultaneous neutron tomography of multiple cells, capturing in situ 3D images of dead Li accumulation, particularly near graphite edges. These heterogeneous deposits and disconnected Li clusters suggest localized current density hotspots during XFC.

25 ENERGY STORAGE↗
DOE OSTI2024

Molecular Understanding of Nitrogen Oxide Fixation of Water-Lean Carbon Capture Solvents by Atomistic Modeling

Nitrogen oxides, present in flue gas, can cause negative impacts on amine carbon capture solvents by the formation of heat-stable salts and suspected carcinogens. Thus, to maximize the performance of water-lean solvents, a better understanding of this process in these systems is necessary. Here, a computational study for the fixation of the CO 2 capture solvent N-(2-ethoxyethyl)-3-morpholinopropan-1-amine (EEMPA) to nitramine/nitrosamine was conducted. The first step involves the dissociation of the NH bond of EEMPA, in which the homolytic mechanism is energetically more favorable than the heterolytic mechanism. The second step involves radical recombination to form N–N bonds. While NO 2 directly reacts with EEMPA, NO has almost no effect. However, in the presence of O 2 , fixation of EEMPA by NO is enhanced via the formation of N 2 O 4 species. Finally, low reaction energies indicate that the formation of nitramine/nitrosamine may be a reversible process, suggesting that EEMPA could be recovered under thermal stripping conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2010

Constraints and Opportunities in GCM Model Development

Over the past 30 years climate models have evolved from relatively simple representations of a few atmospheric processes to complex multi-disciplinary system models which incorporate physics from bottom of the ocean to the mesopause and are used for seasonal to multi-million year timescales. Computer infrastructure over that period has gone from punchcard mainframes to modern parallel clusters. Constraints of working within an ever evolving research code mean that most software changes must be incremental so as not to disrupt scientific throughput. Unfortunately, programming methodologies have generally not kept pace with these challenges, and existing implementations now present a heavy and growing burden on further model development as well as limiting flexibility and reliability. Opportunely, advances in software engineering from other disciplines (e.g. the commercial software industry) as well as new generations of powerful development tools can be incorporated by the model developers to incrementally and systematically improve underlying implementations and reverse the long term trend of increasing development overhead. However, these methodologies cannot be applied blindly, but rather must be carefully tailored to the unique characteristics of scientific software development. We will discuss the need for close integration of software engineers and climate scientists to find the optimal processes for climate modeling.

Schmidt, Gavin↗