Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational geometry”

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 397 records · Page 22

An Immersed Boundary Method for Solving the Compressible Navier-Stokes Equations with Fluid Structure Interaction

An immersed boundary method for the compressible Navier-Stokes equation and the additional infrastructure that is needed to solve moving boundary problems and fully coupled fluid-structure interaction is described. All the methods described in this paper were implemented in NASA's LAVA solver framework. The underlying immersed boundary method is based on the locally stabilized immersed boundary method that was previously introduced by the authors. In the present paper this method is extended to account for all aspects that are involved for fluid structure interaction simulations, such as fast geometry queries and stencil computations, the treatment of freshly cleared cells, and the coupling of the computational fluid dynamics solver with a linear structural finite element method. The current approach is validated for moving boundary problems with prescribed body motion and fully coupled fluid structure interaction problems in 2D and 3D. As part of the validation procedure, results from the second AIAA aeroelastic prediction workshop are also presented. The current paper is regarded as a proof of concept study, while more advanced methods for fluid structure interaction are currently being investigated, such as geometric and material nonlinearities, and advanced coupling approaches.

Equations↗

Flyover-Noise Reduction of Commercial Aircraft via Reduced Slat Deflections and Modified Flight Procedures

Current regulation and operation of commercial aircraft result in flyover noise that generates public protest from residents well outside the recognized noise footprint of modern airports. This noise issue was recently exacerbated by the implementation of NextGen air-traffic control which uses GPS-navigation for predictable flight paths and improved efficiency. A side-effect of the narrow flight corridors is a perceived increase in noise level and annoyance experienced by the residents directly under the flight paths. As a possible solution to the community noise issue, the feasibility of a lower slat deflections and modified flight procedures were explored. To document the maximum potential noise reduction, field measurements were conducted on the ground at locations before and after slat deployment. For safety considerations, the stall margin of a lower slat deflection was computed for a realistic geometry based on a modern commercial aircraft. The logistics of setting lower slat angles on current aircraft were also investigated through discussions with the FAA (Federal Aviation Administration) and aircraft manufacturers. The results suggest that community flyover noise could be reduced via lower slat deflections during early approach, but that modifying current aircraft would be difficult due to existing control systems and hurdles in certification. Alternatively, modifications to existing flight procedures could reduce the perceived noise of the NextGen flight patterns. For next generation commercial aircraft, reduced slat deflections could be considered if mandated by FAA noise requirements during early approach.

Storms, Bruce L.↗

Simulation and Analysis of Opportunistic MSPA for Multiple Cubesat Deployments

