Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation decomposition”

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 37 records · Page 2

Arcjet nozzle area ratio effects

An experimental investigation was conducted to determine the effect of nozzle area ratio on the operating characteristics and performance of a low power dc arcjet thruster. Conical thoriated tungsten nozzle inserts were tested in a modular laboratory arcjet thruster run on hydrogen/nitrogen mixtures simulating the decomposition products of hydrazine. The converging and diverging sides of the inserts had half angles of 30 and 20 degrees, respectively, similar to a flight type unit currently under development. The length of the diverging side was varied to change the area ratio. The nozzle inserts were run over a wide range of specific power. Current, voltage, mass flow rate, and thrust were monitored to provide accurate comparisons between tests. While small differences in performance were observed between the two nozzle inserts, it was determined that for each nozzle insert, arcjet performance improved with increasing nozzle area ratio to the highest area ratio tested and that the losses become very pronounced for area ratios below 50. These trends are somewhat different than those obtained in previous experimental and analytical studies of low Re number nozzles. It appears that arcjet performance can be enhanced via area ratio optimization.

Curran, Francis M.↗

Experimental simulations of oxidizing conditions and organic decomposition on the surface of Mars

One important scientific objective of a Mars Rover Sample Return mission would be to look for traces of living and extinct life on Mars. An instrument to search for organic carbon may be the simplest instrument that could screen samples which are interesting from a biological point of view. An experimental program is described which would help to understand the nature of the oxidizing soil on Mars and the mechanism responsible for organic degradation on the Martian surface. This is approached by lab simulations of the actual conditions that occur on Mars, particularly the oxidant production by atmospheric photochemistry, and the combined effects of UV light and oxidants in decomposing organic compounds. The results will be used to formulate models of the photochemistry of the atmospheric, the atmosphere-soil interaction, and the diffusion of reactive compounds into the soils. This information will provide insights and constraints on the design of a sampling strategy to search for organic compounds on Mars.

Stoker, C. R.↗

Real-time simulation of dynamically deformable finite element models using modal analysis and spectral Lanczos decomposition methods

We present two efficient methods for simulating real-time behavior of a dynamically deformable 3D object modeled by finite element equations. The first method is based on modal analysis, which uses the most significant vibration modes of the object to compute the deformations in real-time for applied forces. The second method uses the spectral Lanczos decomposition to obtain the explicit solutions of the finite element equations.

medical↗

Distributed simulation, no special tools required

A tool kit of C language functions that can be linked with SIMSCRIPT programs to provide the data communication primitives necessary for distributed simulation is presented. A test case is discussed and some timing data are presented. Additionally some metrics, developed to determine the applicability of the server model decomposition for particular simulations, are discussed.

Paterra, Frank↗

Quantifying and Characterizing Striping of Microwave Humidity Sounder with Observation and Simulation

Striping has been observed in the MetOp-A Microwave Humidity Sounder (MHS) data since its degradation in November 2018. However, accurate striping quantification and characterization remain challenging due to the large scene dynamics observed at W/G bands of MHS. Here we have developed a set of novel algorithms for striping quantification, decomposition, characterization, and simulation. Our algorithm extracts striping from the warm-load and cold-space scenes that are relatively stable. We break down the striping into two parts of thermal and 1/f noise and quantify their absolute magnitude and relative ratio. We found a significant increase in striping at 157 GHz that has more than quadrupled by October 2019 relative to its normal level. Regardless of the degradation, the ratio of thermal and 1/f noise remains the same. Our simulation reproduces all the characteristics of striping against observation. It is shown that 1/f noise generates sharp, non-periodic stripes, while thermal noise also generates stripes but with smoother band features. The latter is due to the periodic calibration that has a chopping effect. The striping percentage, defined as the ratio of 1/f to total noise, shows no dependence on the scene temperature. Striping is pronounced not only in 157 GHz but also in 89 and190 GHz with the striping percentage over 50%, while lower in183 GHz of 20%. The results provide insights for quantifying and understanding striping. Our algorithm can be applied to other radiometers and to simulate striping for evaluating its impact in data assimilation and science products.

