Engineering Papers⌕ Search

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.

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

At least 883 records · Page 49

Multiobjective Multidisciplinary Optimization of Low-Boom Supersonic Transports Using Multifidelity Models

A multidisciplinary optimization (MDO) method has been developed to design a computational fluid dynamics (CFD) based low-boom configuration that can be obtained from a Pareto solution of a low-fidelity multiobjective MDO problem with mission constraints. This paper refines the developed MDO method using multifidelity models for CFD-based multiobjective MDO. The refined MDO method can generate a low-boom configuration that satisfies the mission requirements, has the lowest takeoff gross weight and the longest range for the low-boom mission, trims the low-boom cruise flight with fuel redistributions, and has a reversed equivalent area distribution closely matching a low-boom target with ground noise level below 70 PLdB. The validity of the refined MDO method is demonstrated by a design study of a low-boom supersonic transport that carries 40 passengers, flies a low-boom mission with cruise Mach of 1.7 and range of 3500 nm, and cruises overwater at Mach 1.8 with range of 3882 nm. Moreover, the refined MDO method eliminates the difference between the assumed cruise weight for CFD-based low-boom inverse design optimization and the estimated cruise weight of the optimal inverse design solution with respect to the mission requirements.

multidisciplinary optimization↗

Multiobjective Multidisciplinary Optimization of Low-Boom Supersonic Transports Using Multifidelity Models

A multidisciplinary optimization (MDO) method has been developed to design a computational fluid dynamics (CFD) based low-boom configuration that can be obtained from a Pareto solution of a low-fidelity multiobjective MDO problem with mission constraints. This paper refines the developed MDO method by using multifidelity models for CFD-based multiobjective MDO and a better method for the system-level trade between the target low boom ground noise level and the overland range. The refined MDO method can generate a low-boom configuration that satisfies the mission requirements, has the lowest takeoff gross weight and the longest range for the low-boom mission, trims the low-boom cruise flight with fuel redistributions, and has a reversed equivalent area distribution closely matching a low-boom target with ground noise level below 70 PLdB. The validity of the refined MDO method is demonstrated by a design study of a low-boom supersonic transport that carries 40 passengers, flies a low-boom mission with cruise Mach of 1.7 and range of 3500 nm, and cruises overwater at Mach 1.8 with range of 3882 nm. Moreover, the refined MDO method eliminates the difference between the assumed cruise weight for CFD-based low-boom inverse design optimization and the estimated cruise weight of the optimal inverse design solution with respect to the mission requirements.

MDO↗

FUN3D Analyses of the 5th AIAA Propulsion Aerodynamics Workshop Inlet Test Case

This work presents the results of FUN3D analyses that were performed for the 5th AIAA Propulsion Aerodynamics Workshop. The workshop was separated into two sections that focused on inlet and nozzle flows. This paper focuses on the inlet section of the workshop, which considered an axisymmetric inlet in close-proximity to the ground and subjected to a crosswind. Four levels of grid refinement were provided by the workshop committee. Participants were asked to perform simulations for three crosswind velocities on a minimum of three levels of grid refinement, for a total of nine required cases. This work employed an updated set of grids provided by the committee after completion of the workshop based on feedback from participants. The updated grids featured a larger computational domain and increased refinement near the ground vortex region. This work employed the three coarsest grid levels from the updated grid set. Note that an additional extension downstream of the inlet was added to the provided grids for this work to mitigate observed convergence issues. Additionally, a limited turbulence model study was performed for the highest value of crosswind velocity and the finest grid. The results show that grid convergence was only achieved for the case with the lowest value of crosswind velocity. However, the results were generally observed to improve, relative to the experimental data, with increasing grid refinement. The FUN3D predictions were observed to agree qualitatively with the provided experimental data. However, FUN3D exhibited a tendency to underpredict the minimum value of the radially averaged total pressure at the Aerodynamic Interface Plane, which is consistent with the general findings from the workshop. Overall, FUN3D was found to perform well relative to its peers for the PAW5 workshop inlet test case.

PAI↗

FUN3D Analyses of the 5th AIAA Propulsion Aerodynamics Workshop Inlet Test Case