We describe a software approach for simultaneous demodulation and decoding of multiple frequency- multiplexed spacecraft across multiple ground stations. The approach involves the use of a single ground antenna with wide enough aperture to receive multiple angularly adjacent spacecraft simultaneously. This technique uses a wide-band RF signal digitizer coupled with easily and cheaply duplicated software-receiver modules to independently process the frequency-channelized downlink signals from the spacecraft. This approach relaxes the need for realizable modems at each antenna, for example in the Deep-SpaceNetwork (DSN), and can also function as a delayed data retrieval method for cubesat missions with routine science data return. Thus, we hope that this solution can enable more efficient utilization of DSN assets. We concentrate on a simulation similar to the proposed Exploration Mission 1 (EM-1) mission from a geometric perspective. To that end, the simulated scenario involves ten secondary cubesats deployed and tracked for a span of four days (the modeled motion over the first four days of the mission does not include any Trajectory Correction Maneuvers (TCMs) that may be necessary as the cubesats approach the moon). Transmissions from the cubesats may be received through a combination of DSN ground sites based on visibility. The cubesats are assumed to utilize typical cubesat transmit powers and a waveform similar to that of the Iris radio. One cubesat is chosen as the “target” and is considered tracked by all ground-stations throughout the simulation (i.e., it is consistently at the center of the main lobe of each ground station antenna when the ground station is in view of the cubesat in terms of elevation angle). The simulation effort involves synthesizing a wideband signal that includes the ten cubesats across a large bandwidth due to each cubesat having its own center-frequency in X-band, near 8.4 GHz. Each signal is characterized by its own Doppler-frequency shift and free-space path-loss computed through the underlying geometry of the simulation as well as the cubesat transmit power. Finally, each cubesat signal experiences a unique antenna gain at each ground station due to the underlying antenna pattern and spacecraft-ground station geometry. We establish two interesting findings: First, the link-budget for the EM-1-like scenario is almost completely limited by the angle between the center of the antenna main beam and the cubesat, which means that the signal-to-noise ratio is wholly adequate for demodulation and decoding at the simulated bitrates unless the cubesat exits the main beam. Secondly, and fortunately due to the use of a software radio architecture, we show that it is possible to successfully receive signals from most cubesats for the entire 4-day simulation. Due to the signal-to-noise ratio being sufficient for demodulation even after the side-lobe’s 17dB reduction in SNR, the software radio can still demodulate such cubesats as long as the radio is capable of re-establishing carrier lock as the cubesats leave the main beam and enter the side-lobe. Considering that the target application utilizes offline processing, a loss of lock can be detected and lock can be re-established by iterating over the data. Extensive simulations demonstrate these results.

Abraham, Douglas S.↗

Rapid Hypersonic Simulations using US3D and Pointwise

For hypersonic simulations, unstructured flow solvers typically have problems predicting surface heat fluxes when strong shocks are present. To address these issues, this paper outlines a workflow that applies best practices developed for structured grids to unstructured meshes. In addition, unstructured grid generation can significantly reduce the time required to create quality grids for complex geometries. Several examples are computed using DPLR, a structured grid flow solver, and US3D flow, an unstructured mesh solver. Results from the two codes are compared, and they show excellent agreement. Overall, the unstructured grid workflow offers a viable and attractive alternative for hypersonic simulations.

C. Tang↗

Rapid Hypersonic Simulations using US3D and Pointwise

For hypersonic simulations, unstructured flow solvers typically have problems predicting surface heat fluxes when strong shocks are present. To address these issues, this paper outlines a workflow that applies best practices developed for structured grids to unstructured meshes. In addition, unstructured grid generation can significantly reduce the time required to create quality grids for complex geometries. Several examples are computed using DPLR, a structured grid flow solver, and US3D flow, an unstructured mesh solver. Results from the two codes are compared, and they show excellent agreement. Overall, the unstructured grid workflow offers a viable and attractive alternative for hypersonic simulations.

Chun Tang↗

Hypersonic Simulations with US3D using Unstructured Grids from Fidelity Pointwise

For hypersonic simulations, unstructured flow solvers typically have problems predicting surface heat fluxes when strong shocks are present. This article outlines a workflow that applies best practices for structured and unstructured grids. In addition, unstructured grid generation can significantly reduce the time required to create quality grids for complex geometries. Several examples are computed using DPLR, a structured grid flow solver, and US3D flow, an unstructured flow solver. Results from the two codes are compared and they show excellent agreement. The unstructured grid workflow offers a viable and attractive alternative for hypersonic simulations.

C. Tang↗

Plug20 Test Report: Particle Image Velocimetry (PIV) Measurements

