Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geometric adaptivity”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Efficient Berry phase calculation via adaptive variational quantum computing approach

We present an adaptive variational quantum algorithm to estimate the Berry phase accumulated by a nondegenerate ground state under cyclic, adiabatic evolution of a time-dependent Hamiltonian. Our method leverages cyclic adiabatic evolution of the Hamiltonian and employs adaptive variational quantum algorithms for state preparation and evolution, optimizing circuit efficiency while maintaining high accuracy. We benchmark our approach on dimerized Fermi–Hubbard chains with four sites, demonstrating precise Berry phase simulations in both noninteracting and interacting regimes. Our results show that circuit depths reach up to 106 layers for noninteracting systems and increase to 279 layers for interacting systems due to added complexity. In addition, we demonstrate the robustness of our scheme across a wide range of parameters governing adiabatic evolution and variational algorithms. These findings highlight the potential of adaptive variational quantum algorithms for advancing quantum simulations of topological materials and computing geometric phases in strongly correlated systems.

Mootz, Martin [Ames Laboratory (AMES), Ames, IA (U↗
DOE OSTI2024

Robust Containment Queries over Collections of Rational Parametric Curves via Generalized Winding Numbers

Point containment queries for regions bound by watertight geometric surfaces, i.e., closed and without self-intersections, can be evaluated straightforwardly with a number of well-studied algorithms. When this assumption on domain geometry is not met, such methods are either unusable, or prone to misclassifications that can lead to cascading errors in downstream applications. More robust point classification schemes based on generalized winding numbers have been proposed, as they are indifferent to these imperfections. However, existing algorithms are limited to point clouds and collections of linear elements. We extend this methodology to encompass more general curved shapes with an algorithm that evaluates the winding number scalar field over unstructured collections of rational parametric curves. In particular, we evaluate the winding number for each curve independently, making the derived containment query robust to how the curves are arranged. We ensure geometric fidelity in our queries by treating each curve as equivalent to an adaptively constructed polyline that provably has the same generalized winding number at the point of interest. Our algorithm is numerically stable for points that are arbitrarily close to the model, and explicitly treats points that are coincident with curves. We demonstrate the improvements in computational performance granted by this method over conventional techniques as well as the robustness induced by its application.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Camellia: A Rapid Development Framework for Finite Element Solvers

We report the discontinuous Petrov–Galerkin (DPG) methodology of Demkowicz and Gopalakrishnan guarantees the optimality of the finite element solution in a user-controllable energy norm, and provides several features supporting adaptive schemes. The approach provides stability automatically; there is no need for carefully derived numerical fluxes (as in DG schemes) or for mesh-dependent stabilization terms (as in stabilized methods). In this paper, we focus on features of Camellia that facilitate implementation of new DPG formulations; chief among these is a rich set of features in support of symbolic manipulation, which allow, e.g., bilinear formulations in the code to appear much as they would on paper. Many of these features are general in the sense that they can also be used in the implementation of other finite element formulations. In fact, because DPG’s requirements are essentially a superset of those of other finite element methods, Camellia provides built-in support for most common methods. We believe, however, that the combination of an essentially “hands-free” finite element methodology as found in DPG with the rapid development features of Camellia are particularly winsome, so we focus on use cases in this class. In addition to the symbolic manipulation features mentioned above, Camellia offers support for one-irregular adaptive meshes in 1D, 2D, 3D, and space-time. It provides a geometric multigrid preconditioner particularly suited for DPG problems, and supports distributed parallel execution using MPI. For its load balancing and distributed data structures, Camellia relies on packages from the Trilinos project, which simplifies interfacing with other computational science packages. Camellia also allows loading of standard mesh formats through an interface with the MOAB package. Camellia includes support for static condensation to eliminate element-interior degrees of freedom locally, usually resulting in substantial reduction of the cost of the global problem. Finally, we include a discussion of the variational formulations built into Camellia, with references to those formulations in the literature, as well as an MPI performance study.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2002

Triangle Geometry Processing for Surface Modeling and Cartesian Grid Generation

Cartesian mesh generation is accomplished for component based geometries, by intersecting components subject to mesh generation to extract wetted surfaces with a geometry engine using adaptive precision arithmetic in a system which automatically breaks ties with respect to geometric degeneracies. During volume mesh generation, intersected surface triangulations are received to enable mesh generation with cell division of an initially coarse grid. The hexagonal cells are resolved, preserving the ability to directionally divide cells which are locally well aligned.

Aftosmis, Michael J.↗
DOE OSTI2015

Automatic mesh adaptivity for hybrid Monte Carlo/deterministic neutronics modeling of difficult shielding problems

The CADIS and FW-CADIS hybrid Monte Carlo/deterministic techniques dramatically increase the efficiency of neutronics modeling, but their use in the accurate design analysis of very large and geometrically complex nuclear systems has been limited by the large number of processors and memory requirements for their preliminary deterministic calculations and final Monte Carlo calculation. Three mesh adaptivity algorithms were developed to reduce the memory requirements of CADIS and FW-CADIS without sacrificing their efficiency improvement. First, a macromaterial approach enhances the fidelity of the deterministic models without changing the mesh. Second, a deterministic mesh refinement algorithm generates meshes that capture as much geometric detail as possible without exceeding a specified maximum number of mesh elements. Finally, a weight window coarsening algorithm decouples the weight window mesh and energy bins from the mesh and energy group structure of the deterministic calculations in order to remove the memory constraint of the weight window map from the deterministic mesh resolution. The three algorithms were used to enhance an FW-CADIS calculation of the prompt dose rate throughout the ITER experimental facility. Using these algorithms resulted in a 23.3% increase in the number of mesh tally elements in which the dose rates were calculated in a 10-day Monte Carlo calculation and, additionally, increased the efficiency of the Monte Carlo simulation by a factor of at least 3.4. The three algorithms enabled this difficult calculation to be accurately solved using an FW-CADIS simulation on a regular computer cluster, eliminating the need for a world-class super computer.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

Recycling facility periodic safety review: Application of a methodology used for the nuclear criticality safety review of geometrically safe or favorable equipment - 14610

In France, for all fuel cycle facilities, the reference safety frame is continuously adapted to take into account technical developments such as changes in process and/or maintenance operations. Every ten years according to French law, the facility operator is required to launch a process of Periodic Safety Review (PSR), whose conclusions are sent to the Nuclear Safety Authority (NSA). This PSR is split in two parts: first, a review of the facility to its safety documentation and authorization, second, a safety reassessment. Actually, this PSR takes account for topics (actual aging for instance) which require several years of feedback from experience to be integrated in safety evaluations, or which cannot be done on a daily basis (documentation reorganization). It also incorporates the new regulatory safety requirements. A specific methodology was developed for the first PSR of UP3-A since the publication of the law. UP3-A is a used fuel reprocessing facility of La Hague plant. This methodology, adapted for large plants with numerous process and fissile media, is based on a two-step approach intending to focus the analysis on the real hard spots. First, 'sensitive' systems are identified. Second, a specific review is performed on these systems. After a reminder of the methodology, this paper presents the results of its application on the UP3-A scope. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2013

Automatic mesh adaptivity for CADIS and FW-CADIS neutronics modeling of difficult shielding problems

The CADIS and FW-CADIS hybrid Monte Carlo/deterministic techniques dramatically increase the efficiency of neutronics modeling, but their use in the accurate design analysis of very large and geometrically complex nuclear systems has been limited by the large number of processors and memory requirements for their preliminary deterministic calculations and final Monte Carlo calculation. Three mesh adaptivity algorithms were developed to reduce the memory requirements of CADIS and FW-CADIS without sacrificing their efficiency improvement. First, a macro-material approach enhances the fidelity of the deterministic models without changing the mesh. Second, a deterministic mesh refinement algorithm generates meshes that capture as much geometric detail as possible without exceeding a specified maximum number of mesh elements. Finally, a weight window coarsening algorithm de-couples the weight window mesh and energy bins from the mesh and energy group structure of the deterministic calculations in order to remove the memory constraint of the weight window map from the deterministic mesh resolution. The three algorithms were used to enhance an FW-CADIS calculation of the prompt dose rate throughout the ITER experimental facility. Using these algorithms resulted in a 23.3% increase in the number of mesh tally elements in which the dose rates were calculated in a 10-day Monte Carlo calculation and, additionally, increased the efficiency of the Monte Carlo simulation by a factor of at least 3.4. The three algorithms enabled this difficult calculation to be accurately solved using an FW-CADIS simulation on a regular computer cluster, obviating the need for a world-class super computer. (authors)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2004

Intelligent Control for the BEES Flyer

This paper describes the effort to provide a preliminary capability analysis and a neural network based adaptive flight control system for the JPL-led BEES aircraft project. The BEES flyer was envisioned to be a small, autonomous platform with sensing and control systems mimicking those of biological systems for the purpose of scientific exploration on the surface of Mars. The platform is physically tightly constrained by the necessity of efficient packing within rockets for the trip to Mars. Given the physical constraints, the system is not an ideal configuration for aerodynamics or stability and control. The objectives of this effort are to evaluate the aerodynamics characteristics of the existing design, to make recommendaaons as to potential improvements and to provide a control system that stabilizes the existing aircraft for nominal flight and damaged conditions. Towards this several questions are raised and analyses are presented to arrive at answers to some of the questions raised. CART3D, a high-fidelity inviscid analysis package for conceptual and preliminary aerodynamic design, was used to compute a parametric set of solutions over the expected flight domain. Stability and control derivatives were extracted from the database and integrated with the neural flight control system. The Integrated Vehicle Modeling Environment (IVME) was also used for estimating aircraft geometric, inertial, and aerodynamic characteristics. A generic neural flight control system is used to provide adaptive control without the requirement for extensive gain scheduling or explicit system identification. The neural flight control system uses reference models to specify desired handling qualities in the roll, pitch, and yaw axes, and incorporates both pre-trained and on-line learning neural networks in the inverse model portion of the controller. Results are presented for the BEES aircraft in the subsonic regime for terrestrial and Martian environments.

Krishnakumar, K.↗
DOE OSTI2025

Seismicity-constrained fault detection and characterization with a multitask machine learning model

Geological fault detection and characterization are crucial for understanding subsurface dynamics across scales. While methods for fault delineation based on either seismicity location analysis or seismic image reflector discontinuity are well-established, a systematic approach that integrates both data types remains absent. We develop a novel machine learning model that unifies seismic reflector images and seismicity location information to automatically identify geological faults and characterize their geometrical properties. The model encodes a seismic image and a seismicity location image separately, and fuses the encoded features with a spatial-channel attention fusion module to improve the learning of important features in both inputs. We design an automated strategy to generate high-quality synthetic training data and labels. To improve the realism of the seismicity location image, we include random seismicity noise and missing seismicity location associated with some of the faults. We validate the model’s efficacy and accuracy using synthetic data examples and two field data examples. Moreover, we show that fine-tuning the trained model with a small, domain-specific dataset enhances its fidelity for field data applications. The results demonstrate that integrating seismicity location and seismic images into a unified framework allows the end-to-end neural network to achieve higher fidelity and accuracy in delineating subsurface faults and their geometrical properties compared with image-only fault detection methods. Our approach offers an adaptive data-driven tool for geological fault characterization and seismic hazard mitigation, bridging the gap between seismicity location and image-based fault detection methods.

58 GEOSCIENCES↗
DOE OSTI2013

Clinical Implementation of Intrafraction Cone Beam Computed Tomography Imaging During Lung Tumor Stereotactic Ablative Radiation Therapy

Purpose: To develop and clinically evaluate a volumetric imaging technique for assessing intrafraction geometric and dosimetric accuracy of stereotactic ablative radiation therapy (SABR). Methods and Materials: Twenty patients received SABR for lung tumors using volumetric modulated arc therapy (VMAT). At the beginning of each fraction, pretreatment cone beam computed tomography (CBCT) was used to align the soft-tissue tumor position with that in the planning CT. Concurrent with dose delivery, we acquired fluoroscopic radiograph projections during VMAT using the Varian on-board imaging system. Those kilovolt projections acquired during millivolt beam-on were automatically extracted, and intrafraction CBCT images were reconstructed using the filtered backprojection technique. We determined the time-averaged target shift during VMAT by calculating the center of mass of the tumor target in the intrafraction CBCT relative to the planning CT. To estimate the dosimetric impact of the target shift during treatment, we recalculated the dose to the GTV after shifting the entire patient anatomy according to the time-averaged target shift determined earlier. Results: The mean target shift from intrafraction CBCT to planning CT was 1.6, 1.0, and 1.5 mm; the 95th percentile shift was 5.2, 3.1, 3.6 mm; and the maximum shift was 5.7, 3.6, and 4.9 mm along the anterior-posterior, left-right, and superior-inferior directions. Thus, the time-averaged intrafraction gross tumor volume (GTV) position was always within the planning target volume. We observed some degree of target blurring in the intrafraction CBCT, indicating imperfect breath-hold reproducibility or residual motion of the GTV during treatment. By our estimated dose recalculation, the GTV was consistently covered by the prescription dose (PD), that is, V100% above 0.97 for all patients, and minimum dose to GTV >100% PD for 18 patients and >95% PD for all patients. Conclusions: Intrafraction CBCT during VMAT can provide geometric and dosimetric verification of SABR valuable for quality assurance and potentially for treatment adaptation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

CAMELLIA

Camellia offers a framework for development of solvers for the discontinuous Petrov-Galerkin methodology of Demkowicz and Gopalakrishnan as well as other finite element methodologies. Spatially 1D-3D meshes are supported, and space-time elements are supported for all spatial topologies. Both h and p adaptivity are supported. A general interface for direct and iterative solvers is provided; a geometric multigrid preconditioner is also provided. The core computational kernels take advantage of MPI to execute in parallel.

Roberts, Nathan↗
DOE OSTI2024

Preserving Superconvergence of Spectral Elements for Curved Domains [Slides]

Finite Element Methods (FEM) and Spectral Element Methods (SEM) are crucial for solving partial differential equations (PDEs) on complex geometries. SEM offers superior accuracy due to potential superconvergence for simple domains. Challenges persist for domains with curved boundaries, restricting SEM’s advantages in real-world applications. A proposed solution is the introduction of a novel strategy to enhance accuracy and maintain superconvergence of SEM in curved domains. The strategy includes a mesh-generation procedure with geometrically refined elements near curved boundaries and a post-processing phase using the Adaptive Extended Stencil Finite Element Method (AES-FEM). The method, named AES-FEM post-processed Spectral Element Method (ApSEM), aligns the accuracy of non-tensor-product elements with superconvergent spectral elements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Price equation captures the role of drug interactions and collateral effects in the evolution of multidrug resistance

Bacterial adaptation to antibiotic combinations depends on the joint inhibitory effects of the two drugs (drug interaction [DI]) and how resistance to one drug impacts resistance to the other (collateral effects [CE]). Here we model these evolutionary dynamics on two-dimensional phenotype spaces that leverage scaling relations between the drug-response surfaces of drug-sensitive (ancestral) and drug-resistant (mutant) populations. We show that evolved resistance to the component drugs – and in turn, the adaptation of growth rate – is governed by a Price equation whose covariance terms encode geometric features of both the two-drug-response surface (DI) in ancestral cells and the correlations between resistance levels to those drugs (CE). Within this framework, mean evolutionary trajectories reduce to a type of weighted gradient dynamics, with the drug interaction dictating the shape of the underlying landscape and the collateral effects constraining the motion on those landscapes. We also demonstrate how constraints on available mutational pathways can be incorporated into the framework, adding a third key driver of evolution. Our results clarify the complex relationship between drug interactions and collateral effects in multidrug environments and illustrate how specific dosage combinations can shift the weighting of these two effects, leading to different and temporally explicit selective outcomes.

Gjini, Erida↗
NASA NTRS1973

Digital processing of Mariner 9 television data.

The digital image processing performed by the Image Processing Laboratory (IPL) at JPL in support of the Mariner 9 mission is summarized. The support is divided into the general categories of image decalibration (the removal of photometric and geometric distortions from returned imagery), computer cartographic projections in support of mapping activities, and adaptive experimenter support (flexible support to provide qualitative digital enhancements and quantitative data reduction of returned imagery). Among the tasks performed were the production of maximum discriminability versions of several hundred frames to support generation of a geodetic control net for Mars, and special enhancements supporting analysis of Phobos and Deimos images.

Green, W. B.↗
NASA NTRS1984

Using a commercial CAD system for simultaneous input to theoretical aerodynamic programs and wind-tunnel model construction

The Computer Aided Design (CAD) system's common geometry database was used to generate input for theoretical programs and numerically controlled (NC) tool paths for wind tunnel part fabrication. This eliminates the duplication of work in generating separate geometry databases for each type of analysis. Another advantage is that it reduces the uncertainty due to geometric differences when comparing theoretical aerodynamic data with wind tunnel data. The system was adapted to aerodynamic research by developing programs written in Design Analysis Language (DAL). These programs reduced the amount of time required to construct complex geometries and to generate input for theoretical programs. Certain shortcomings of the Design, Drafting, and Manufacturing (DDM) software limited the effectiveness of these programs and some of the Calma NC software. The complexity of aircraft configurations suggests that more types of surface and curve geometry should be added to the system. Some of these shortcomings may be eliminated as improved versions of DDM are made available.

Enomoto, F.↗
NASA NTRS1996

Solution-Adaptive Cartesian Cell Approach for Viscous and Inviscid Flows

A Cartesian cell-based approach for adaptively refined solutions of the Euler and Navier-Stokes equations in two dimensions is presented. Grids about geometrically complicated bodies are generated automatically, by the recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, polygonal cut cells are created using modified polygon-clipping algorithms. The grid is stored in a binary tree data structure that provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite volume formulation. The convective terms are upwinded: A linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The results of a study comparing the accuracy and positivity of two classes of cell-centered, viscous gradient reconstruction procedures is briefly summarized. Adaptively refined solutions of the Navier-Stokes equations are shown using the more robust of these gradient reconstruction procedures, where the results computed by the Cartesian approach are compared to theory, experiment, and other accepted computational results for a series of low and moderate Reynolds number flows.

Coirier, William J.↗
NASA NTRS1998

The AVIRIS Low Altitude Option-An Approach to Increase Geometric Resolution and Improve Operational Flexibility Simultaneously

From 1987 through 1997 the Airborne Visible-InfraRed Imaging Spectrometer has matured into a remote sensing instrument capable of producing prodigious amounts of high quality data. Using the NASA/Ames ER-2 high altitude aircraft platform, flight operations have become very reliable as well. Being exclusively dependent on the ER-2, however, has limitations: the ER-2 has a narrow cruise envelope which fixes the AVIRIS ground pixel at 20 meters; it requires a significant support infrastructure; and it has a very limited number of bases it can operate from. In the coming years, the ER-2 will also become less available for AVIRIS flights as NASA Earth Observing System satellite underflights increase. Adapting AVIRIS to lower altitude, less specialized aircraft will create a much broader envelope for data acquisition, i.e., higher ground geometric resolution while maintaining nearly the ideal spatial sampling. This approach will also greatly enhance flexibility while decreasing the overall cost of flight operations and field support. Successful adaptation is expected to culminate with a one-month period of demonstration flights.

Sarture, Charles M.↗
DOE OSTI2019

Dependence of Lamellar Eutectic Growth with Convection on Boundary Conditions and Geometric Confinement: A Phase-Field Lattice-Boltzmann Study

Simulating eutectic growth with convection is challenging because of the enormous computing demand resulting from the required domain size compared with the eutectic scale. In this study, a phase-field lattice-Boltzmann approach was employed to reproduce the eutectic growth with convection (both natural and forced). To enable and accelerate large-scale simulations, a parallel-adaptive mesh refinement algorithm was developed to combine with this numerical approach. The effects of the boundary conditions and the geometric confinement on the eutectic growth with convection were systematically investigated in the Al-Cu eutectic alloy. Results showed that both boundary conditions and the geometric confinement altered the eutectic evolution by affecting the distributions of the solute and velocity ahead of the eutectic solid/liquid interface. As the first attempt to explore the two factors on the eutectic growth with convection, our thorough investigation lays a foundation for modeling the coupled convection–solute eutectic evolution.

36 MATERIALS SCIENCE↗