This work presents the results of FUN3D analyses that were performed for the 5th AIAA Propulsion Aerodynamics Workshop. The workshop was separated into two sections that focused on inlet and nozzle flows. This paper focuses on the inlet section of the workshop, which considered an axisymmetric inlet in close-proximity to the ground and subjected to a crosswind. Four levels of grid refinement were provided by the workshop committee. Participants were asked to perform simulations for three crosswind velocities on a minimum of three levels of grid refinement, for a total of nine required cases. This work employed an updated set of grids provided by the committee after completion of the workshop based on feedback from participants. The updated grids featured a larger computational domain and increased refinement near the ground vortex region. This work employed the three coarsest grid levels from the updated grid set. Note that an additional extension downstream of the inlet was added to the provided grids for this work to mitigate observed convergence issues. Additionally, a limited turbulence model study was performed for the highest value of crosswind velocity and the finest grid. The results show that grid convergence was only achieved for the case with the lowest value of crosswind velocity. However, the results were generally observed to improve, relative to the experimental data, with increasing grid refinement. The FUN3D predictions were observed to agree qualitatively with the provided experimental data. However, FUN3D exhibited a tendency to underpredict the minimum value of the radially averaged total pressure at the Aerodynamic Interface Plane, which is consistent with the general findings from the workshop. Overall, FUN3D was found to perform well relative to its peers for the PAW5 workshop inlet test case.

PAI↗

Mechanically robust high magnetic-performance Sm-Co sintered magnets through microstructure engineering

Samarium-cobalt (Sm-Co) sintered magnets have high magnetic energy densities, great resistance to demagnetization and corrosion, and excellent thermal stability in a wide temperature range (–50–550 °C). However, the utilization of these magnets is restricted by their brittleness. Based on micromechanical and the Zener pinning model, Sm-Co sintered magnets with improved mechanical properties have been designed and fabricated via microstructure engineering. A small amount of fine Sm 2 O 3 particulates (0–3 wt%) has been incorporated into Sm 2 (CoFeCuZr) 17 sintered magnets to refine the grain size by up to approximately 50% (from 45 to 22 µm) and narrow the grain size distribution. Doping with 3 wt% Sm 2 O 3 increased the flexural strength by 62% while maintaining magnetic performance. Both grain-refined unimodal microstructure and heterogeneous laminated coarse/fine grain microstructure were formed by strategically designed assemblies of Sm 2 O 3 -added Sm-Co powder feedstock mixtures. The technology is compatible with existing magnet manufacturing processes. Numerical micromechanics simulation indicates that the fracture is dominated by intragranular mode. The mechanical strength is mainly enhanced by the additive-induced grain refinement, while the small amount of Sm 2 O 3 addition has a small direct positive contribution. Finally, these magnets will be more cost-effective, efficient, and robust for various functional applications.

36 MATERIALS SCIENCE↗

Toughening Sm–Co sintered magnets via microstructure modification with additives

We report the mechanical properties of the brittle Smsingle bondCo permanent magnets are of great practical significance. However, studies on the magnets have mostly focused on their magnetic properties. This paper reports the modified microstructure and refined unimodal grain size, enhanced flexural strength, and magnetic properties of Sm 2 (Co,Fe,Cu,Zr) 17 sintered magnets doped with La 2 O 3 , MgO, and CaF 2 fine particulates. The correlations between microstructure, phase composition, and mechanical and magnetic properties were studied. Doping of a small amount (e.g., 0.5–3 wt%) of La 2 O 3 , MgO, or CaF 2 fine particulates could significantly refine the unimodal grain sizes of the Smsingle bondCo magnets via the Zener pinning effect. Moreover, doping significantly improved the flexural strengths σ of the magnets. For example, the σ values of the magnets with 0.5 wt% MgO, 1 wt% La 2 O 3 , and 1 wt% CaF 2 were approximately 65%, 63%, and 42% higher than that of the reference magnet, respectively. Micromechanical simulations revealed that the fine particles of La 2 O 3 could deflect crack growth, while the CaF 2 particles could attract or arrest cracks during the fracture process. The mechanical strengthening effect was mainly due to grain size refinement. The Smsingle bondCo magnets with 0.5–1.5 wt% CaF2 and 0.5 wt% La 2 O 3 exhibited excellent magnetic properties while doping 1–3 wt% La2O3 and 0.5–3 wt% MgO deteriorated magnetic performance. The rational design of CaF 2 - or La 2 O 3 -doped microstructure can be an economical and effective method for producing toughened Sm–Co sintered magnets with high magnetic performance.