This report provides turbulent flow statistics made using particle image velocimetry (PIV) for several internally mixed exhaust systems with external plugs (Plug20 designs). The tests were conducted at the NASA Glenn Research Center’s Aero-Acoustic Propulsion Laboratory (AAPL) in May and June 2022. The corresponding acoustic data for the configurations and flows presented in this report can be found in an earlier Technical Memorandum, NASA/TM-20210010291, “Plug20 Test Report.” The Plug20 nozzle family was designed to represent commercial supersonic vehicles in the 2025 to 2035 timeframe with a bypass ratio around 4, having an internal mixer, with and without an external plug to improve nacelle boattail drag and shock. The flow conditions measured represent engine cycles for commercial supersonic aircraft at takeoff, having pressure ratios of 2.0 and heated fan and core streams, measured both with and without a flight stream surrounding the plume. A few unheated flow conditions were also tested to tie results to other baseline test results. The PIV measurements in both streamwise and cross-stream planes are described, along with brief analysis to visually catalog the contents and confirm their validity. In this analysis, the streamwise PIV is found to match benchmark datasets within its experimental error. However, the cross-stream turbulent quantities measured by two systems do not match each other or the streamwise measurements in amplitude—caution is required when using the cross-stream data. Data acquired and included in this report are (1) velocity statistics of the plumes in both a streamwise and several cross-stream planes, and (2) the hot nozzle geometry. Tabular forms of the data and computer-aided design (CAD) files of the hot geometry are presented as supplemental files described in Appendix A (TM-20230006779SUPPL1.zip) and Appendix B (TM-20230006779SUPPL2.zip).

Turbulence↗

Empirical Noise Modeling of Internally Mixed Exhaust Systems

There appear to be no non-proprietary methods to predict the noise of internally mixed exhaust systems, and no guidance for how to adapt known jet noise models for these configurations. This paper surveys literature and historical databases acquired at NASA Glenn's Aero-Acoustic Propulsion Lab to give such guidance. The core premise is that an exhaust system with well-designed mixer produces noise that is to first approximation the same as a fully-mixed jet flow. Refinement of what is meant by a “fully mixed jet” can lead to more accurate prediction of the main jet noise. Additional noise is often generated within the nozzle, typically at high frequencies, whose source mechanism(s) are not obvious. However, a noise prediction method can be established that captures some aspects of the excess noise and provide an estimate of the total jet noise. Explorations of source mechanisms associated with the internal mixer have led to a new noise model which includes the effects of having an external plug nozzle, a feature desirable for near-term supersonic aircraft. Statistical analysis of historical data is provided to estimate the uncertainty in using this method given the variations found that cannot be directly computed without detailed mixer geometry.

Noise Prediction↗

Empirical Noise Modeling of Internally Mixed Exhaust Systems

There appears to be no non-proprietary methods to predict the noise of internally mixed exhaust systems, and no guidance for how to adapt known jet noise models for these configurations. This paper surveys literature and historical databases acquired at NASA Glenn's Aero-Acoustic Propulsion Lab to give such guidance. The core premise is that an exhaust system with well-designed mixer produces noise that is to first approximation the same as a fully-mixed jet flow. Refinement of what is meant by a “fully mixed jet” can lead to more accurate prediction of the main jet noise. Additional noise is often generated within the nozzle, typically at high frequencies, whose source mechanism(s) are not obvious. However, a noise prediction method can be established that captures some aspects of the excess noise and provide an estimate of the total jet noise. Explorations of source mechanisms associated with the internal mixer have led to a new noise model which includes the effects of having an external plug nozzle, a feature desirable for near-term supersonic aircraft. Statistical analysis of historical data is provided to estimate the uncertainty in using this method given the variations found that cannot be directly computed without detailed mixer geometry.

Noise Prediction↗

Aerodynamic Design Optimization for Natural Laminar Flow Airfoils