John Xun Yang↗

Kinematic domain decomposition for boundary-motion-induced flow simulations

A method is developed to solve the unsteady Navier-Stokes equations on a composite grid, which consists of subdomain grids moving with respect to each other. These subdomains are structured grids with different topologies. This method eliminates assuming the moving components to be instantaneously stationary, where deciding on the particular frozen instants is difficult and affects the solution adversely. Moreover, this method captures the boundary-motion-induced flow component. The method is demonstrated through a transonic flow past an airfoil, which experiences a combined motion of pitching and plunging. An O-grid around the airfoil is overlapped on a fine Cartesian grid, which is zonally embedded in a coarse Cartesian grid. The coarse grid is stationary but the other two grids are plunging. Only the O-grid is also sinusoidally pitching. The results are compared successfully with the experimental data.

Baysal, Oktay↗

Marangoni effect and its potential utilization in supporting lunar habitats and other extraterrestrial endeavors

Tears of wine is a phenomenon frequently observed as a ring of wine formed near the top of the glass generates droplets that fall back into the wine. This phenomenon can be explained by the Marangoni effect driven by the surface tension created via gradient on concentration and temperature along with the interface between two phases. Besides being observed while drinking wine, the Marangoni effect is crucial for welding metals, manufacturing integrated circuits, and growing crystals; the Marangoni effect might be also crucial in supporting lunar habitats and other extraterrestrial endeavors as vacuum and reduced gravity are expected to augment the Marangoni effect on extraterrestrial molten soil leading to sustainable extraterrestrial in-situ resources utilization (ISRU). Preliminary Marangoni effect assessment on molten lunar regolith was conducted via JSC-1A lunar regolith simulant under vacuum yielding a spontaneous upwards migration of a uniform molten thin-film front that climbed the crucible wall covering the entire wall’s surface and reaching the crucible’s top end. A temperature gradient within the melt’s bulk and along the crucible’s wall generates a surface tension large enough not only to form a meniscus within the three-phase junction (crucible’s wall, bulk melt, and vacuum) but also to supersede the gravitational force to create the upwards migration of a thin film front originated on the meniscus. The observed extensive wettability of the melt with the crucible surface, an additional key factor besides surface tension necessary to sustain upwards migration, is atypical to what it has been witnessed on molten JSC-1A lunar regolith simulant under non-vacuum conditions and smaller sample sizes. The unique formation and self-migration of a thin-film front could enable key ISRU processes such as thin-film coating, crystal growth, and feedstock generation for additive 3D printing. Further assessment on the higher-temperature vacuum-driven decomposition of the uniform thin-film front as it climbed the crucible’s wall revealed a significant decomposition of various simulant’s metal-oxide components into their respective oxygen and metal elements strengthening the feasibility of extraterrestrial generation and fractional separation of oxygen and metals/alloys with none or minimal need for terrestrial precursors. The Marangoni effect assessment on molten JSC-1A lunar regolith simulant clearly opens novel pathways in future extraterrestrial ISRU endeavors including manufacturing and mining. Extraterrestrial mining will require an operational approach completely different from the terrestrial one that heavily relies on not only mechanical excavation and transportation but also post-processing; an in-situ all-in-one process might be feasible via these two coupled phenomena, self-migration and decomposition.

J A Dominguez↗

Simulations of bypass transition for spatially evolving disturbances

