SEARCH · Engineering Papers

Results for “Refining”

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.

SourceResearch connectionFrom yearThrough year

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

At least 811 records · Page 45

All sources NASA NTRS DOE OSTI

NASA NTRS1989

An adaptive embedded mesh procedure for leading-edge vortex flows

A procedure for solving the conical Euler equations on an adaptively refined mesh is presented, along with a method for determining which cells to refine. The solution procedure is a central-difference cell-vertex scheme. The adaptation procedure is made up of a parameter on which the refinement decision is based, and a method for choosing a threshold value of the parameter. The refinement parameter is a measure of mesh-convergence, constructed by comparison of locally coarse- and fine-grid solutions. The threshold for the refinement parameter is based on the curvature of the curve relating the number of cells flagged for refinement to the value of the refinement threshold. Results for three test cases are presented. The test problem is that of a delta wing at angle of attack in a supersonic free-stream. The resulting vortices and shocks are captured efficiently by the adaptive code.

Powell, Kenneth G.↗

NASA NTRS1998

Accurate Modeling of the Postbuckling Response of Composite Panels with Skewed Stiffeners

Modeling procedures which predict accurately the postbuckling response of composite panels with skewed stiffeners are presented. A particular experimental program is considered and detailed finite element models are created to represent the conditions in the experiment. First, a simple baseline model with minimal detail is defined. Then, a number of modeling refinements that improve the accuracy of the model are identified. These refinements are applied individually to the baseline model to assess the effect of each refinement on the predicted response. The individual refinements are then combined to create an advanced finite element model. In addition to model refinements which affect the overall panel response, a modeling procedure for simulating nonuniform initial load introduction that is specific to an particular experiment is described and demonstrated. By implementing a number of modeling refinements, a finite element model is developed which improves the correlation between the measured and predicted responses. The skewed stiffener detail and anisotropic skin in the selected problem create axial-shear stiffness coupling which makes this problem sensitive to the constraint on the shear displacement. The remainder of the modeling refinements considered are applicable to general structural configurations. The results presented may be used to guide future efforts to improve correlation between predicted and measured responses for post buckled stiffened panels.

Young, Richard D.↗

NASA NTRS2016

Uncertainty Reduction using Bayesian Inference and Sensitivity Analysis: A Sequential Approach to the NASA Langley Uncertainty Quantification Challenge

This paper presents a computational framework for uncertainty characterization and propagation, and sensitivity analysis under the presence of aleatory and epistemic un- certainty, and develops a rigorous methodology for efficient refinement of epistemic un- certainty by identifying important epistemic variables that significantly affect the overall performance of an engineering system. The proposed methodology is illustrated using the NASA Langley Uncertainty Quantification Challenge (NASA-LUQC) problem that deals with uncertainty analysis of a generic transport model (GTM). First, Bayesian inference is used to infer subsystem-level epistemic quantities using the subsystem-level model and corresponding data. Second, tools of variance-based global sensitivity analysis are used to identify four important epistemic variables (this limitation specified in the NASA-LUQC is reflective of practical engineering situations where not all epistemic variables can be refined due to time/budget constraints) that significantly affect system-level performance. The most significant contribution of this paper is the development of the sequential refine- ment methodology, where epistemic variables for refinement are not identified all-at-once. Instead, only one variable is first identified, and then, Bayesian inference and global sensi- tivity calculations are repeated to identify the next important variable. This procedure is continued until all 4 variables are identified and the refinement in the system-level perfor- mance is computed. The advantages of the proposed sequential refinement methodology over the all-at-once uncertainty refinement approach are explained, and then applied to the NASA Langley Uncertainty Quantification Challenge problem.

Uncertainty↗

DOE OSTI2023

POWTEX visits POWGEN