Natural laminar flow technology is a passive laminar flow control (LFC) strategy that seeks to delay the onset of boundary-layer transition (BLT) through shape optimization to reduce the drag of the aerodynamic vehicle. Adjoint-based design optimization for LFC is proposed in an integrated multidisciplinary framework, which includes the computational fluid dynamics (CFD), geometry and grid deformation, and linear stability analysis (LSA) for transition prediction. In particular, the BLT location is predicted using the dual N-factor method that is based on a linear stability theory (LST) eigenvalue problem. The dual N-factor criterion accounts for the amplification of planar Tollmien-Schlichting (TS) and stationary crossflow (CF) boundary-layer instabilities to predict the transition location in three-dimensional boundary-layer flows. The adjoint-based shape optimization procedure is based on an iteratively coupled CFD and LSA methodology to converge the transition location and flow solutions, as well as to calculate the sensitivities of the aerodynamic metrics of interest with respect to the flow and shape design parameters. The RAE 2822 airfoil at 0 and 30 degrees yaw angles, an angle of attack of 0.72 degrees, and subsonic conditions (M∞ = 0.19, Rec = 5.6 × 106 ) are used as baseline configurations for design optimization. The angle of attack and the vertical displacement of free-form-deformation control points are used as design variables to reduce the drag coefficient while reaching a specified lift coefficient. The optimized unswept airfoil designs achieve a 30% drag reduction accompanied by a downstream shift of the transition locations over both suction and pressure sides of the airfoil. The initial design iterations for the swept case also show a favorable trend in the drag reduction with transition delay over both sides.

Transition↗

Comparison of Geant4-DNA and RITRACKS/RITCARD: Microdosimetry, Nanodosimetry and DNA Damage Predictions for Low to High LET Ions

Low linear energy transfer (LET) ionizing radiations induce homogeneously distributed DNA damage at the cellular scale, while high LET ionizing radiations are well known to create complex biological damage including clustered double strand breaks (DSB). Clustered DSB are critical events that are difficult to repair, leading to increased biological outcomes compared to low LET radiation exposure. Recent studies have shown that DNA compaction, which can be actively modified by the cell to allow DNA repair, influences the yield of DSB [1]. The goal of this work is to determine the impact of DNA compaction on chromosome aberration formation, a biomarker of radiation-induced cancer risk. The RITRACKS track structure code [2], together with the RITCARD code [3], allow one to simulate the transport of high LET ions, calculate DNA damage and repair, and chromosome aberration formation, but does not include detailed geometrical models of DNA. The Geant4-DNA toolkit provides such detailed models with different levels of DNA compaction [1,4], but does not model chromosome aberration formation. In this work, these highly detailed DNA geometries are used to compute DNA breaks with Geant4-DNA and combine the results with the RITCARD tool to compute chromosome aberrations. As a first step, we compared the transport models used in RITRACKS and Geant4-DNA and newly released heavy ion model transport (G4DNARuddIonisationExtendedModel Geant4 v11.2) to ensure consistency between the physical stage of the two radiation transport codes. We compared energy deposition at the micrometric and nanometric scales for different ions (H 250 MeV, H 150 Me, He 250 MeV/n, C 290 MeV/n, O 350 MeV/n, O 55 MeV/n, Si 170MeV/n, Ti 300 MeV/n, Ti 300 MeV/n, Fe 600 MeV/n, Fe 450 MeV/n and Fe 300 MeV/n) with LET ranging from ~0.4 keV/μm up to 235 keV/μm. Excellent agreement is found for both microdosimetric and nanodosimetric spectra for all ion types and energies. The calculation of DNA breaks and their complexity with the Geant4-DNA is an ongoing work, and comparison with RITCARD will be presented based on result availability at the time of the conference. [1] Tang et al. (2019). Med. Phys., 46(3), 1501-1511. [2] Plante et al. (2008). NJP, 10(12), 125020. [3] Plante et al. (2019). Rad. Res., 192(3), 282-298. [4] Incerti et al. (2018). Med. Phys., 45(8), e722-e739.

Floriane Poignant↗

Revisiting Topographic Horizons in the Era of Big Data and Parallel Computing