The spatial evolution of disturbances in plane Poiseuille flow and zero pressure gradient boundary layer flow is considered. For disturbances governed by the linearized equations, potential for significant transient growth of the amplitude is demonstrated. The maximum amplification occurs for disturbances with zero or near zero frequencies. Spatial numerical simulations of the transition scenario involving a pair of oblique waves has been conducted for both flows. A fully spectral solver using a simple but efficient fringe region technique allowed the flows to be computed with high resolution into the fully turbulent domain. A modal decomposition of the simulation results indicates that non-linear excitation of the transient growth is responsible for the rapid emergence of low-frequency structures. Physically, this corresponds to streaky flow structures, as seen from the results of a numerical amplitude expansion. Thus, this spatial transition scenario has been found to be similar to the corresponding temporal one. In the boundary layer simulations the streaks are seen to break down from what appears to be a secondary instability.

Lundbladh, A.↗

Implementation and Characterization of Three-Dimensional Particle-in-Cell Codes on Multiple-Instruction-Multiple-Data Massively Parallel Supercomputers

A three-dimensional electrostatic particle-in-cell (PIC) plasma simulation code has been developed on coarse-grain distributed-memory massively parallel computers with message passing communications. Our implementation is the generalization to three-dimensions of the general concurrent particle-in-cell (GCPIC) algorithm. In the GCPIC algorithm, the particle computation is divided among the processors using a domain decomposition of the simulation domain. In a three-dimensional simulation, the domain can be partitioned into one-, two-, or three-dimensional subdomains ("slabs," "rods," or "cubes") and we investigate the efficiency of the parallel implementation of the push for all three choices. The present implementation runs on the Intel Touchstone Delta machine at Caltech; a multiple-instruction-multiple-data (MIMD) parallel computer with 512 nodes. We find that the parallel efficiency of the push is very high, with the ratio of communication to computation time in the range 0.3%-10.0%. The highest efficiency (> 99%) occurs for a large, scaled problem with 64(sup 3) particles per processing node (approximately 134 million particles of 512 nodes) which has a push time of about 250 ns per particle per time step. We have also developed expressions for the timing of the code which are a function of both code parameters (number of grid points, particles, etc.) and machine-dependent parameters (effective FLOP rate, and the effective interprocessor bandwidths for the communication of particles and grid points). These expressions can be used to estimate the performance of scaled problems--including those with inhomogeneous plasmas--to other parallel machines once the machine-dependent parameters are known.

Lyster, P. M.↗

TPSAS-NF1676L-29139-DND

The Arctic has been warming at a rate outpacing globally-averaged warming by 2-3 times, a phenomenon known as "Arctic Amplification" (AA). AA is evident in observations of Arctic surface temperature over the last century as well as in model predictions. The accuracy of model predictions is crucial to determine our response to climate change, as Arctic changes have consequences for the atmospheric circulation, sea level rise, and the carbon cycle. Unfortunately, models show a wide range of possible futures and intermodel spread in projections of Arctic temperature is larger than for any other region: a 2°C global temperature increase- the limit set by the Paris Climate Agreement- would lead to an Arctic temperature change between 3-7°C according to Coupled Model Intercomparison 5 (CMIP5) models. This intermodel spread represents a key deficiency in the understanding of Arctic feedback processes and how changes in various physical parameters (e.g. sea ice cover, clouds, heat transport and surface fluxes) affect the temperature response. Using a process-oriented decomposition of the simulate Arctic surface energy budget, we attribute the largest contribution to the CMIP5 intermodel spead in projected Arctic warming to the Barents-Kara Sea region. Our results reveal that the increasing seasonal fluxing of ocean heat content is a significant part of the intermodel spread in AA and relates strongly to change in surface turbulent fluxes. Our results indicate that if we are to significantly reduce the intermodel spread in projected AA we must focus on understanding the behavior and physical processes that drive surface energy flows in the B-K sea region.

Patrick Taylor↗

Reduced-Order Modeling of a Heaving Airfoil