The high-intensity time-of-flight (TOF) neutron diffractometer POWTEX for powder and texture analysis is currently being built prior to operation in the eastern guide hall of the research reactor FRM II at Garching close to Munich, Germany. Because of the world-wide 3 He crisis in 2009, the authors promptly initiated the development of 3 He-free detector alternatives that are tailor-made for the requirements of large-area diffractometers. Herein is reported the 2017 enterprise to operate one mounting unit of the final POWTEX detector on the neutron powder diffractometer POWGEN at the Spallation Neutron Source located at Oak Ridge National Laboratory, USA. As a result, presented here are the first angular- and wavelength-dependent data from the POWTEX detector, unfortunately damaged by a 50 g shock but still operating, as well as the efforts made both to characterize the transport damage and to successfully recalibrate the voxel positions in order to yield nonetheless reliable measurements. Also described is the current data reduction process using the PowderReduceP2D algorithm implemented in Mantid [Arnold et al. (2014). Nucl. Instrum. Methods Phys. Res. A , 764 , 156–166]. The final part of the data treatment chain, namely a novel multi-dimensional refinement using a modified version of the GSAS-II software suite [Toby & Von Dreele (2013). J. Appl. Cryst. 46 , 544–549], is compared with a standard data treatment of the same event data conventionally reduced as TOF diffraction patterns and refined with the unmodified version of GSAS-II . This involves both determining the instrumental resolution parameters using POWGEN's powdered diamond standard sample and the refinement of a friendly-user sample, BaZn(NCN) 2 . Although each structural parameter on its own looks similar upon comparing the conventional (1D) and multi-dimensional (2D) treatments, also in terms of precision, a closer view shows small but possibly significant differences. For example, the somewhat suspicious proximity of the a and b lattice parameters of BaZn(NCN) 2 crystallizing in Pbca as resulting from the 1D refinement (0.008 Å) is five times less pronounced in the 2D refinement (0.038 Å). Similar features are found when comparing bond lengths and bond angles, e.g. the two N—C—N units are less differently bent in the 1D results (173 and 175°) than in the 2D results (167 and 173°). The results are of importance not only for POWTEX but also for other neutron TOF diffractometers with large-area detectors, like POWGEN at the SNS or the future DREAM beamline at the European Spallation Source.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Adaptive immersed isogeometric level-set topology optimization

Here, this paper presents for the first time an adaptive immersed approach for level-set topology optimization using higher-order truncated hierarchical B-spline discretizations for design and state variable fields. Boundaries and interfaces are represented implicitly by the iso-contour of one or multiple level-set functions. An immersed finite element method, the eXtended IsoGeometric Analysis, is used to predict the physical response. The proposed optimization framework affords different adaptively refined higher-order B-spline discretizations for individual design and state variable fields. The increased continuity of higher-order B-spline discretizations together with local refinement enables direct control over the accuracy of the representation of each field while simultaneously reducing computational cost compared to uniformly refined discretizations. A flexible mesh adaptation strategy enables local refinement based on geometric measures or physics-based error indicators. These adaptive discretization and analysis approaches are integrated into gradient-based optimization schemes, evaluating the design sensitivities using the adjoint method. Numerical studies illustrate the features of the proposed framework with static, linear elastic, multi-material, two- and three-dimensional problems. The examples provide insight into the effect of refining the design variable field on the optimization result and the convergence rate of the optimization process. Using coarse higher-order B-spline discretizations for level-set fields promotes the development of smooth designs and suppresses the emergence of small features. Moreover, adaptive mesh refinement for state variable fields results in a reduction of overall computational cost. Higher-order B-spline discretizations are especially interesting when evaluating gradients of state variable fields due to their higher inter-element continuity.

36 MATERIALS SCIENCE↗

DOE OSTI2024

The Pele Simulation Suite for Reacting Flows at Exascale

In this work, we present the Pele suite of software tools for compressible and incompressible reacting flows. The Pele suite leverages several different libraries, notably AMReX and SUNDIALS, to achieve performance portability on heterogeneous computing architectures across the supercomputing landscape. The Pele suite is comprised of PeleC, a compressible reacting flow block-structured adaptive mesh refinement solver, PeleLMeX, a low-Mach number reacting flow block-structured adaptive mesh refinement solver, Pele-Physics, a library for transport, thermodynamics, finite rate chemistry, soot, spray and radiation physics. The objective of this paper is (i) to present the code development efforts necessary to achieve highly effective and scalable applications for exascale machines and (ii) to detail the performance results of the Combustion-Pele project applications on Oak Ridge National Laboratory's Frontier. We show good weak and strong scaling results for both PeleC and PeleLMeX up to more than 50 billion cells on more than 4096 Frontier graphics processing unit nodes. We also present a capability demonstration simulation of a dual-fuel pulse compression ignition engine (six adaptive mesh refinement levels, and 60 billion cells or 2.1 trillion degrees of freedom) on Frontier, to date one of the largest simulations performed on the first exascale-class supercomputer.