Widely used to calculate illumination geometry forestimates of solar and emitted longwave radiation, and forcorrecting remotely sensed data for topographic effects, digitalelevation models (DEMs) are now extensive globally at 10–30-mspatial resolution and locally at spatial resolutions down to afew centimeters. Globally, regionally, or locally, elevation datasetshave many grid points. Many software packages calculate gradi-ents over every grid cell or point, but in the mountains, shadingby nearby terrain must also be assessed. Terrain may obscure aslope that would otherwise face the Sun. Four decadesago, a fastmethod to calculate topographic horizons at every point in anelevation grid required computations related only linearly to thesize of the grid, but grids now have so many points that parallelcomputing still provides an advantage. Exploiting parallelismover terrain grids can use alternative strategies: among columnsof a rotated grid, or simultaneously at multiple rotation angles,or on different tiles of a grid. On a multi-processor machine, theimprovement in computing time approaches 2/3 the number ofprocessors deployed,

Jeff Dozier↗

Parametric analysis of swept-wing geometry with sheared wing tips

A computational parameter study is presented of potential reductions in induced drag and increases in lateral-directional stability due to sheared wing tips attached to an untwisted wing of moderate sweep and aspect ratio. Sheared tips are swept and tapered wing-tip devices mounted in the plane of the wing. The induced-drag results are obtained using an inviscid, incompressible surface-panel method that models the nonlinear effects due to the deflected and rolled-up wake behind the lifting surface. The induced-drag results with planar sheared tips are compared to straight-tapered tip extensions and nonplanar winglet geometries. The lateral-directional static-stability characteristics of the wing with sheared tips are estimated using a quasi-vortex-lattice method. For certain combinations of sheared-tip sweep and taper, both the induced efficiency of the wing and the relevant static-stability derivatives are predicted to increase compared to the wing with a straight-tapered tip modification.

Fremaux, C. M.↗

The effect of ionization on the infrared absorption spectra of PAHs: A preliminary report

The emission lines observed in many interstellar IR sources at 3.28, 6.2, 7.7, 8.7, and 11.3 microns are theorized to originate from polycyclic aromatic hydrocarbons (PAHs). These assignments are based on analyses of lab IR spectra of neutral PAHs. However, it is likely that in the interstellar medium that PAHs are ionized, i.e., are positively charged. Besides, as pointed out by Allamandola et al., although the IR emission band spectrum resembles what one might expect from a mixture of PAHs, it does not match in details such as frequency, band profile, or relative intensities predicted from the absorption spectra of any known PAH molecule. One source of more information to test the PAH theory is ab initio molecular orbital theory. It can be used to compute, from first principles, the geometries, vibrational frequencies, and vibrational intensities for model PAH compounds which are difficult to study in the lab. The Gaussian 86 computer program was used to determine the effect of ionization on the infrared absorption spectra of several small PAHs: naphthalene and anthracene. A preliminary report is presented of the results of these calculations.

Defrees, Doug J.↗

A Formally-Verified Decision Procedure for Univariate Polynomial Computation Based on Sturm's Theorem

Sturm's Theorem is a well-known result in real algebraic geometry that provides a function that computes the number of roots of a univariate polynomial in a semiopen interval. This paper presents a formalization of this theorem in the PVS theorem prover, as well as a decision procedure that checks whether a polynomial is always positive, nonnegative, nonzero, negative, or nonpositive on any input interval. The soundness and completeness of the decision procedure is proven in PVS. The procedure and its correctness properties enable the implementation of a PVS strategy for automatically proving existential and universal univariate polynomial inequalities. Since the decision procedure is formally verified in PVS, the soundness of the strategy depends solely on the internal logic of PVS rather than on an external oracle. The procedure itself uses a combination of Sturm's Theorem, an interval bisection procedure, and the fact that a polynomial with exactly one root in a bounded interval is always nonnegative on that interval if and only if it is nonnegative at both endpoints.

Narkawicz, Anthony J.↗

Performance Analysis of an Optimization Algorithm for Metamaterial Design on the Integrated High-Performance Computing and Quantum Systems