36 MATERIALS SCIENCE↗

Effect of heat treatment on the microstructure of medium burn-up U-Mo monolithic fuel foils

Using scanning electron microscopy (SEM), this study evaluates the microstructure evolution of U-Mo fuel foils made with and without heat treatment at medium burn-up (of approximately 5 x 10 21 f/cm 3 ). The impact of annealing treatments on critical microstructural properties of the U-Mo fuel foils, including porosity, grain structure, Mo homogeneity, and fuel interaction with the Zr interlayer, was examined using large area lift outs (LALOs). The heat-treated specimens presented less grain refining at these burnups when compared to the un-heated specimens. Grain refinement was associated with porosities and fission products precipitation. Heat treatment can reduce fuel swelling during irradiation. Chemical inhomogeneity (Mo banding) was found to persist in the un-heated samples but was not present in the heat-treated samples. Thus, heat-treated U-Mo foils allows for more predictable fuel behavior under irradiation with respect to un-heated foils. The U-Mo and Zr interaction layer appears to be thicker and more continuous in the heat-treated sample which has been associated with stronger interface integrity during irradiation, as also observed in previous studies. These observations may indicate an overall improved performance of heat-treated fuel foil in a reactor. Further, the effect of local burn-up on grain size/refinement and porosities in each LALO specimen, sampled from different positions in the fuel foil, was difficult to analyze due to the large standard deviation of these parameters. Finally, evidence of grain refinement by polygonization may be present in these specimens.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Angular-spatial hp -adaptivity for radiative transfer with discontinuous Galerkin spectral element methods

Radiative transfer is important for many science and engineering applications, and numerical simulations of radiative transfer can be challenging. For instance, the radiation field is seven-dimensional – three spatial, two angular, one wavelength, and one temporal – and often features steep gradients. Therefore, memory usage is a key issue. To reduce memory, some past work has investigated the use of adaptive mesh refinement (AMR), typically for either the spatial or angular coordinate, and typically for only h -adaptivity. Here, we propose the use of AMR for the spatial and angular coordinates together, and the use of h - and p -adaptivity together as hp -AMR for the potential for further memory savings. We implemented the proposed method for several test cases in two spatial and one angular dimension, with the discontinuous Galerkin spectral element method. These test cases featured highly anisotropic angular radiation, with or without steep spatial gradients. Our primary findings from these test cases were: (1) Angular hp -adaptivity can deliver the radiation solution with the same accuracy as, and with much less computational memory than, uniform angular h - or p -refinements, or angular h -adaptivity alone. This is most obvious when the incoming radiation is highly anisotropic, in which case the savings can be orders of magnitude. (2) Full spatial-angular hp -adaptivity is more efficient in solution representation, compared to solely spatial or solely angular -adaptivity. This is most evident when steep gradients are present in both the spatial and angular distribution. These results suggest that adaptive spatial- hp angular-refinement may perform well in large-scale seven-dimensional applications.

Adaptive refinement↗

Dynamic recrystallization of a wrought magnesium alloy: Grain size and texture maps and their application for mechanical behavior predictions

Effect of dynamic recrystallization (DRX) on the grain refinement and texture modification was studied by conducting a series of hot compression on AZ31B Mg alloy. Processing-grain size-texture maps were established as a function of temperature, strain rate, strain, and the Zener-Hollomon parameter (Z). Moreover, influence of simultaneous changes in the grain size and texture on tensile yield strength and ductility of the hot-worked Mg alloy was studied using the grain size and texture maps established. The processing-grain size map showed that the DRX grain size decreased with the increase in Z, also revealing various characteristics ranging from a grain growth, grain refinement, to a bimodal distribution of ultrafine grains and partially recrystallized grains. The effect of twinning on the grain refinement was also evident at high Z conditions. The processing-texture map revealed that the initial fiber texture was altered significantly to a shear, off-normal, or extension-twin texture with the increase in Z. The Schmid-factor maps were calculated to identify a dominant deformation mode during a subsequent tensile deformation of the hot-worked samples. Finally, a corresponding Hall-Petch relationship for each dominant deformation mode was used to establish tensile yield-strength maps, which agree well with the measured data.