A reduced-order model of a flapping airfoil is developed using Proper Orthogonal Decomposition (POD). The proper basis functions, developed from snapshots of full Navier-Stokes simulations, are used for a Galerkin projection of the governing equations. The resulting coupled, nonlinear ordinary di.erential equations have a low dimension because the first few basis members capture most of the energy of the flow. The reduced-order model is used to simulate heaving motions that are both similar to and different from the motion(s) used to generate the basis functions, and the errors in the model are quantified. Several methods are used to generate mode sets that can be used over a range of heaving parameters, including snapshots from one, two, and multiple Navier-Stokes simulations. As snapshots from additional simulations are added to the decomposition, the mode sets become richer and can simulate a wider range of parameter space, at some computational cost. Whereas the POD method is fully applicable in three dimensions, the simulation technique based on a body-fixed and body-fitted grid suffers large overhead when extended to three dimensions. To reduce the overhead, an embedding technique is discussed which embeds the solid wing into a fixed Cartesian grid. The wing, which can now have multiple pieces and also be flexible, is represented by a distribution of body forces. This distribution is determined to give exactly the flow around a flapping wing.

Haj-Hariri, H.↗

Counterrotating prop-fan simulations which feature a relative-motion multiblock grid decomposition enabling arbitrary time-steps

Improvements are presented of a computer algorithm developed for the time-accurate flow analysis of rotating machines. The flow model is a finite volume method utilizing a high-resolution approximate Riemann solver for interface flux definitions. The numerical scheme is a block LU implicit iterative-refinement method which possesses apparent unconditional stability. Multiblock composite gridding is used to orderly partition the field into a specified arrangement of blocks exhibiting varying degrees of similarity. Block-block relative motion is achieved using local grid distortion to reduce grid skewness and accommodate arbitrary time step selection. A general high-order numerical scheme is applied to satisfy the geometric conservation law. An even-blade-count counterrotating unducted fan configuration is chosen for a computational study comparing solutions resulting from altering parameters such as time step size and iteration count. The solutions are compared with measured data.

Janus, J. Mark↗

Hydrothermal Decomposition of Amino Acids and Origins of Prebiotic Meteoritic Organic Compounds

The organic compounds found in carbonaceous chondrite meteorites provide insight into primordial solar system chemistry. Evaluating the formation and decomposition mechanisms of meteoritic amino acids may aid our understanding of the origins of life and homochirality on Earth. The amino acid glycine is widespread in meteorites and other extraterrestrial environments; other amino acids, such as isovaline, are found with enantiomeric excesses in some meteorites. The relationship between meteoritic amino acids and other compounds with similar molecular structures, such as aliphatic monoamines and monocarboxylic acids is unclear; experimental results evaluating the decomposition of amino acids have produced inconclusive results about the preferred pathways, reaction intermediates, and if the conditions applied may be compatible with those occurring inside meteoritic parent bodies. In this work, we performed extensive tandem metadynamics, umbrella sampling, and committor analysis to simulate the neutral mild hydrothermal decomposition mechanisms of glycine and isovaline and put them into context for the origins of meteoritic organic compounds. Our ab initio simulations aimed to determine free energy profiles and decomposition pathways for glycine and isovaline. We found that under our modeled conditions, methylammonium, glycolic acid, and sec-butylamine are the most likely decomposition products. These results suggest that meteoritic aliphatic monocarboxylic acids are not produced from decomposition of meteoritic amino acids. Our results also indicate that the decomposition of L-isovaline prefers an enantioselective pathway resulting in the production of (S)-sec-butylamine.

Enhanced sampling↗

Evaluation of the Use of Second Generation Wavelets in the Coherent Vortex Simulation Approach

