Search NASASearch

SEARCH · Search NASA

Results for “Coarse-graining”

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 37 records · Page 2

Computational Insights into the Salt-Induced Modulation of Electron Transporting Conjugated Polyelectrolytes

The morphological and electronic properties of conjugated polyelectrolytes (CPEs) are highly sensitive to their ionic environment and remain poorly understood. To elucidate structure–property relationships in CPEs, we investigate the role of salt concentration on CPE morphology and hole conductivity using a quantum mechanically informed coarse-grained (CG) model coupled with semiclassical rate theory. Under good solvent conditions, high salt concentration induces torsional disorder along the conjugated backbone, decreasing hole delocalization. In contrast, under poor solvent conditions, high salt concentration promotes CPE aggregation, leading to thicker fibers and increased hole mobilities. Collectively, this work characterizes the competing interactions governing CPE assembly and hole transport as a function of salt concentration, highlighting ion engineering as a powerful strategy for tailoring the properties of mixed-conducting polymers.

diseases

Development of an Amine Oxide Polyzwitterion Brush Martini Model with Polarizable Water and Ions

Abstract In order to accurately simulate the fouling process of proteins onto polyzwitterion brushes, models that accurately capture the hydration properties and chain conformations of such brushes must first be established. We developed a Martini coarse-grained (CG) model for amine oxide polyzwitterion (PNOMA) brushes, a promising class of antifouling materials, in polarizable water and ions by fitting to all-atom bond and angle distributions, monomer hydration free energy, monomer–monomer distance potential of mean force (PMF), and monomer–salt radial distribution functions (RDFs). Martini 2.2P was selected for compatibility with the established polarizable water and ion models. For comparison with PNOMA, we also constructed models for conventional sulfobetaine (PSBMA) and phosphorylcholine (PMPC) polyzwitterions and the polycation PMETAC using established nonbonded bead types from the literature and refitting bond and angle potentials. We simulated each polymer brush chemistry for varying grafting density and chain length, validating brush height scaling relations against experimental data. The CG models captured the relative hydration strengths among different polyzwitterion chemistries, and brush heights extrapolated to higher molecular weights are in agreement with experimental ellipsometry data. We find that chain swelling of the superhydrophilic PNOMA brushes lies between that of the traditional polyzwitterions PSBMA/PMPC and the polycation PMETAC. For PNOMA brushes in NaCl solution, simulated brush height decreases with salt concentration due to the selectively strong interactions between amine oxide and sodium ions.

