SEARCH · Search NASA
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.
Los Alamos National Laboratory (LANL) is currently working on 3 new production applications, VPC, xRage, and Pagosa. VPIC was designed to be a 3D relativist, electromagnetic Particle-In-Cell code for plasma simulation. xRage, a 3D AMR mesh amd multi physics hydro code. Pagosa, is a 3D structured mesh and multi physics hydro code.
The height of scattering phase retrieved from InSAR data is considered being correlated with the tree height and the spatial structure of the forest stand. Though some researchers have used simple backscattering models to estimate tree height from the height of scattering center, the effect of forest spatial structure on InSAR data is not well understood yet. A three-dimensional coherent radar backscattering model for forest canopies based on realistic three-dimensional scene was used to investigate the effect in this paper. The realistic spatial structure of forest canopies was established either by field measurements (stem map) or through use of forest growth model. Field measurements or a forest growth model parameterized using local environmental parameters provides information of forest species composition and tree sizes in certain growth phases. A fractal tree model (L-system) was used to simulate individual 3- D tree structure of different ages or heights. Trees were positioned in a stand in certain patterns resulting in a 3-D medium of discrete scatterers. The radar coherent backscatter model took the 3-D forest scene as input and simulates the coherent radar backscattering signature. Interferometric SAR images of 3D scenes were simulated and heights of scattering phase centers were estimated from the simulated InSAR data. The effects of tree height, crown cover, crown depth, and the spatial distribution patterns of trees on the scattering phase center were analyzed. The results will be presented in the paper.
The objective of this work is the development of computational models and analysis of the coupled fluid-structure response of a slat gap filler (SGF) noise treatment applied to the leading-edge-slat component of a high-lift system typical of modern transport aircraft. The representative airframe chosen for this work is NASA’s High-lift Common Research Model (CRM-HL) in a baseline high-lift configuration. Superelastic shape memory alloys (SMAs) have been identified as enabling materials for these structural treatments. Since the technology elements rely upon having a highly reconfigurable structure, designs must be assessed for static aeroelastic deflection as well as dynamic aeroelastic stability using coupled computational fluid dynamics (CFD) and nonlinear computational structural dynamics (NL CSD) tools. The technical approach consists of solving for the flow field around the entire vehicle using a global CFD model, followed by extraction of relevant local subdomain data for CFD and NL-CSD co-simulations. The SGF design is assessed using both 2D and 3D co-simulations to predict quasi-static aeroelastic deformations and to assess dynamic aeroelastic stability.
This final report documents work performed by ATA Engineering, Inc., (ATA) to develop computational models and analyze the coupled fluid-structure response of two types of noise treatments applied to the leading-edge-slat component of a high-lift system typical of modern transport aircraft. The first treatment is a slat-gap filler (SGF), which closes the gap between the suction surfaces of a deployed slat and an aircraft main wing, and the second treatment is a slat-cove filler (SCF), which replaces the recirculating flow on the slat cove with a surface that promotes flow attachment. The representative airframe chosen for this work was NASA’s High-Lift Common Research Model (CRM-HL) in a baseline high-lift configuration. Superelastic shape memory alloys (SMAs) have been identified as enabling materials for these structural treatments. Since the technology elements rely upon having a highly reconfigurable structure, designs were assessed for their static aeroelastic deflection as well as their dynamic aeroelastic stability using coupled computational fluid dynamics (CFD) and nonlinear computational structural dynamics (NL-CSD) tools. Specifically, fluid-structure interaction (FSI) problems were solved computationally using the CFD solver Loci/CHEM and the NL-CSD solver Abaqus. As a part of the overall project, a similar capability was implemented using the CFD solver FUN3D coupled to Abaqus, although that work is documented in a separate report and that FSI framework was not used to analyze any of NASA’s SGF and SCF configurations. The technical approach consisted of solving for the flow field around the entire vehicle using a global CFD model, followed by extraction of relevant local subdomain data for CFD and NL-CSD co-simulations. The SGF design was analyzed using both 2D and 3D co-simulations to predict quasi-static aeroelastic deformations and to assess dynamic aeroelastic stability, whereas the SCF was analyzed in 2D only. SGF static aeroelastic response predictions focused on characterizing the deformed shape, with maximum displacements predicted to be on the order of magnitude of the technology element panel thickness. SGF dynamic aeroelastic response predictions used Partial Floquet analysis of the temporal evolution of selected nodal displacements to quantify the sign and magnitude of aeroelastic damping. Results suggest that the CRM-HL operating conditions would result in a dynamically stable response. The simulated dynamic pressure was also increased up to a factor of about four, and resulting responses suggest that predicted dynamic stability would be achieved with some margin.
Distortion and buildup of residual stresses are common design problems in structures manufactured using additive manufacturing techniques, popularly known as 3D printing. These problems lead to poor product quality, which often is improved by trial and error, where the sensitivity of different manufacturing configurations to product quality are experimentally determined. With the introduction of Big Area Additive Manufacturing (BAAM), where significant amount of material is deposited per unit time, trial and error philosophy reduces the economical efficiency of the process. In this work, a numerical simulation methodology based on Finite Element Method (FEM), multi-scale damage mechanics and fracture mechanics is introduced to simulate the BAAM process for determining product quality in terms of distortions, material damage and interface fracture due to manufacturing. In addition, this tool can be used to obtain the sensitivities of the manufacturing configurations with respect to product quality in an economically efficient manner compared to the trial and error philosophy. Computational tools providing a fully coupled thermo-mechanical solution of polymer additive manufacturing with chopped fiber reinforced plastics are presented.
This project expands upon some of the FY23 IRAD funded "Full Lunar Surface Simulation Platform" (FLS-SP) derived technologies. The results have demonstrated a low-cost, compact, immersive technology that supports training, interaction, mission planning, and other man-in-loop activities in simulated environments called “3D Stereo Immersion Environment Utilizing Mobile Devices” (3D SIEUMD). These activities will become crucial for safety and mission, planning, and assurance in the Artemis program, as well as future commercial crew programs. The FLS-SP enabled the visualization of simulated lunar surface operations for rovers, landers, habitation, and EVA. The 3D SIEUMD technology allows for NASA S&MA experts, tasked with defining safety protocols and procedures, to interact in an immersive environment without expensive immersion gear. Reduction of the overall safety risk by immersing safety concepts into the design and engineering phase for very small cost is accomplished even before hardware is deployed to the lunar surface.
In this study, a new simulation model for tokamak boundary plasma, SOLT3D, is implemented in the BOUT++ framework. The simulation model includes a set of dynamic equations describing collisional boundary plasma and neutral gas in the tokamak scrape-off layer and divertor region. The model is verified against standard linear plasma instabilities and available nonlinear results. For L-mode like conditions, SOLT3D reproduces characteristics of boundary plasma turbulent fluctuations that are similar to published experimental data, in terms of the amplitude and spatial dependence of the fluctuations. It also reproduces realistic plasma fluxes on material surfaces and expected Bohm-like effective radial transport. Plasma fluctuations inferred from the simulations imply inevitably a significant level of intrinsic 'noise' for axisymmetric tokamak plasma transport modeling, introducing errors on the order of unity. In particular, the toroidally averaged atomic rates below 5–10 eV are strongly modified by turbulent plasma fluctuations, which should significantly affect the standard axisymmetric modeling of the tokamak edge plasma and divertor.
Purpose: Cone-beam CT (CBCT) is a widely used imaging modality for image-guided radiotherapy. Most vendors provide CBCT systems that are mounted on a linac gantry. Thus, CBCT can be used to estimate the actual 3-dimensional (3D) position of moving respiratory targets in the thoracic/abdominal region using 2D projection images. The authors have developed a method for estimating the 3D trajectory of respiratory-induced target motion from CBCT projection images using interdimensional correlation modeling. Methods: Because the superior–inferior (SI) motion of a target can be easily analyzed on projection images of a gantry-mounted CBCT system, the authors investigated the interdimensional correlation of the SI motion with left–right and anterior–posterior (AP) movements while the gantry is rotating. A simple linear model and a state-augmented model were implemented and applied to the interdimensional correlation analysis, and their performance was compared. The parameters of the interdimensional correlation models were determined by least-square estimation of the 2D error between the actual and estimated projected target position. The method was validated using 160 3D tumor trajectories from 46 thoracic/abdominal cancer patients obtained during CyberKnife treatment. The authors’ simulations assumed two application scenarios: (1) retrospective estimation for the purpose of moving tumor setup used just after volumetric matching with CBCT; and (2) on-the-fly estimation for the purpose of real-time target position estimation during gating or tracking delivery, either for full-rotation volumetric-modulated arc therapy (VMAT) in 60 s or a stationary six-field intensity-modulated radiation therapy (IMRT) with a beam delivery time of 20 s. Results: For the retrospective CBCT simulations, the mean 3D root-mean-square error (RMSE) for all 4893 trajectory segments was 0.41 mm (simple linear model) and 0.35 mm (state-augmented model). In the on-the-fly simulations, prior projections over more than 60° appear to be necessary for reliable estimations. The mean 3D RMSE during beam delivery after the simple linear model had established with a prior 90° projection data was 0.42 mm for VMAT and 0.45 mm for IMRT. Conclusions: The proposed method does not require any internal/external correlation or statistical modeling to estimate the target trajectory and can be used for both retrospective image-guided radiotherapy with CBCT projection images and real-time target position monitoring for respiratory gating or tracking.
We have introduced in a simple and efficient manner quantum mechanical corrections in our 3D 'atomistic' MOSFET simulator using the density gradient formalism. We have studied in comparison with classical simulations the effect of the quantum mechanical corrections on the simulation of random dopant induced threshold voltage fluctuations, the effect of the single charge trapping on interface states and the effect of the oxide thickness fluctuations in decanano MOSFETs with ultrathin gate oxides. The introduction of quantum corrections enhances the threshold voltage fluctuations but does not affect significantly the amplitude of the random telegraph noise associated with single carrier trapping. The importance of the quantum corrections for proper simulation of oxide thickness fluctuation effects has also been demonstrated.
Although remote sensing (RS) of solar-induced chlorophyll fluorescence (SIF) is increasingly used as a valuable source of information about vegetation photosynthetic activity, the RS SIF observations are significantly influenced by canopy-specific structural features (i.e., canopy architecture including leaf area index and presence of woody components), atmospheric conditions during their acquisition (e.g., proportion of direct and diffuse irradiance) and observational geometric configurations (e.g., sun and viewing directions). Radiative transfer (RT) models have the potential to provide a better understanding of the canopy structural effects on the SIF emission and RS signals. Here, we used the DART model to assess the daily influence, from morning to evening, of forest 3D architecture on SIF nadir radiance, emission, escape factor and nadir yield of eight 100 m × 100 m forest study plots established in a temperate deciduous forest of the Smithsonian Environmental Research Center (Edgewater, MD, USA). The 3D architecture of each plot was derived from airborne LiDAR. DART simulations of these 3D forest plots and their 1D (i.e., vertical profile of sun-adapted and shade-adapted leaves) and 0D (i.e., homogeneous layer of sun-adapted leaves above an homogeneous layer of shade-adapted leaves) abstractions were compared to assess the relative errors (ε1D−3D and ε0D−3D) associated with horizontal and vertical structural heterogeneity, respectively. Forest 3D structure, especially horizontal heterogeneity, had a great influence on forest nadir SIF radiance, resulting in ε1D−3D up to 55% at 8:00 and 18:00 (i.e., for oblique sun directions). The key indicators of this impact, in the descending order of importance, were the SIF escape factor (ε1D−3D up to 40%), the attenuation of incident photosynthetically active radiation (ε1D−3D less than 5%), and the SIF emission yield (ε1D−3D less than 2%). The influence of forest architecture on the nadir SIF escape factor and SIF yield (ε1D−3D up to 40%) varied over time, with differences in forest stand structure, and per spectral domain, being always larger between 640 and 700 nm than between 700 and 850 nm. In addition, woody elements demonstrated a large influence on forest SIF radiance due to their “shading” effect (ε up to 17%) and their “blocking” effect (ε ≈ 10%), both of them higher for far-red than for red SIF. These results underline the importance of 3D forest canopy architecture, especially 2D heterogeneity, and inclusion of woody elements in RT modeling used for interpretation of the RS SIF signal, and subsequently for the estimation of gross primary production and detection of vegetation stress.
Massive stars can end their lives with a successful supernova explosion (leaving behind a neutron star or, more rarely, a black hole), or a failed explosion that leaves behind a black hole. The density structure of the pre-collapse progenitor star already encodes much of the information regarding the outcome and properties of the explosion. However, the complexity of the collapse and subsequent shock expansion phases prevents drawing a straightforward connection between the pre-collapse and post-explosion properties. In order to derive such a connection several explodability studies have been performed in recent years. However, different studies can predict different explosion outcomes. In this article, we show how compactness, which is related to the average density of the star’s core, has an important role in determining the efficiency of neutrino heating, and therefore the outcome of the explosion. Commonly, high-compactness progenitors are assumed to yield failed explosions, due to their large mass accretion rates, preventing the shock from expanding. We show by analysing ~150 2D flash and F ornax simulations and 20 3D F ornax simulations that this is not the case. Instead, due to the rapid increase of neutrino heating with compactness, high-compactness progenitors lead to successful shock revival. We also show that 1D+ simulations that include v-driven convection using a mixing-length theory approach correctly reproduce this trend. Finally, we compare 1D+ models, which we show can reproduce some aspects of multi-D simulations with reasonable accuracy, with other widely used 1D models in the literature.
Relativistic jets are considered to be generated by magnetic fields in a rotating black hole with accretion disk. Consequently, resulting outflows contain magnetic fields in them and control the propagation of jets. We have performed 3D relativistic MHD simulations to investigate the stability and structure of precessed MHD jets with large Lorentz factor by using a newly developed 3D GRMHD code. We have performed simulations of supermagnetosonic jets surrounded by a fast wind. The simulation results reveal complex pressure structure inside the RMHD jet. The structure is produced by a combination of the helical surface and body modes excited by the precession as predicted theoretically. The wavelength of the body mode which occurs in an internal helical twist is much shorter than that of the helical twist surface mode. We will present some comparisons between the RMHD simulations and theoretical predictions, and potential observables and discuss the effect of wind.
Based on our global 3D magnetospheric MHD simulation model, we investigate the phenomena and physical mechanism of the B{sub y} component of the interplanetary magnetic field (IMF) penetrating into the magnetotail. We find that the dayside reconnected magnetic field lines move to the magnetotail, get added to the lobe fields, and are dragged in the IMF direction. However, the B{sub y} component in the plasma sheet mainly originates from the tilt and relative slippage of the south and north lobes caused by plasma convection, which results in the original B{sub z} component in the plasma sheet rotating into a B{sub y} component. Our research also shows that the penetration effect of plasma sheet B{sub y} from the IMF B{sub y} during periods of northward IMF is larger than that during periods of southward IMF.
Purpose: To suggest a new guidance for the evaluation of 3D patient specific QA, a structure-specific risk-index (RI) method was designed and implemented. Methods: A new algorithm was designed to assign the score of Pass, Fail or Pass with Risk to all 3D voxels in each structure by improving a conventional Gamma Index (GI) algorithm, which implied the degree of the risk of under-dose to the treatment target or over-dose to the organ at risks (OAR). Structure-specific distance to agreement (DTOA), dose difference and minimum checkable dose were applied to the GI algorithm, and additional parameters such as dose gradient factor and dose limit of structures were used to the RI method. Maximum passing rate (PR) and minimum PR were designed and calculated for each structure with the RI method. 3D doses were acquired from a spine SBRT plan by simulating the shift of beam iso-center, and tested to show the feasibility of the suggested method. Results: When the iso-center was shifted by 1 mm, 2 mm, and 3 mm, the PR of conventional GI method between shifted and non-shifted 3D doses were 99.9%, 97.4%, and 89.7% for PTV, 99.8%, 84.8%, and 63.2% for spinal cord, and 100%, 99.5%, 91.7% for right lung. The minimum PRs from the RI method were 98.9%, 96.9%, and 89.5% for PTV, and 96.1%, 79.3%, 57.5% for spinal cord, and 92.5%, 92.0%, 84.4% for right lung, respectively. The maximum PRs from the RI method were equal or less than the PRs from the conventional GI evaluation. Conclusion: Designed 3D RI method showed more strict acceptance level than the conventional GI method, especially for OARs. The RI method is expected to give the degrees of risks in the delivered doses, as well as the degrees of agreements between calculated 3D doses and measured (or simulated) 3D doses.
A mixture of n-dodecane and m-xylene is investigated as a diesel fuel surrogate for compression ignition engine applications. Compared to neat n-dodecane, this binary mixture is more representative of diesel fuel because it contains an alkyl-benzene which represents an important chemical class present in diesel fuels. A detailed multi-component mechanism for n-dodecane and m-xylene was developed by combining a previously developed n-dodecane mechanism with a recently developed mechanism for xylenes. The xylene mechanism is shown to reproduce experimental ignition data from a rapid compression machine and shock tube, speciation data from the jet stirred reactor and flame speed data. This combined mechanism was validated by comparing predictions from the model with experimental data for ignition in shock tubes and for reactivity in a flow reactor. The combined mechanism, consisting of 2885 species and 11754 reactions, was reduced to a skeletal mechanism consisting 163 species and 887 reactions for 3D diesel engine simulations. The mechanism reduction was performed using directed relation graph (DRG) with expert knowledge (DRG-X) and DRG-aided sensitivity analysis (DRGASA) at a fixed fuel composition of 77% of n-dodecane and 23% m-xylene by volume. The sample space for the reduction covered pressure of 1 – 80 bar, equivalence ratio of 0.5 – 2.0, and initial temperature of 700 – 1600 K for ignition. The skeletal mechanism was compared with the detailed mechanism for ignition and flow reactor predictions. Finally, the skeletal mechanism was validated against a spray flame dataset under diesel engine conditions documented on the Engine Combustion Network (ECN) website. These multidimensional simulations were performed using a Representative Interactive Flame (RIF) turbulent combustion model. In conclusion, encouraging results were obtained compared to the experiments with regards to the predictions of ignition delay and lift-off length at different ambient temperatures.
The jet composition and energy dissipation mechanism of gamma-ray bursts (GRBs) and blazars are fundamental questions that remain not fully understood. One plausible model is to interpret the γ-ray emission of GRBs and optical emission of blazars as synchrotron radiation of electrons accelerated from the collision-induced magnetic dissipation regions in Poynting-flux-dominated jets. The polarization observation is an important and independent information to test this model. Based on our recent 3D relativistic MHD simulations of collision-induced magnetic dissipation of magnetically dominated blobs, here we perform calculations of the polarization properties of the emission in the dissipation region and apply the results to model the polarization observational data of GRB prompt emission and blazar optical emission. In this article, we show that the same numerical model with different input parameters can reproduce well the observational data of both GRBs and blazars, especially the 90° polarization angle (PA) change in GRB 100826A and the 180° PA swing in blazar 3C279. This supports a unified model for GRB and blazar jets, suggesting that collision-induced magnetic reconnection is a common physical mechanism to power the relativistic jet emission from events with very different black hole masses.
Real clouds and cloud systems are inherently three-dimensional (3D). Because of the limitations in computer resources, however, most cloud-resolving models (CRMs) today are still two-dimensional (2D, see a brief review by Tao 2003). Only recently have 3D experiments been performed for multi-day periods for tropical cloud systems with large horizontal domains at the National Center for Atmospheric Research, at NOAA GFDL, at the U. K. Met. Office, at Colorado State University and at NASA Goddard Space Flight Center (Tao 2003). At Goddard, a 3D Goddard Cumulus Ensemble (GCE) model was used to simulate periods during TOGA COARE (December 19-27, 1992), GATE (September 1-7, 1974), SCSMEX (June 2-11, 1998), ARM (June 26-30, 1997) and KWAJEX (August 7-13, August 18-21, and August 29-September 12, 1999) using a 512 km domain (with 2-kilometer resolution). The results indicate that surface precipitation and latent heating profiles are similar between the 2D and 3D GCE model simulations. However, there are difference in radiation, surface fluxes and precipitation characteristics. The 2D GCE model was used to perform a long-term integration on ARM/GCSS case 4 (22 days at the ARM southern Great Plains site in March 2000). Preliminary results showed a large temperature bias in the upper troposphere that had not been seen in previous tropical cases. The major objectives of this paper are: (1) to determine the sensitivities to model configuration (ie., 2D in west-east, south-north or 3D), (2) to identify the differences and similarities in the organization and entrainment rates of convection between 2D- and 3D-simulated ARM cloud systems, and (3) assess the impact of upper tropospheric forcing on tropical and ARM case 4 cases.
Real clouds and cloud systems are inherently three-dimensional (3D). Because of the limitations in computer resources, however, most cloud-resolving models (CRMs) today are still two-dimensional (2D, see a brief review by Tao 2003). Only recently have 3D experiments been performed for multi-day periods for tropical cloud systems with large horizontal domains at the National Center for Atmospheric Research, at NOAA GFDL, at the U. K. Met. Office, at Colorado State University and at NASA Goddard Space Flight Center (Tao 2003). At Goddard, a 3D Goddard Cumulus Ensemble (GCE) model was used to simulate periods during TOGA COARE (December 19-27, 1992), GATE (September 1-7, 1974), SCSMEX (June 2-11, 1998), ARM (June 26-30, 1997) and KWAJEX (August 7-13, August 18-21, and August 29-September 12, 1999) using a 512 by 512 km domain (with 2-km resolution). The results indicate that surface precipitation and latent heating profiles are similar between the 2D and 3D GCE model simulations. However, there are difference in radiation, surface fluxes and precipitation characteristics. The 2D GCE model was used to perform a long-term integration on ARM/GCSS case 4 (22 days at the ARM Southern Great Plains site in March 2000). Preliminary results showed a large temperature bias in the upper troposphere that had not been seen in previous tropical cases. The major objectives of this paper are: (1) to determine the sensitivities to model configuration (i.e., 2D in west-east, south-north or 3D), (2) to identify the differences and similarities in the organization and entrainment rates of convection between 2D- and 3D-simulated ARM cloud systems, and (3) assess the impact of upper tropospheric forcing on tropical and ARM case 4 cases.