36 MATERIALS SCIENCE↗

pyDiSCaMB : enabling the use of multipolar scattering factors in Phenix

Multipolar scattering models, such as the transferable aspherical atom model, account for atomic chemical interactions and provide a more accurate representation of experimental data. However, the simpler independent atom model (IAM), which assumes non-interacting atoms, is the only model available in the most widely used macromolecular refinement programs. This is primarily because IAM offers a hard-to-beat combination of computational efficiency and modelling power at typical macromolecular resolutions. By contrast, more accurate multipolar modelling has historically been limited due to its computational cost and the absence of an interface between software capable of calculating structure factors and gradients based on multipolar models and software designed for macromolecular refinement. This work introduces pyDiSCaMB , a Python software package designed to integrate between the computational crystallography toolbox ( cctbx ) and the quantum crystallography library DiSCaMB ( Densities in Structural Chemistry and Molecular Biology ), thus enabling multipolar scattering models in Phenix 's toolkit. The implementation, features and capabilities of pyDiSCaMB are presented, the runtimes for the calculation of structure factor and target gradients with respect to atomic parameters are explored, and Fourier images of electrostatic potential, electron density and deformation maps are computed as illustrative examples. The pyDiSCaMB library will make multipolar modelling widely available to the structural biology community, potentially transforming refinement and model-building for both crystallography and cryogenic electron microscopy (cryoEM).

MATTS data bank↗

Transferable Hirshfeld atom model for rapid evaluation of aspherical atomic form factors

Form factors based on aspherical models of atomic electron density have brought great improvement in the accuracies of hydrogen atom parameters derived from X-ray crystal structure refinement. Today, two main groups of such models are available, the banks of transferable atomic densities parametrized using the Hansen–Coppens multipole model which allows for rapid evaluation of atomic form factors and Hirshfeld atom refinement (HAR)-related methods which are usually more accurate but also slower. In this work, a model that combines the ideas utilized in the two approaches is tested. It uses atomic electron densities based on Hirshfeld partitions of electron densities, which are precalculated and stored in a databank. This model was also applied during the refinement of the structures of five small molecules. A comparison of the resulting hydrogen atom parameters with those derived from neutron diffraction data indicates that they are more accurate than those obtained with the Hansen–Coppens based databank, and only slightly less accurate than those obtained with a version of HAR that neglects the crystal environment. The advantage of using HAR becomes more noticeable when the effects of the environment are included. To speed up calculations, atomic densities were represented by multipole expansion with spherical harmonics up to l = 7, which used numerical radial functions (a different approach to that applied in the Hansen–Coppens model). Calculations of atomic form factors for the small protein crambin (at 0.73 Å resolution) took only 68 s using 12 CPU cores.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Arginine off-kilter: guanidinium is not as planar as restraints denote

Crystallographic refinement of macromolecular structures relies on stereochemical restraints to mitigate the typically poor data-to-parameter ratio. For proteins, each amino acid has a unique set of geometry restraints which represent stereochemical information such as bond lengths, valence angles, torsion angles, dihedrals and planes. It has been shown that the geometry in refined structures can differ significantly from that present in libraries; for example, it was recently reported that the guanidinium moiety in arginine is not symmetric. In this work, the asymmetry of the N ∊ —C ζ —N η1 and N ∊ —C ζ —N η2 valence angles in the guanidinium moiety is confirmed. In addition, it was found that the C δ atom can deviate significantly (more than 20°) from the guanidinium plane. This requires the relaxation of the planar restraint for the C δ atom, as it otherwise causes the other atoms in the group to compensate by distorting the guanidinium core plane. A new set of restraints for the arginine side chain have therefore been formulated, and are available in the software package Phenix , that take into account the asymmetry of the group and the planar deviation of the C δ atom. This is an example of the need to regularly revisit the geometric restraint libraries used in macromolecular refinement so that they reflect the best knowledge of the structural chemistry of their components available at the time.