adaptive mesh refinement↗

DOE OSTI2026

Hydride and Seek: Comparing Crystallographic Hydride Placement Techniques with an Open-Shell Cobalt Complex

Locating hydrides is crucial in organometallic chemistry but difficult to do accurately using X-ray diffraction. Electron diffraction has been proposed as a way to overcome this problem but has not been systematically compared to neutron diffraction and to quantum crystallography (Hirshfeld atom refinement, HAR) to test this hypothesis. Here, we present a comparative analysis of methods for a terminal cobalt hydride complex by comparing a single-crystal neutron diffraction reference structure to results from single-crystal X-ray diffraction with and without Hirshfeld atom refinement (HAR, NoSpherA2), density functional theory (DFT), and electron diffraction (3D-ED/MicroED) refined under kinematical and dynamical formalisms. Conventional X-ray diffraction gives lower precision than neutron diffraction as expected. Despite expected improvements, HAR gives systematic deviation from the neutron benchmark. Interestingly, optimized DFT equilibrium geometries are closer to the neutron value than the value from HAR. On the other hand, electron diffraction with a high-quality data set coupled with dynamical refinement localizes the hydride in difference maps and gives excellent agreement with the neutron data. Dynamical refinement is crucial, as kinematical refinement does not allow assignment of a hydride peak. This cross-modal comparison defines the conditions under which 3D-ED/MicroED delivers high-precision metal–hydride distances for this open-shell cobalt hydride.

anions↗

DOE OSTI2024

Technical and Economic Assessment and Gap Analysis of Advanced Nuclear Reactor Integration with a Reference Oil Refinery

Efforts to identify the most-economic methods to decarbonize several sectors of the U.S. economy are underway. Industrial processes such as crude-oil refining rely heavily on energy-dense and easily stored and transported fossil fuels for powering their operations. Refineries use large amounts of energy, primarily derived from fossil sources to separate crude-oil components, break down heavier hydrocarbons into lighter compounds, remove impurities, reform hydrocarbon molecules, and generate steam and electricity for pumps and compressors and other various auxiliary systems. Crude-oil refining operations such as distillation, cracking, desulfurization, reforming, utilities systems and some offsite facilities collectively account for most of the energy consumption. Other operations such as hydrocracking or hydrotreating also require hydrogen for developing hydrogenation reactions which involve substantial heating to keep the reactors at high-temperature and pressure levels. All heat and energy demands are typically provided by natural gas (NG), oil, or other fuels, which makes refinery industry one of the most-difficult sectors to decarbonize. Nuclear power is a viable and energy-dense source of clean electricity, heat, and hydrogen to provide the large, sustainable energy supply that the refining industry demands. The U.S. Department of Energy’s (DOE’s) Integrated Energy Systems (IES) program is working to perform research and development, design, economic siting, and risk analysis. This state-of-the-art work will enable the first on-site demonstrations and commercial deployments of advanced small modular nuclear reactors (SMNRs) integrated with industries such as chemical production, refining, iron and steel making, and more. IES seeks to demonstrate the ability of advanced nuclear reactors to meet the heat and power demands of these industries while reducing carbon emissions in a sustainable and cost-competitive way. The primary objective of this research effort is to analyze industrial-scale SMNR integration intended to decarbonize refining facilities. The foreseen outcome is the provision of reliable, cost-competitive, and sustainable clean energy, alongside a reduction of carbon emissions. Specifically, the focus of this work lies on meeting the reference facilities’ heat and electricity demands with nuclear power while also supplying clean hydrogen via integrated high-temperature steam electrolysis (HTSE). This report presents a comprehensive technical and economic assessment of the integration of advanced nuclear reactors into a reference refinery, leveraging financial incentives from the Inflation Reduction Act (IRA). The evaluation aims to explore the potential economic benefits and challenges associated with incorporating advanced nuclear reactors into refinery operations, particularly in terms of energy efficiency, economic implications and environmental impact. By examining both the technical feasibility and economic viability, this analysis seeks to identify existing gaps and propose solutions for successful nuclear integration implementation. The findings are intended to provide valuable insights for stakeholders considering the adoption of advanced nuclear reactors in the refining sector. A refinery reference-plant was developed, using an open-source refinery model, Petroleum Refinery Lifecycle Inventory Model (PRELIM) and expert assessment, as a base case for comparison with various nuclear integration options. The capacity of 100 kbd/day (KBD) of heavy crude-oil feed was selected to represent a general coking-type refinery with deep conversion capabilities (incorporating heavy-oil upgrading with FCC, coking, and associated hydrotreating process units), using a heavy crude-oil feed, which represents about 70% of U.S. refineries configurations. A summary of all cases considered in this study is shown in Table 1.