Optimizing metamaterials with complex geometries is a big challenge. Although an active learning algorithm, combining machine learning (ML), quantum computing, and optical simulation, has emerged as an efficient optimization tool, it still faces difficulties in optimizing complex structures that have potentially high performance. In this work, we comprehensively analyze the performance of an optimization algorithm for metamaterial design on the integrated HPC and quantum systems. We demonstrate significant time advantages through message-passing interface (MPI) parallelization on the high-performance computing (HPC) system showing approximately 54% faster ML tasks and 67 times faster optical simulation against serial workloads. Furthermore, we analyze the performance of a quantum algorithm designed for optimization, which runs with various quantum simulators on a local computer or HPC-quantum system. Results showcase ~24 times speedup when executing the optimization algorithm on the HPC-quantum hybrid system. This study paves a way to optimize complex metamaterials using the integrated HPC-quantum system.

Kim, Seongmin↗

Towards Accurate Prediction of Turbulent, Three-Dimensional, Recirculating Flows with the NCC

The National Combustion Code (NCC) was used to calculate the steady state, nonreacting flow field of a prototype Lean Direct Injection (LDI) swirler. This configuration used nine groups of eight holes drilled at a thirty-five degree angle to induce swirl. These nine groups created swirl in the same direction, or a corotating pattern. The static pressure drop across the holes was fixed at approximately four percent. Computations were performed on one quarter of the geometry, because the geometry is considered rotationally periodic every ninety degrees. The final computational grid used was approximately 2.26 million tetrahedral cells, and a cubic nonlinear k - epsilon model was used to model turbulence. The NCC results were then compared to time averaged Laser Doppler Velocimetry (LDV) data. The LDV measurements were performed on the full geometry, but four ninths of the geometry was measured. One-, two-, and three-dimensional representations of both flow fields are presented. The NCC computations compare both qualitatively and quantitatively well to the LDV data, but differences exist downstream. The comparison is encouraging, and shows that NCC can be used for future injector design studies. To improve the flow prediction accuracy of turbulent, three-dimensional, recirculating flow fields with the NCC, recommendations are given.

Iannetti, A.↗

The implementation of the graphics of program EAGLE: A numerical grid generation code on NASA Langley SNS computer system

Program EAGLE (Eglin Arbitrary Geometry Implicit Euler) Numerical Grid Generation System is a composite (multi-block) algebraic or elliptic grid generation system designed to discretize the domain in and/or around any arbitrarily shaped three dimensional regions. This system combines a boundary conforming surface generation scheme and includes plotting routines designed to take full advantage of the DISSPLA Graphics Package (Version 9.0). Program EAGLE is written to compile and execute efficiently on any Cray machine with or without solid state disk (SSD) devices. Also, the code uses namelist inputs which are supported by all Cray machines using the FORTRAN compiler CFT77. The namelist inputs makes it easier for the user to understand the inputs and operation of Program EAGLE. EAGLE's numerical grid generator is constructed in the following form: main program, EGG (executive routine); subroutine SURFAC (surface generation routine); subroutine GRID (grid generation routine); and subroutine GRDPLOT (grid plotting routines). The EAGLE code was modified to use on the NASA-LaRC SNS computer (Cray 2S) system. During the modification a conversion program was developed for the output data of EAGLE's subroutine GRID to permit the data to be graphically displayed by IRIS workstations, using Plot3D. The code of program EAGLE was modified to make operational subroutine GRDPLOT (using DI-3000 Graphics Software Packages) on the NASA-LaRC SNS Computer System. How to implement graphically, the output data of subroutine GRID was determined on any NASA-LaRC graphics terminal that has access to the SNS Computer System DI-300 Graphics Software Packages. A Quick Reference User Guide was developed for the use of program EAGLE on the NASA-LaRC SNS Computer System. One or more application program(s) was illustrated using program EAGLE on the NASA LaRC SNS Computer System, with emphasis on graphics illustrations.

Houston, Johnny L.↗