59 BASIC BIOLOGICAL SCIENCES↗

What Is the Limit of Quantification for the Minor Phase in Time-of-Flight Neutron Diffraction? A Case Study on Fe and Ni Powder Mixtures at VULCAN

A phase present in small quantities within materials may not simply serve as a secondary component; it can play a crucial role in determining the integrity, properties, and performance of the material. These minor but important phases usually draw attention in material design and processing for fundamental understanding as well as material quality control. Accurately quantifying a minor phase amid a majority phase, especially at extremely low fractions, remains a challenging task. Time-of-flight neutron diffraction, coupled with advanced pattern analysis techniques like Rietveld refinement, is a powerful tool for crystal structure identification and phase quantification. The deep penetrating capability of neutrons enables the detection and quantification of trace phases within materials. In this study, the quantification limits of time-of-flight neutron diffraction were explored using the VULCAN diffractometer at the Spallation Neutron Source, using Fe–Ni powder mixtures as a sample system. By comparing the refinement results to the known weighed values, it was determined that the reliable quantification of a minor Ni phase is achievable down to about 0.1 wt% while a Ni fraction as low as 0.02 wt% is difficult to trace. Effective control of the refinement parameters, especially the profile function parameters, are found to significantly influence the convergence of fittings and the accuracy of phase quantification.

Rietveld refinement↗

The arbitrary‐order virtual element method for linear elastodynamics models: convergence, stability and dispersion‐dissipation analysis

Abstract We design the conforming virtual element method for the numerical approximation of the two‐dimensional elastodynamics problem. We prove stability and convergence of the semidiscrete approximation and derive optimal error estimates under h ‐ and p ‐refinement in both the energy and the L 2 norms. The performance of the proposed virtual element method is assessed on a set of different computational meshes, including nonconvex cells up to order four in the h ‐refinement setting. Exponential convergence is also experimentally observed under p ‐refinement. Finally, we present a dispersion‐dissipation analysis for both the semidiscrete and fully discrete schemes, showing that polygonal meshes behave as classical simplicial/quadrilateral grids in terms of dispersion‐dissipation properties.

Antonietti, Paola F.↗

Modeling Protein–Protein and Protein–Ligand Interactions by the ClusPro Team in CASP16

ABSTRACT In the CASP16 experiment, our team employed hybrid computational strategies to predict both protein–protein and protein–ligand complex structures. For protein–protein docking, we combined physics‐based sampling—using ClusPro FFT docking and molecular dynamics—with AlphaFold (AF)‐based sampling, followed by AF‐based refinement. Our method produced numerous high‐accuracy complex models, including cases where AF alone failed, underscoring the critical role of physics‐based sampling alongside deep learning‐based refinement. For protein–ligand docking, we integrated the ClusPro LigTBM template‐based approach with a machine learning‐based confidence model for rescoring. The method preserves conserved interaction fragments derived from homologous complexes, followed by local resampling using physics‐based sampling and a diffusion model. Our template‐based strategy achieved a mean lDDT‐PLI of 0.69 across 233 targets, which was highly competitive. These results demonstrate that combining physics‐based modeling with AI‐driven refinement can significantly enhance the accuracy of both protein–protein and protein–ligand structure predictions.