13 HYDRO ENERGY↗

DOE OSTI2021

Greenhouse Gas Impact of Algal Bio-Crude Production for a Range of CO2 Supply Scenarios

Refined bio-crude production from hydrothermal liquefaction of algae holds the potential to replace fossil-based conventional liquid fuels. The microalgae act as natural carbon sequestrators by consuming CO2. However, this absorbed CO2 is released to the atmosphere during the combustion of the bio-crude. Thus, the life-cycle greenhouse gas (GHG) emissions of refined bio-crude are linked to the production and supply of the materials involved and the process energy demands. One prominent raw material is CO2, which is the main source of carbon for algae and the subsequent products. The emissions associated with the supply of CO2 can have a considerable impact on the sustainability of the algae-based refined bio-crude production process. Furthermore, the diurnal algae growth cycle complicates the CO2 supply scenarios. Traditionally, studies have relied on CO2 supplied from existing power plants. However, there is potential for building natural gas or biomass-based power plants with the primary aim of supplying CO2 to the biorefinery. Alternately, a direct air capture (DAC) process can extract CO2 directly from the air. The life-cycle GHG emissions associated with the production of refined bio-crude through hydrothermal liquefaction of algae are presented in this study. Different CO2 supply scenarios, including existing fossil fuel power plants and purpose-built CO2 sources, are compared. The integration of the CO2 sources with the algal biorefinery is also presented. The CO2 supply from biomass-based power plants has the highest potential for GHG reduction, with a GHG footprint of −57 g CO2 eq./MJ refined bio-crude. The CO2 supply from the DAC process has a GHG footprint of 49 CO2 eq./MJ refined bio-crude, which is very similar to the scenario that considers the supply of CO2 from an existing conventional natural gas-based plant and takes credit for the carbon utilization.

36 MATERIALS SCIENCE↗

NASA NTRS2000

Modal Analysis and Model Correlation of the Mir Space Station

This paper will discuss on-orbit dynamic tests, modal analysis, and model refinement studies performed as part of the Mir Structural Dynamics Experiment (MiSDE). Mir is the Russian permanently manned Space Station whose construction first started in 1986. The MiSDE was sponsored by the NASA International Space Station (ISS) Phase 1 Office and was part of the Shuttle-Mir Risk Mitigation Experiment (RME). One of the main objectives for MiSDE is to demonstrate the feasibility of performing on-orbit modal testing on large space structures to extract modal parameters that will be used to correlate mathematical models. The experiment was performed over a one-year span on the Mir-alone and Mir with a Shuttle docked. A total of 45 test sessions were performed including: Shuttle and Mir thruster firings, Shuttle-Mir and Progress-Mir dockings, crew exercise and pushoffs, and ambient noise during night-to-day and day-to-night orbital transitions. Test data were recorded with a variety of existing and new instrumentation systems that included: the MiSDE Mir Auxiliary Sensor Unit (MASU), the Space Acceleration Measurement System (SAMS), the Russian Mir Structural Dynamic Measurement System (SDMS), the Mir and Shuttle Inertial Measurement Units (IMUs), and the Shuttle payload bay video cameras. Modal analysis was performed on the collected test data to extract modal parameters, i.e. frequencies, damping factors, and mode shapes. A special time-domain modal identification procedure was used on free-decay structural responses. The results from this study show that modal testing and analysis of large space structures is feasible within operational constraints. Model refinements were performed on both the Mir alone and the Shuttle-Mir mated configurations. The design sensitivity approach was used for refinement, which adjusts structural properties in order to match analytical and test modal parameters. To verify the refinement results, the analytical responses calculated using original and refined math models were compared with the measured responses. The results from this study show that the refined models predict more accurately the dynamics of the actual structure. The MiSDE test and analysis results provided the information and experience on test design, flight testing, and data analysis for the ISS and other future spacecraft, which are critical to the verification of analytical models and structural loads.