The objective of this study is to investigate the use of the second generation bi-orthogonal wavelet transform for the field decomposition in the Coherent Vortex Simulation of turbulent flows. The performances of the bi-orthogonal second generation wavelet transform and the orthogonal wavelet transform using Daubechies wavelets with the same number of vanishing moments are compared in a priori tests using a spectral direct numerical simulation (DNS) database of isotropic turbulence fields: 256(exp 3) and 512(exp 3) DNS of forced homogeneous turbulence (Re(sub lambda) = 168) and 256(exp 3) and 512(exp 3) DNS of decaying homogeneous turbulence (Re(sub lambda) = 55). It is found that bi-orthogonal second generation wavelets can be used for coherent vortex extraction. The results of a priori tests indicate that second generation wavelets have better compression and the residual field is closer to Gaussian. However, it was found that the use of second generation wavelets results in an integral length scale for the incoherent part that is larger than that derived from orthogonal wavelets. A way of dealing with this difficulty is suggested.

Goldstein, D. E.↗

The helical decomposition and the instability assumption

Direct numerical simulations show that the triadic transfer function T(k,p,q) peaks sharply when q (or p) is much smaller than k. The triadic transfer function T(k,p,q) gives the rate of energy input into wave number k from all interactions with modes of wave number p and q, where k, p, q form a triangle. This observation was thought to suggest that energy is cascaded downscale through non-local interactions with local transfer and that there was a strong connection between large and small scales. Both suggestions were in contradiction with the classical Kolmogorov picture of the energy cascade. The helical decomposition was found useful in distinguishing between kinematically independent interactions. That analysis has gone beyond the question of non-local interaction with local transfer. In particular, an assumption about the statistical direction of triadic energy transfer in any kinematically independent interaction was introduced (the instability assumption). That assumption is not necessary for the conclusions about non-local interactions with local transfer recalled above. In the case of turbulence under rapid rotation, the instability assumption leads to the prediction that energy is transferred in spectral space from the poles of the rotation axis toward the equator. The instability assumption is thought to be of general validity for any type of triad interactions (e.g. internal waves). The helical decomposition and the instability assumption offer detailed information about the homogeneous statistical dynamics of the Navier-Stokes equations. The objective was to explore the validity of the instability assumption and to study the contributions of the various types of helical interactions to the energy cascade and the subgrid-scale eddy-viscosity. This was done in the context of spectral closures of the Direct Interaction or Quasi-Normal type.

Waleffe, Fabian A.↗

Army-NASA aircrew/aircraft integration program (A3I) software detailed design document, phase 3

The capabilities and design approach of the MIDAS (Man-machine Integration Design and Analysis System) computer-aided engineering (CAE) workstation under development by the Army-NASA Aircrew/Aircraft Integration Program is detailed. This workstation uses graphic, symbolic, and numeric prototyping tools and human performance models as part of an integrated design/analysis environment for crewstation human engineering. Developed incrementally, the requirements and design for Phase 3 (Dec. 1987 to Jun. 1989) are described. Software tools/models developed or significantly modified during this phase included: an interactive 3-D graphic cockpit design editor; multiple-perspective graphic views to observe simulation scenarios; symbolic methods to model the mission decomposition, equipment functions, pilot tasking and loading, as well as control the simulation; a 3-D dynamic anthropometric model; an intermachine communications package; and a training assessment component. These components were successfully used during Phase 3 to demonstrate the complex interactions and human engineering findings involved with a proposed cockpit communications design change in a simulated AH-64A Apache helicopter/mission that maps to empirical data from a similar study and AH-1 Cobra flight test.

Banda, Carolyn↗

Accuracy and resolution of Earth Radiation Budget Estimates

A numerical filter inversion technique that reduces wide-angle satellite measurements to top-of-the-atmosphere radiant exitances has been proposed for the Earth Radiation Budget Experiment (ERBE). The matrix formulation of this technique is presented, and the design of the numerical filter is discussed. The filter is smoothed with a singular value decomposition. The inversion process is simulated by generating synthetic measurements from a 24 degree spherical harmonic radiation field derived from Nimbus 6 ERB data. The numerical filter is applied to these measurements after they are corrupted with instrument error. The results are curves of expected error versus resolution area.

Green, R. N.↗