Ashizawa, Ryota [Department of Applied Mathematics↗

Influence of neutron cross-section resonances on organ/tissue equivalent and effective dose coefficients for the ICRP voxel phantoms

The materials which compose the ICRP Voxel phantoms used in the computation of conversion coefficients involve neutron interaction cross-sections that have resonances at specific energies. Depending on the energy bin structure used in the computations, these cross-section resonances may occur at energies that fall between energies at which dose coefficients are computed, thus their effects may not be completely accounted for in the reported coefficients. Here In this study, a highly refined energy grid that closely follows the resonance structure in the phantom material cross-sections was identified and used to calculate dose coefficients. Both the equivalent organ/tissue doses for male and female voxel phantoms were computed as well as their summation to obtain the effective dose coefficients. The used refined energy grid tracks very closely the cross-sections in the vicinity of the resonances. The resulting refined energy grid coefficients are compared to coefficients for the coarser energy grid used in ICRP Publication 116. Additionally, reference spectra have been folded with both the fine and coarse sets of conversion coefficients. The resulting total effective doses for these reference spectra are used to assess the adequacy of the dose coefficients calculated on the original ICRP 116 energy grid. The dose coefficients were similarly computed for the local skin dose on the trunk of the body using the ICRU Report 95 phantom. The overall impact of the resonances on the organ/tissue equivalent dose, the effective dose, and the local skin dose are presented and discussed. In general, it was found that resonances can impact neutron dose coefficients, but in most cases the wide range of neutron energies encountered minimized this effect. The impact of resonances was further limited when computing effective dose due to organ/tissue summing and sex-averaging. For the neutron fields studied here, the impact was below 5%.

36 MATERIALS SCIENCE↗

Structural determination of a full-length plant cellulose synthase informed by experimental and in silico methods

Three-dimensional structure determination and prediction of proteins with intrinsically disordered regions, unstructured regions, conformational flexibility, and lacking homologous structures are challenging. We previously predicted and refined an in silico structure of a plant cellulose synthase from cotton (GhCESA1), and more recently, cryo-electron microscopy (cryo-EM) has resolved a majority of the lengths of two CESA structures from poplar (PttCESA8) and cotton (GhCESA7). However, 26–30% of these cryo-EM structures remain unresolved, including the N-terminal domain, half of the class-specific region, the gating loop region, and the C-terminal domain. Here, we describe the generation and evaluation of a full-length hybrid GhCESA1 model based on this cryo-EM PttCESA8 structure, with unresolved regions completed using this in silico refined GhCESA1 model. All-atom molecular dynamics simulations and subsequent energy minimizations were performed for the in silico and hybrid GhCESA1 models in a lipid bilayer-water-ion environment, and structural stability, dynamics, energetics, contacts, and quality were evaluated. The unresolved regions were found to be the most dynamic, in agreement with their poor electron density with cryo-EM. The hybrid model exhibited a higher total secondary structure content, more favorable intra-protein and protein-lipid interaction energies, and improved quality metrics. Moreover, hydrogen bonding was revealed to be a primary mechanism for intra-protein and protein-lipid contacts. These results demonstrate that in silico structure prediction and refinement may be useful to augment experimental structure determination, especially for disordered and unstructured regions. Furthermore, this hybrid model can serve as a steppingstone to derive full-length homology models of other CESAs found in more experimentally tractable organisms.

59 BASIC BIOLOGICAL SCIENCES↗

Impact of artificial topological changes on flow and transport through fractured media due to mesh resolution

Abstract We performed a set of numerical simulations to characterize the interplay of fracture network topology, upscaling, and mesh refinement on flow and transport properties in fractured porous media. We generated a set of generic three-dimensional discrete fracture networks at various densities, where the radii of the fractures were sampled from a truncated power-law distribution, and whose parameters were loosely based on field site characterizations. We also considered five network densities, which were defined using a dimensionless version of density based on percolation theory. Once the networks were generated, we upscaled them into a single continuum model using the upscaled discrete fracture matrix model presented by Sweeney et al. (2019). We considered steady, isothermal pressure-driven flow through each domain and then simulated conservative, decaying, and adsorbing tracers using a pulse injection into the domain. For each simulation, we calculated the effective permeability and solute breakthrough curves as quantities of interest to compare between network realizations. We found that selecting a mesh resolution such that the global topology of the upscaled mesh matches the fracture network is essential. If the upscaled mesh has a connected pathway of fracture (higher permeability) cells but the fracture network does not, then the estimates for effective permeability and solute breakthrough will be incorrect. False connections cannot be eliminated entirely, but they can be managed by choosing appropriate mesh resolution and refinement for a given network. Adopting octree meshing to obtain sufficient levels of refinement leads to fewer computational cells (up to a 90% reduction in overall cell count) when compared to using a uniform resolution grid and can result in a more accurate continuum representation of the true fracture network.

58 GEOSCIENCES↗