Kim, Hyoung M.↗

DOE OSTI2022

Formal verification and validation of run-to-completion style state charts using Event-B

State chart notations with ‘run to completion’ semantics are popular with engineers for designing controllers that react to environment events with a sequence of state transitions but lack formal refinement and rigorous verification methods. State chart models are typically used to design complex control systems that respond to environmental triggers with a sequential process. The model is usually constructed at a concrete level and verified and validated using animation techniques relying on human judgement. Event-B, on the other hand, is based on refinement from an initial abstraction and is designed to make formal verification by automatic theorem provers feasible. Abstraction and formal verification provide greater assurance that critical (e.g. safety or security) properties are not violated by the control system. In this paper, we introduce a notion of refinement into a ‘run to completion’ state chart modelling notation and leverage Event-B’s tool support for theorem proving. We describe the difficulties in translating ‘run to completion’ semantics into Event-B refinements and suggest a solution. We illustrate our approach and show how models can be validated at different refinement levels using our scenario checker animation tools. We show how critical invariant properties can be verified by proof despite the reactive nature of the system and how behavioural aspects of the system can be verified by testing the expected reactions using a temporal logic, model checking approach. To verify liveness, we outline a proof that the run to completion is deadlock-free and converges to complete the run.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Overall protein structure quality assessment using hydrogen-bonding parameters

Atomic model refinement at low resolution is often a challenging task. This is mostly because the experimental data are not sufficiently detailed to be described by atomic models. To make refinement practical and ensure that a refined atomic model is geometrically meaningful, additional information needs to be used such as restraints on Ramachandran plot distributions or residue side-chain rotameric states. However, using Ramachandran plots or rotameric states as refinement targets diminishes the validating power of these tools. Therefore, finding additional model-validation criteria that are not used or are difficult to use as refinement goals is desirable. Hydrogen bonds are one of the important noncovalent interactions that shape and maintain protein structure. These interactions can be characterized by a specific geometry of hydrogen donor and acceptor atoms. Systematic analysis of these geometries performed for quality-filtered high-resolution models of proteins from the Protein Data Bank shows that they have a distinct and a conserved distribution. Here, it is demonstrated how this information can be used for atomic model validation.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

ReaxFF Parameter Set for Boron Clusters and Icosahedral Boron Crystals: Comparison with Density Functional Theory and Machine-Learning Potentials

Icosahedral boron materials, which include regular icosahedra of 12 boron atoms have gained increasing attention due to their potential applications as superhard materials, semiconductors, and energy storage media. However, the synthesis of high quality crystals of these materials has been a major barrier to the development of these applications. To enable computational prediction of synthesis conditions yielding high-quality icosahedral boron crystals, herein we tested and refined a set of ReaxFF parameters for the nucleation and growth of such crystals. We focused on matching the relative energies of small boron clusters obtained by density functional theory since such small clusters and similar motifs are likely present in crystal nuclei and at the interface of growing crystals. Using a training set of B 80 clusters, including a low-energy core–shell structure containing a B 12 icosahedron core and a high-energy single-shell structure produced in preliminary ReaxFF simulations, the ReaxFF parameter set was refined to better reproduce energies calculated by density functional theory (DFT). Among existing ReaxFF parameter sets and the machine-learning interatomic potentials MACE-MP-0, MACE-MP-0b3, MACE-MPA-0, PFP v7.0.0, and SevenNet-MF-ompa, only our new parameter set and PFP v7.0.0 correctly ranked these B 80 clusters. This refinement led to improved agreement with DFT for a test set of 58 clusters consisting of 8–103 boron atoms. Furthermore, our refined parameter set yielded greater local icosahedral structure than the previously existing ReaxFF parameter set for larger scale simulations of crystallization from supercooled liquid boron. Additionally, simulations of solid boron in contact with molten nickel using our refined ReaxFF parameters yielded a boron solubility value that agrees moderately well with experimental expectations, while the previous boron parameters gave a value that was much too low.