Walker, Christopher C. [Oak Ridge National Laborat

Molecular Mechanisms Underlying Surfactant-Based Plastics De-Inking

Surfactant-mediated binder removal is critical for de-inking pretreatment in plastics recycling, yet the molecular mechanisms governing surfactant performance remain poorly understood. We used all-atom (AA) and coarse-grained (CG) molecular dynamics (MD) simulations alongside alkaline surfactant washing experiments to investigate interactions between a series of surfactants and a polyether urethane (PEU) binder in solution and on a polyethylene (PE) surface. Experiments reveal a range of de-inking efficiencies, ranging from <25 to >95% depending upon surfactant headgroup charge and tail length. AA simulations reveal that charged surfactants reach stable levels of surfactant coverage, while nonionic surfactants aggregate on the binder. CG umbrella sampling calculations quantify the thermodynamics of binder desorption in water. In ∼0.25 M surfactant solutions, up to a 52% reduction in the free energy barrier is computed, with trends in good agreement (R 2 = 0.92, Pearson’s r = –0.96, Spearman’s ρ = –0.83) with experimental de-inking efficiencies. We find that charged surfactants are more effective than nonionic surfactants for de-inking and propose three regimes of surfactant de-inking processes: good de-inking occurs in surfactants that promote PEU desorption with a low radius of gyration (R g ); moderate de-inking occurs when surfactants stabilize PEU but increase R g ; poor de-inking occurs in surfactants that aggregate on the binder and promote extensive anchoring to the surface. Overall, these molecular-level insights have the potential to guide the design of surfactant formulations for plastics recycling applications.

ink removal

Influence of Rigidity–Hydration Coupling on Size-Dependent Diffusion in Hydrated Polymer Membranes

Selective ion transport in polymer membranes depends critically on how penetrant motion couples to polymer dynamics and hydration. Yet, the mechanistic interplay between polymer rigidity, water content, and penetrant size remains poorly understood, especially in the regime where the penetrant diameter, polymer Kuhn length, and correlation length are comparable. Here, we employ coarse-grained molecular dynamics simulations to systematically investigate penetrant diffusion in hydrated polymer networks across a broad range of water volume fractions, chain rigidities, and penetrant sizes. The results reveal a transition from a decoupled regime, where small penetrants diffuse nearly independently of polymer relaxation, to a coupled regime in which large penetrants require cooperative polymer motion for transport. Increasing polymer rigidity amplifies the sensitivity of diffusivity to hydration, particularly at low water content, leading to pronounced deviations from Stokes−Einstein scaling. Comparison with scaling theories and free-volume models shows that classical nanoparticle-based frameworks fail to capture this intermediate regime. To address this gap, we extend the Yasuda model to incorporate polymer rigidity through a single parameter that quantifies the dynamic contribution of chain stiffness to free-volume fluctuations. The resulting model collapses diffusivity data across all sizes, water contents, and rigidities, providing a unified description of penetrant transport in hydrated polymer matrices. Furthermore, these findings establish polymer rigidity as a key, tunable determinant of diffusion and offer a framework for interpreting size-dependent transport in ion-selective membranes.

diffusion

A Robotic High-Throughput Grid-Search Platform for Mapping Phase Behavior in Triblock Copolymer–Homopolymer Blends

We present a high-throughput experimental investigation of the phase behavior in triblock copolymers (PS-b-PB-b-PS and PS-b-PI-b-PS) and polystyrene (PS) homopolymer blends as a function of homopolymer molecular weight (MW) and blend ratio. Using a robotic thin-film processing platform (NOVA) integrated with Grazing Incidence Small-Angle X-ray Scattering (GISAXS) and Atomic Force Microscopy (AFM), we systematically mapped the order–disorder transition (ODT) boundaries and domain spacing evolution across a broad MW range (4.0–101.3 kDa) with varying homopolymer loadings (10% to 90%). The results reveal three distinct regimes: low-MW homopolymers, corresponding to the wet-brush regime produced only gradual domain swelling before disordering at high blend ratios (weight fraction); medium-MW homopolymers, corresponding to thedry-brush regime induced significant domain spacing increase up to 80% followed by earlier disordering, while high-MW homopolymers led to macrophase separation with minimal changes in domain spacing. Additionally, coarse-grained molecular dynamics simulations confirmed our experimental finding that in the low-MW region, the PS homopolymer uniformly distributed in the PS domain. These findings demonstrate that homopolymer molecular weight critically governs both the extent of domain swelling and the onset of disorder in triblock copolymer systems. This high-throughput platform enables the rapid mapping of composition–morphology relationships and can be integrated with AI/ML tools for designing next-generation nanostructured polymers.

36 MATERIALS SCIENCE

Insights into the Structure of Ultrasmall Fluorescent Core–Shell Silica Nanoparticles

Ultrasmall fluorescent core–shell nanoparticles (NPs) with a silica core and poly(ethylene glycol) ligand shell are the earliest example of hybrid NPs that have received U.S. investigational new drug FDA approval. They are among only a few inorganic NPs translated to safety, diagnostic, and therapeutic human clinical trials. Despite these achievements, little is known about the exact structure of their 3–4 nm sized silica cores. We report the surprising discovery of a well-defined pentagonal bipyramidal core structure preferentially formed in the aqueous synthesis built from seven primary silica NPs. A combination of reverse-phase high-performance liquid chromatography, cryogenic transmission electron microscopy, and coarse-grained simulations provides fundamental insights into this magic-size cluster formation and its unusual stability. Here, results rationalize the successful NP synthesis scale-up from 1 mL to 50 L, provide clues to the recent discovery of their self-therapeutic properties in oncology via ferroptosis, an iron-dependent cell death mechanism, and promise improved control of particle size distribution via chromatographic separations.

cluster chemistry

Structural Heterogeneity and Hydrodynamics of an Intrinsically Disordered Protein Condensate

Biology demonstrates precise control over the free-energy landscape through the selective partitioning of biomacromolecules into membraneless organelles, enabling essential functions such as biochemical transformations, signaling cascades, and mechanical reinforcement. Although the function of these condensates depends on their underlying structure and hydrodynamics, molecular-scale information on these systems remains sparse. Here, in this study, neutron scattering is used to probe the organization and dynamics of the intrinsically disordered N-terminal domain of Galectin-3, an extracellular lectin responsible for facilitating liquid–liquid phase separation on the cellular surface, in both dilute and condensed phases. Dilute solutions contain isolated protein chains in equilibrium with mesoscopic clusters, whereas the condensed phase adopts a bicontinuous, microemulsion-like morphology. The dilute phase behavior is quantitatively described by coarse-grained polymer models from soft-matter physics, demonstrating their predictive power for complex biological proteins. At elevated concentrations, the proteins self-assemble akin to block copolymers, microphase separating through the aggregation of hydrophobic domains along the protein contour. The resulting condensate remains fluid-like despite a 25-fold increase in concentration; its internal hydrodynamics slow by only a factor of 3 relative to dilute protein chains. These results provide a molecular-level framework for how disordered proteins achieve both the structural complexity and dynamic fluidity of biomolecular condensates.

Carrick, Brian R. [Massachusetts Inst. of Technolo

Why Seeding Works When Nucleation Barriers Vanish

Crystallization is a process governed by the interplay between nucleation and growth. While crystalline seeds are known to reduce nucleation barriers and accelerate crystallization under nucleation-limited conditions, their influence when nucleation is not the limiting step remains poorly understood. This creates a mechanistic puzzle in systems where nucleation barriers are already negligible, yet seeding still accelerates crystallization. Here, the synthesis of zeolites is a quintessential example of growth-limited crystallization in which seeds accelerate the process despite negligible homogeneous nucleation barriers. Using coarse-grained molecular dynamics simulations─validated across two zeolites and the unrelated case of ice crystallization─we establish that growth-limited crystallization produces many small, misoriented crystallites whose slow coarsening into larger domains controls the emergence of X-ray-detectable crystallinity. Local crystalline order, structural coherence, and X-ray detectability are therefore kinetically decoupled milestones: the first can be reached rapidly while the latter two lag significantly. Seeds resolve this lag by imposing a common orientational registry on nascent crystallites, enabling their coherent coalescence into large seed-bound domains. This early coherence-building step produces a crystallite-size asymmetry that accelerates subsequent coarsening, advancing the onset of X-ray-detectable crystallinity without necessarily increasing the nucleated fraction. We conclude that under growth-limited conditions, the apparent induction period observed in powder X-ray diffraction reflects the time required to build long-range coherence, not the time to form crystalline material. Accordingly, seeds function not primarily by reducing nucleation barriers but by enforcing spatial coherence, thereby shortening the time required to develop long-range order detectable by X-ray diffraction.

Growth-Limited Crystallization

Interfacial Inversion of Stealth Surfactants

Amphiphilic macromolecular surfactants segregate to liquid–liquid interfaces, thereby reducing the interfacial tension and free energy. Here, we investigated “stealth surfactants” in the form of core–shell bottlebrush polymers comprised of pH-responsive diblock copolymer side chains forming a hydrophilic core and a hydrophobic shell, enabling solubility in oil. At liquid–liquid interfaces, these polymers undergo a structural “inversion”, with hydrophilic blocks segregating into the aqueous phase and hydrophobic blocks residing in the oil phase. The reconfiguration kinetics and surfactant properties are influenced by multiple factors, including the molecular weights of the backbone and side chain components, the hydrophilic-to-hydrophobic balance of the side chains, and the pH of the aqueous phase. An observed nonmonotonic dependence of interfacial tension with time is attributed to a progressive structural inversion, where the projected area of the macromolecule onto the interface decreases. To validate this inversion hypothesis, interfacial properties were characterized by sum-frequency generation vibrational spectroscopy, which revealed configurational changes of the core–shell bottlebrush polymers at the fluid interface and revealed a pH-dependent interfacial coverage. Coarse-grained molecular dynamics simulations supported these experimental findings, showing that the pH-responsive core and hydrophobic shell assume a time-averaged configuration with orientations parallel and perpendicular to the plane of the interface, respectively. These findings open routes to design multistimuli-responsive polymeric surfactants and compatibilizers, expanding their potential applications in advanced interfacial systems.

Stealth surfactants

Investigating PVC polymer–plasticizer interactions with atomistic MD simulations and potential of mean force calculations

For this work, atomistic molecular dynamics (MD) simulations coupled with potential of mean force (PMF) calculations were employed to investigate the interactions between PVC polymer chains containing 6–20 repeating units and various plasticizers, with the goal of identifying potential replacements for the toxic plasticizer di(2-ethylhexyl)phthalate (DEHP) used in blood bags. The selected plasticizers belong to various chemical families, such as orthophthalates, citrates, adipates, and the terephthalate DEHT. Both the polymer and the plasticizers lack ionizable groups, and their interactions are primarily governed by van der Waals and electrostatic forces. A model correlating PMF profiles with interaction forces was developed and validated across polyvinyl chloride (PVC) polymers of different lengths and all investigated plasticizers. This model provides insight into how structural variations in plasticizers influence their respective PMF values. The study ranks the investigated plasticizers based on binding affinity, identifying TOTM (tris(2-ethylhexyl) trimellitate, TEHTM), BTHC (butyryl trihexyl citrate, Citroflex B-6), ATHC (acetyl trihexyl citrate, Citroflex A-6, CA-6), and DEHT (bis(2-ethylhexyl)terephthalate, DOTP/DEHTP) as promising alternatives for further investigation. This comprehensive study encompasses PVC polymers of four different lengths and 14 plasticizers, with all data averaged over 30 independent simulations. The approach provides a deeper understanding of molecular interactions, enabling the tailoring of polymer–plasticizer systems for diverse applications, including extractables and leachables, with relevance spanning materials science to biomedical engineering, and also serves as a basis for developing coarse-grained simulation protocols. Overall, this study provides valuable insights for designing safer and more efficient plasticizer substitutes and is well supported by other studies.

Shet, Sai Athmeeya G. [Sri Sathya Sai Institute of

Molecular dynamics based study on the effects of cation size on the local structure and diffusion in polymerized ionic liquids

We have used coarse-grained molecular dynamics simulations to understand the effect of cation size on ion diffusion in polymerized ionic liquids at temperatures well above the glass transition temperature of the polymers. We investigated dependencies of the diffusion constant on the cation radius and static dielectric constant and interpreted these results in terms of underlying structural changes, and decoupling phenomenon. We have found non-monotonic effects of the cation radius on the diffusion constant, with a maximum at an intermediate radius resulting from two different size dependent effects. Changes in the radial distribution functions characterizing spatial distribution of cations with respect to anions result in a monotonic increase in the coordination number. This yields an increase in the effective interaction energy as a function of cation radius, which can explain the decrease in the diffusion constant. This increase is counteracted at small cation radii by a corresponding increase in decoupling of cation and anion dynamics. In addition, the diffusion constant is found to increase on increasing the static dielectric constant for all radii of the cations. These results highlight the significant role of decoupling in designing polymerized ionic liquids with an enhanced diffusion constant of cations.

Gillespie, Colin [ORNL] (ORCID:0000000320287080)

High-throughput small-angle X-ray scattering reveals effective structure factor transitions linked to high-concentration antibody viscosity

High-concentration monoclonal antibody (mAb) formulations are often constrained by elevated viscosity, largely driven by protein–protein interactions, which complicates manufacturing and limits subcutaneous delivery. Early viscosity risk assessment is essential during discovery, yet traditional measurements require large sample volumes, and lack high-throughput capability. Here, we develop a high-throughput small-angle X-ray scattering (SAXS) protocol to detect mAb self-association at dilute concentrations, enabling early predictive insights into high-concentration viscosity. Synchrotron SAXS measurements were conducted for 21 mAbs formulated in a histidine buffer at pH 6.0. An initial subset of 10 mAbs analyzed across 1–150 mg/mL revealed that effective structure factor transitions in the low-q region, indicative of interparticle interactions, consistently emerged below 25 mg/mL. Subsequently, 11 additional mAbs were analyzed at 1–25 mg/mL using automated liquid handling and flow cells to enable high-throughput screening. High-viscosity mAbs exhibited detectable low-q upturns at concentrations ≤10 mg/mL, whereas low-viscosity mAbs showed downturns. A classification criterion based on effective structure factor transitions accurately classified all high- and low-viscosity mAbs at 150 mg/mL, offering a scalable, sample-efficient alternative to conventional methods. These results extend recent findings on the concentration-dependent sensitivity of SAXS to short-range attractions, demonstrating that they can emerge at lower concentrations than previously reported. This study presents the most comprehensive and diverse SAXS dataset for mAbs reported to date within a single formulation, providing a valuable resource for developing and validating coarse-grained models that can more accurately capture intermolecular interactions governing high-concentration solution behavior, thereby enabling rational antibody engineering and improved developability.

36 MATERIALS SCIENCE

Analog and symbolic computation through the Koopman framework

We develop a Koopman operator framework for studying the computational structure of dynamical systems. Specifically, we show that the resolvent of the Koopman operator provides a natural abstraction of halting, yielding a ‘Koopman halting problem’ that is recursively enumerable in general. For symbolic systems, such as those defined on Cantor space, this operator formulation captures reachability between clopen sets, while for equicontinuous systems we prove that the Koopman halting problem is decidable. Our framework demonstrates that absorbing (halting) states in coarse-grained finite automata correspond to Koopman eigenfunctions with eigenvalue one, while cycles in the transition graph impose spectral constraints associated with periodic dynamics. These results provide a unifying perspective on computation in symbolic and analog systems, showing how computational universality is reflected in operator spectra, invariant subspaces, and algebraic structures. Beyond symbolic dynamics, this operator-theoretic lens opens pathways to analyze the computational properties of a broader class of dynamical systems, including polynomial and analog models, and suggests that computational hardness may admit dynamical signatures in terms of Koopman spectral structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC

Seamlessly joining length scales: From atomistic thermal graphs to anisotropic continuum conductivity

Thermal transport in complex solids is governed by local structure, defects, and anisotropy, yet most continuum models still rely on oversimplified and homogenized conductivities. Here, we bridge atomistic and continuum descriptions by building finite element (FE) models directly from the site-projected thermal conductivity (SPTC), an atomic-level decomposition of the Green–Kubo thermal conductivity. We introduce a toolkit, the “Simulator Collection for Atomic-to-Continuum Scales (SCACS)”, which uses a graph neural network to predict SPTC on large atomic structures, coarse-grains these fields into anisotropic conductivity tensors, and embeds them into the heat-flow FE equation with a customized, anisotropy-aware adaptive mesh refinement scheme. Applied to silicon nanostructures, the resulting FE models act as representative volume elements, reproduce bulk conductivities, and capture interfacial and defect-driven anisotropy while maintaining thermodynamic consistency. Additionally, SCACS predicts experimental conductance trends and fields. This work demonstrates a general route for transferring atomistic transport information into device-scale thermal simulations with physics-based approximations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Complex orders and chirality in the classical Kitaev-Γ model

It is well recognized that the low-energy physics of many Kitaev materials is governed by two dominant energy scales, the Ising-type Kitaev coupling 𝐾 and the symmetric off-diagonal Γ coupling. An understanding of the interplay between these two scales is therefore the natural starting point toward a quantitative description that includes subdominant perturbations that are inevitably present in real materials. This study focuses on the classical 𝐾−Γ model on the honeycomb lattice, with a specific emphasis on the region 𝐾< 0 and Γ > 0 , which is the most relevant for the available materials and which remains enigmatic in both quantum and classical limits, despite much effort. We employ large-scale Monte Carlo simulations on specially designed finite-size clusters and unravel the presence of a complex multisublattice magnetic order in a wide region of the phase diagram, whose structure is characterized in detail. We show that this order can be quantified in terms of a coarse-grained scalar-chirality order, featuring a counterrotating modulation on the two spin sublattices. Here, we also provide a comparison to previous studies and discuss the impact of quantum fluctuations on the phase diagram.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

iMMD

Iterative Multiscale Molecular Dynamics in OpenMM: Accelerating Conformational Sampling of Biomolecular Systems by Iterating All-Atom and Coarse-Grained Molecular Dynamics Simulations

Do, Hung Nguyen [Los Alamos National Laboratory]

Why Seeding Works When Nucleation Barriers Vanish

This paper explains why seeds can speed up crystallization even when the formation of new crystal nuclei is not the slow step. In zeolite synthesis, local crystalline order can form relatively early, but the material may still appear poorly crystalline by X-ray diffraction because many small crystallites are misoriented and lack long-range coherence. Using coarse-grained molecular simulations, the study showed that seeds help solve this coherence problem. Rather than primarily lowering a nucleation barrier, seeds impose a common orientation on nearby growing crystallites, allowing them to merge into larger, better aligned domains. This accelerates the development of X-ray-detectable crystallinity even if the total amount of locally ordered material has not changed much. The work changes how induction times in zeolite synthesis should be interpreted: the delay before X-ray crystallinity appears can reflect the time needed to build long-range structural coherence, not simply the time needed to form crystalline material. This provides a new explanation for why seeding is effective in growth-limited crystallization and offers design principles for controlling crystallization rate and crystal quality.

36 MATERIALS SCIENCE

High-Order Methods in NASA’s Next Generation of Computational Fluid Dynamics Tools

The missions of the National Aeronautics and Space Administration (NASA) routinely produce unique requirements and challenges for development and application of Computational Fluid Dynamics (CFD) methods. NASA presently embodies four distinct Mission Directorates: Aeronautics Research, Exploration Systems, Science, and Space Operations. These missions generate requirements for systems that operate in a wide variety of environments. They range from the high-speed flight of aerodynamically optimized vehicles operating in the earth’s atmosphere to spacecraft designed for missions that don’t favor aerodynamic optimization, some operating in the atmosphere of planets and planetary moons such as Mars and Venus or Saturn’s moon Titan. Systems supporting these vehicles, such as rocket and jet propulsion, reaction control systems, fluid and thermal transfer systems, etc. can also generate their own unique set of flow phenomena that challenge today’s CFD methodology. Through the NASA Engineering and Safety Center (NESC), NASA annually conducts state-of-the-discipline assessments in fifteen distinct engineering disciplines. These assessments are performed by the NASA Technical Fellows that lead Technical Discipline Teams (TDT) of recognized experts in these fifteen areas. In the Aerosciences discipline, three topics have been identified as the top challenges for the discipline: aero-plume interaction prediction, unsteady separated flows, and aerothermodynamic prediction. These challenge areas are defined by the Agency’s high-risk projects and problems on which the NESC is requested to perform independent testing, analysis, and assessments. When viewed as a whole, these tests, analyses, and assessments provide a clear view of the recurring technical challenges facing Agency engineers and researchers and can be used to guide future research and technology development. The present state-of-the-art in the application of CFD at NASA is the use of Reynolds-Averaged Navier- Stokes (RANS) solvers, primarily executed in a steady-state mode of operation. In isolated cases, Unsteady RANS (URANS) solvers have been employed when steady RANS solutions produce poorly converging or oscillating results or in cases, such as aeroelastic analysis, which require unsteady aerodynamic simulation. For most traditional external and internal aerodynamic flows, structured overset grids or unstructured grids are employed to minimize geometric modeling and grid generation times. Grid adaptation, primarily as a series of coarse-grain intermediate processing steps is also seeing use on particularly complex flow problems and configurations. In the case of aerothermodynamic flows, engineers have been forced to continue to employ structured grid techniques as the present unstructured grid methodology has proven inadequate in the prediction of surface heating. In the area of aero-plume interaction modeling, two-gas, frozen chemistry simulation is generally the state-of-the- art, with some production solvers capable of predicting flows with only a single gas component. Prediction of flows falling into the afore-mentioned top Aerosciences technical challenges have severely stressed the present state-of-the-art in CFD prediction and for some problems, such as unsteady separated flows and aero-plume interaction cases, engineers have begun employing Large Eddy Simulation (LES) and Hybrid RANS/LES techniques. In some isolated aero-propulsion interaction cases, chemically reacting flow simulations have been applied. These methods are highly evolutionary and engineers have little experience in their application, so they cannot be heavily relied upon in today’s application environment. Therefore, this leads one to muse over which numerical technologies will be included in the CFD tools that will be employed 30 years in the future. This presentation will describe specific technical problems that have stressed NASA’s traditional CFD methods to their breaking point and will link these issues to the Agency’s top Aerosciences technical challenges. The discussion will then shift to the characteristics of future CFD solvers that will be required to attack these challenges and how these characteristics differ from the present state-of-the- art. High-order methods certainly appear to have a place in the development of future CFD tools and some of the physical characteristics of our most challenging problems suggest that high-order methods are the only way to effectively solve them. But there are some relatively severe implementation issues that face these methods, particularly in the area of general applicability and robust operation as an engineering tool. Desired characteristics of next-generation CFD solvers will be discussed and the author’s view of which emerging numerical technologies might be employed to address these attributes will also be presented

David M Schuster