boron↗

DOE OSTI2020

Influences of North Pacific Ocean Domain Extent on the Western U.S. Winter Hydroclimatology in Variable-Resolution CESM

Variable-resolution global climate models (VRGCMs) are a dynamical downscaling method that can reach spatiotemporal scales needed for regional climate assessments. Over the years, several users of VRGCMs have assumed where the location and extent of the refinement domain should be based on knowledge of the prevailing storm tracks and resolution dependence of important regional climate processes (e.g., atmospheric rivers [ARs] and orographic uplift), but the effect of high-resolution domain size and extent on the simulation of downstream hydroclimatic phenomena has not been systematically evaluated. Here, we use variable resolution in the Community Earth System Model (VR-CESM) to perform such a test. To do this, three VR-CESM grids were generated that span the entire, two thirds, and one third of the North Pacific and evaluated for a 30-year climatology using Atmospheric Model Intercomparison Project protocols. Simulations are compared with reanalysis products offshore (fifth-generation of the European Centre for Medium-Range Weather Forecasts atmospheric reanalysis [ERA5]) and onshore (Livneh, 2015, https://doi.org/10.1038/sdata.2015.42, and Parameter-elevation Regressions on Independent Slopes Model [PRISM]) of the western United States. The westward expansion of refinement domain influenced integrated vapor transport (IVT), which was generally high biased but minimally impacted AR characteristics. Due to slight differences in landfalling AR counts in the western United States, California winter precipitation generally improved with westward expansion of the refinement domain. Western U.S. mountain snowpack and surface temperatures were insensitive to refinement domain size and were more influenced by changes in topographic resolution and/or land surface model version. Here, given minimal dependence of simulated western U.S. hydroclimate on refinement domain size over the North Pacific, we advise future VR-CESM studies to focus grid resolution on better resolving land surface heterogeneity.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Adaptive Grid Redistribution for a 1D Model of Turbulence and Clouds

In global atmospheric models, resolving stratocumulus (Sc) in the vertical is computationally expensive. However, Sc appear only under special meteorological conditions. Therefore, there is motivation to refine the vertical grid levels adaptively. In order to facilitate the possibility of parallelization on graphical processing units, our grid adaptation method prescribes the number of vertical levels a priori. Then grid levels are relocated toward altitude ranges in need of refinement. Because the method relocates existing grid levels, rather than adding extra levels, there is a risk of creating regions with overly coarse grid spacing, that is, voids in the grid mesh. To prevent such voids from forming, a simple method is developed to impose a maximum grid spacing. To decide where to place enhanced resolution, the authors develop an empirical mesh refinement criterion. It refines grid spacing near the ground, near strong temperature gradients, and within clouds. Our grid adaptation method is implemented in a single-column model and evaluated on four test cases: decaying stratocumulus, developing shallow cumulus, a quasi-stationary stratocumulus deck, and the diurnal cycle of a dry boundary layer. In the stratocumulus cases, mesh refinement leads to improvements in both the time evolution of fields and their time averages. The other two cases show smaller differences.

Carstensen, Steffen [Univ. of Wisconsin, Milwaukee↗

DOE OSTI2020

Two-way coupled particle-turbulence interaction: Effect of numerics and resolution on fluid and particle statistics

Euler-Lagrange point-particle simulation has emerged as a premier methodology for studying dispersed particle-laden flows. This method's popularity stems from its ability to resolve fine-scale fluid structures while also tracking individual particles at reduced cost using an appropriate particle acceleration model. However, the point-particle model has known convergence issues in that refinement of the fluid grid can lead to changes in the predicted statistics. The reasons for non-convergence are twofold: the point-particle two-way coupling force in the Navier-Stokes equations requires a numerical regularization and, without careful implementation, yields a singular force on the fluid with grid refinement. The second factor that yields grid-dependent statistics is that the point-particle force model in general depends on the undisturbed fluid velocity. When the undisturbed fluid velocity is not robustly modeled in a grid-insensitive way, the calculated force for both particles and fluid will be grid-dependent, contaminating their respective statistics. While the first issue regarding regularizing the point-particle source term has received attention in the literature, the consequences of robustly modeling the undisturbed velocity in the context of grid refinement of turbulence has received little attention. Here, we consider decaying homogeneous isotropic turbulence laden with particles at different Stokes numbers. For a given Stokes number, we systematically refine the grid and demonstrate that explicitly modeling the undisturbed fluid velocity yields relative grid insensitivity for the energy of the particle and fluid phases, as well as acceleration of the particles. We also demonstrate that an appropriately defined dissipation rate is also grid-insensitive when an undisturbed fluid velocity correction is used. In contrast, when the undisturbed fluid velocity is modeled using the conventional approach of interpolating the local fluid velocity to the particle location, we show this procedure yields divergent statistics with grid refinement. In particular, we show that higher-order interpolation of the fluid velocity in two-way coupled problems is worse than lower-order interpolation, in the absence of a correction procedure to estimate the undisturbed fluid velocity. We also examine velocity derivative statistics of the fluid phase and demonstrate that these statistics are not in general convergent even when the undisturbed fluid velocity is explicitly modeled. Collectively, the observations in this paper are used to present a philosophy on the types of questions which are answerable with point-particle methods.

42 ENGINEERING↗

NASA NTRS1998

Purification and Crystal Growth of Lead Iodide by Physical Vapor Transport Method

Lead iodide (PbI2) is a layered compound semiconductor being developed as room temperature x- and gamma-ray detector. Compared to the more studied material, mercuric iodide, PbI2 has a higher melting temperature and no phase transition until liquid phase which are indications of better mechanical properties. In this study, the source material was purified by the zone-refining process, and the purest section was extracted from center of the the zone-refined ingot to be grown by physical vapor transport (PVT) method. The zone-refined material and as-grown crystals were characterized by optical microscopy and differential scanning calorimetry (DSC) to reveal the surface morphology, purity and stoichiometry. The results shows that both materials are near-stoichiometric composition, with the purity of the as-grown crystals higher than zone-refined materials. The resistivity of the as-grown crystal (10" Omega-cm) was derived from current-voltage (I-V) measurement, and is 10 times higher than the zone-refined materials. Detail results will be presented and discussed.

Wright, G. W.↗

NASA NTRS2005

Application of Parallel Adjoint-Based Error Estimation and Anisotropic Grid Adaptation for Three-Dimensional Aerospace Configurations

This paper demonstrates the extension of error estimation and adaptation methods to parallel computations enabling larger, more realistic aerospace applications and the quantification of discretization errors for complex 3-D solutions. Results were shown for an inviscid sonic-boom prediction about a double-cone configuration and a wing/body segmented leading edge (SLE) configuration where the output function of the adjoint was pressure integrated over a part of the cylinder in the near field. After multiple cycles of error estimation and surface/field adaptation, a significant improvement in the inviscid solution for the sonic boom signature of the double cone was observed. Although the double-cone adaptation was initiated from a very coarse mesh, the near-field pressure signature from the final adapted mesh compared very well with the wind-tunnel data which illustrates that the adjoint-based error estimation and adaptation process requires no a priori refinement of the mesh. Similarly, the near-field pressure signature for the SLE wing/body sonic boom configuration showed a significant improvement from the initial coarse mesh to the final adapted mesh in comparison with the wind tunnel results. Error estimation and field adaptation results were also presented for the viscous transonic drag prediction of the DLR-F6 wing/body configuration, and results were compared to a series of globally refined meshes. Two of these globally refined meshes were used as a starting point for the error estimation and field-adaptation process where the output function for the adjoint was the total drag. The field-adapted results showed an improvement in the prediction of the drag in comparison with the finest globally refined mesh and a reduction in the estimate of the remaining drag error. The adjoint-based adaptation parameter showed a need for increased resolution in the surface of the wing/body as well as a need for wake resolution downstream of the fuselage and wing trailing edge in order to achieve the requested drag tolerance. Although further adaptation was required to meet the requested tolerance, no further cycles were computed in order to avoid large discrepancies between the surface mesh spacing and the refined field spacing.

Lee-Rausch, E. M.↗