Search NASA⌕ Search

SEARCH · Search NASA

Results for “Modeling & simulation”

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 865 records · Page 48

SUNBIRD : a simulation-based model for full-shape density-split clustering

Combining galaxy clustering information from regions of different environmental densities can help break cosmological parameter degeneracies and access non-Gaussian information from the density field that is not readily captured by the standard two-point correlation function (2PCF) analyses. However, modelling these density-dependent statistics down to the non-linear regime has so far remained challenging. We present a simulation-based model that is able to capture the cosmological dependence of the full shape of the density-split clustering (DSC) statistics down to intra-halo scales. Our models are based on neural-network emulators that are trained on high-fidelity mock galaxy catalogues within an extended-ΛCDM framework, incorporating the effects of redshift-space, Alcock–Paczynski distortions, and models of the halo–galaxy connection. Our models reach sub-percent level accuracy down to $1 \, h^{-1}\text{Mpc}$ and are robust against different choices of galaxy–halo connection modelling. When combined with the galaxy 2PCF, DSC can tighten the constraints on ω cdm , σ 8 , and n s by factors of 2.9, 1.9, and 2.1, respectively, compared to a 2PCF-only analysis. DSC additionally puts strong constraints on environment-based assembly bias parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

A model for simulating random atmospheres as a function of latitude, season, and time

An empirical stochastic computer model was developed with the capability of generating random thermodynamic profiles of the atmosphere below an altitude of 99 km which are characteristic of any given season, latitude, and time of day. Samples of temperature, density, and pressure profiles generated by the model are statistically similar to measured profiles in a data base of over 6000 rocket and high-altitude atmospheric soundings; that is, means and standard deviations of modeled profiles and their vertical gradients are in close agreement with data. Model-generated samples can be used for Monte Carlo simulations of aircraft or spacecraft trajectories to predict or account for the effects on a vehicle's performance of atmospheric variability. Other potential uses for the model are in simulating pollutant dispersion patterns, variations in sound propagation, and other phenomena which are dependent on atmospheric properties, and in developing data-reduction software for satellite monitoring systems.

Campbell, J. W.↗

Application of system identification to analytic rotor modeling from simulated and wind tunnel dynamic test data, part 2

An introduction to aircraft state and parameter identification methods is presented. A simplified form of the maximum likelihood method is selected to extract analytical aeroelastic rotor models from simulated and dynamic wind tunnel test results for accelerated cyclic pitch stirring excitation. The dynamic inflow characteristics for forward flight conditions from the blade flapping responses without direct inflow measurements were examined. The rotor blades are essentially rigid for inplane bending and for torsion within the frequency range of study, but flexible in out-of-plane bending. Reverse flow effects are considered for high rotor advance ratios. Two inflow models are studied; the first is based on an equivalent blade Lock number, the second is based on a time delayed momentum inflow. In addition to the inflow parameters, basic rotor parameters like the blade natural frequency and the actual blade Lock number are identified together with measurement bias values. The effect of the theoretical dynamic inflow on the rotor eigenvalues is evaluated.

Hohenemser, K. H.↗

Blueprints for Training Information Bottlenecks for Collider Analyses

Dimensionality reduction is a crucial aspect of data analysis in high energy physics, even if accompanied by information loss. Several methods, including histogram- and kernel-based analyses, are only computationally feasible for low-dimensional data. Furthermore, simulation models used in HEP can often only be validated for low-dimensional data. We provide several blueprints for using machine learning to create low-dimensional data representations (continuous event variables and discrete classification labels) for use in signal discovery and parameter estimation tasks. We also describe how to design the learned representation to facilitate a) searches with unknown model parameters and b) validation of simulation models in data control regions.

43 PARTICLE ACCELERATORS↗

A Eulerian-Lagrangian Model to Simulate Two-Phase/Particulate Flows

Figure 1 shows a snapshot of liquid fuel spray coming out of an injector nozzle in a realistic gas-turbine combustor. Here the spray atomization was simulated using a stochastic secondary breakup model (Apte et al. 2003a) with point-particle approximation for the droplets. Very close to the injector, it is observed that the spray density is large and the droplets cannot be treated as point-particles. The volume displaced by the liquid in this region is significant and can alter the gas-phase ow and spray evolution. In order to address this issue, one can compute the dense spray regime by an Eulerian-Lagrangian technique using advanced interface tracking/level-set methods (Sussman et al. 1994; Tryggvason et al. 2001; Herrmann 2003). This, however, is computationally intensive and may not be viable in realistic complex configurations. We therefore plan to develop a methodology based on Eulerian-Lagrangian technique which will allow us to capture the essential features of primary atomization using models to capture interactions between the fluid and droplets and which can be directly applied to the standard atomization models used in practice. The numerical scheme for unstructured grids developed by Mahesh et al. (2003) for incompressible flows is modified to take into account the droplet volume fraction. The numerical framework is directly applicable to realistic combustor geometries. Our main objectives in this work are: Develop a numerical formulation based on Eulerian-Lagrangian techniques with models for interaction terms between the fluid and particles to capture the Kelvin- Helmholtz type instabilities observed during primary atomization. Validate this technique for various two-phase and particulate flows. Assess its applicability to capture primary atomization of liquid jets in conjunction with secondary atomization models.

Apte, S. V.↗

Rapid Prototyping of an Aircraft Model in an Object-Oriented Simulation

A team was created to participate in the Mars Scout Opportunity. Trade studies determined that an aircraft provided the best opportunity to complete the science objectives of the team. A high fidelity six degree of freedom flight simulation was required to provide credible evidence that the aircraft design fulfilled mission objectives and to support the aircraft design process by providing performance evaluations. The team created the simulation using the Langley Standard Real-Time Simulation in C++ (LaSRS++) application framework. A rapid prototyping approach was necessary because the team had only three months to both develop the aircraft simulation model and evaluate aircraft performance as the design and mission parameters matured. The design of LaSRS++ enabled rapid-prototyping in several ways. First, the framework allowed component models to be designed, implemented, unit-tested, and integrated quickly. Next, the framework provides a highly reusable infrastructure that allowed developers to maximize code reuse while concentrating on aircraft and mission specific features. Finally, the framework reduces risk by providing reusable components that allow developers to build a quality product with a compressed testing cycle that relies heavily on unit testing of new components.

Kenney, P. Sean↗

Ground Contact Modeling for the Morpheus Test Vehicle Simulation

The Morpheus vertical test vehicle is an autonomous robotic lander being developed at Johnson Space Center (JSC) to test hazard detection technology. Because the initial ground contact simulation model was not very realistic, it was decided to improve the model without making it too computationally expensive. The first development cycle added capability to define vehicle attachment points (AP) and to keep track of their states in the lander reference frame (LFRAME). These states are used with a spring damper model to compute an AP contact force. The lateral force is then overwritten, if necessary, by the Coulomb static or kinetic friction force. The second development cycle added capability to use the PolySurface class as the contact surface. The class can load CAD data in STL (Stereo Lithography) format, and use the data to compute line of sight (LOS) intercepts. A polygon frame (PFRAME) is computed from the facet intercept normal and used to convert the AP state to PFRAME. Three flat plane tests validate the transitions from kinetic to static, static to kinetic, and vertical impact. The hazardous terrain test will be used to test for visual reasonableness. The improved model is numerically inexpensive, robust, and produces results that are reasonable.

Cordova, Luis↗

Ground Contact Modeling for the Morpheus Test Vehicle Simulation

The Morpheus vertical test vehicle is an autonomous robotic lander being developed at Johnson Space Center (JSC) to test hazard detection technology. Because the initial ground contact simulation model was not very realistic, it was decided to improve the model without making it too computationally expensive. The first development cycle added capability to define vehicle attachment points (AP) and to keep track of their states in the lander reference frame (LFRAME). These states are used with a spring damper model to compute an AP contact force. The lateral force is then overwritten, if necessary, by the Coulomb static or kinetic friction force. The second development cycle added capability to use the PolySurface class as the contact surface. The class can load CAD data in STL (Stereo Lithography) format, and use the data to compute line of sight (LOS) intercepts. A polygon frame (PFRAME) is computed from the facet intercept normal and used to convert the AP state to PFRAME. Three flat plane tests validate the transitions from kinetic to static, static to kinetic, and vertical impact. The hazardous terrain test will be used to test for visual reasonableness. The improved model is numerically inexpensive, robust, and produces results that are reasonable.

Cordova, Luis↗

Changes in Anthropogenic Precursor Emissions Drive Shifts in the Ozone Seasonal Cycle Throughout the Northern Midlatitude Troposphere

Simulations by six Coupled Model Intercomparison Project Phase 6 (CMIP6) Earth system models indicate that the seasonal cycle of baseline tropospheric ozone at northern midlatitudes has been shifting since the mid-20th century. Beginning in ∼ 1940, the magnitude of the seasonal cycle increased by ∼10 ppb (measured from seasonal minimum to maximum), and the seasonal maximum shifted to later in the year by about 3 weeks. This shift maximized in the mid-1980s, followed by a reversal – the seasonal cycle decreased in amplitude and the maximum shifted back to earlier in the year. Similar changes are seen in measurements collected from the 1970s to the present. The timing of the seasonal cycle changes is generally concurrent with the rise and fall of anthropogenic emissions that followed industrialization and the subsequent implementation of air quality emission controls. A quantitative comparison of the temporal changes in the ozone seasonal cycle at sites in both Europe and North America with the temporal changes in ozone precursor emissions across the northern midlatitudes found a high degree of similarity between these two temporal patterns. We hypothesize that changing precursor emissions are responsible for the shift in the ozone seasonal cycle; this is supported by the absence of such seasonal shifts in southern midlatitudes where anthropogenic emissions are much smaller. We also suggest a mechanism by which changing emissions drive the changing seasonal cycle: increasing emissions of NOx allow summertime photochemical production of ozone to become more important than ozone transported from the stratosphere, and increasing volatile organic compounds (VOCs) lead to progressively greater photochemical ozone production in the summer months, thereby increasing the amplitude of the seasonal ozone cycle. Decreasing emissions of both precursor classes then reverse these changes. The quantitative parameter values that characterize the seasonal shifts provide useful benchmarks for evaluating model simulations, both against observations and between models.

Ozone↗

Interannual to Decadal Variability in Subarctic Atlantic Ocean Primary Productivity Historical Variability, Dynamics, and Response to Warming

Phytoplankton primary productivity in the Subarctic Atlantic Ocean sustains valuable fisheries and exports anthropogenic atmospheric carbon dioxide, but multiple Earth system models show that this productivity declines substantially in 21st century global warming scenarios. Prior studies have investigated the spatial and temporal dynamics of the seasonal cycle of the phytoplankton, which exhibit a prominent but spatially variable spring bloom that is caused by the closely related seasonal cycles in solar radiation and vertical mixing that modulate both phytoplankton growth and loss rates. However, recent work indicates that multi-decade declines in this spring bloom with global warming in multiple Earth system models are driven by reductions in circulation of nutrients from the subtropical and tropical thermocline, rather than local changes in vertical mixing or shortwave radiation in the Subarctic. It is unclear if current models of how phytoplankton productivity varies seasonally or over a century in global warming can be interpolated to explain the interannual to decadal variability. In this poster, I will quantify the interannual to decadal variability of phytoplankton chlorophyll using historical satellite ocean color observations and model simulations with the Community Earth System Model. I will then use the observations and model simulations to elucidate the mechanisms of interannual to decadal variability and how that variability is restructured in 21st century global warming scenarios.

nterannual↗

Compressible homogeneous shear: Simulation and modeling

Compressibility effects were studied on turbulence by direct numerical simulation of homogeneous shear flow. A primary observation is that the growth of the turbulent kinetic energy decreases with increasing turbulent Mach number. The sinks provided by compressible dissipation and the pressure dilatation, along with reduced Reynolds shear stress, are shown to contribute to the reduced growth of kinetic energy. Models are proposed for these dilatational terms and verified by direct comparison with the simulations. The differences between the incompressible and compressible fields are brought out by the examination of spectra, statistical moments, and structure of the rate of strain tensor.

Sarkar, S.↗

Compressible homogeneous shear - Simulation and modeling

Compressibility effects were studied on turbulence by direct numerical simulation of homogeneous shear flow. A primary observation is that the growth of the turbulent kinetic energy decreases with increasing turbulent Mach number. The sinks provided by compressible dissipation and the pressure dilatation, along with reduced Reynolds shear stress, are shown to contribute to the reduced growth of kinetic energy. Models are proposed for these dilatational terms and verified by direct comparison with the simulations. The differences between the incompressible and compressible fields are brought out by the examination of spectra, statistical moments, and structure of the rate of strain tensor.

Sarkar, S.↗

All-Atom Modeling and Simulation of Biopolymer Interface: Dual Role of Antifouling Polymer Brushes

Antifouling polymer brushes are well-known for their exceptional resistance to unwanted protein adsorption. While experimental studies have extensively characterized protein–polymer brush interactions, computational investigations remain limited, largely due to the challenges in accurately modeling and integrating these complex, multicomponent systems without resorting to oversimplification. To address this challenge, this study presents one of the most comprehensive and realistic model systems to date, comprising the substrate, grafted polymers, and proteins. Our work interprets the interactions between polymer brush and protein based on realistic modeling without simplification. In particular, this study utilizes molecular modeling and simulation of polycationic and polyzwitterionic brushes─poly(dimethylaminoethyl methacrylate) (PDMAEMA), poly(2-(N-oxide-N,N-dimethylamino)ethyl methacrylate) (PNOMA), and poly(2-(N-3-sulfopropyl-N,N-dimethylammonium)ethyl methacrylate) (PSBMA)─grafted onto α-quartz substrates via polymerization initiator linkers. The brush models were developed to closely replicate experimentally synthesized samples and to provide detailed insights into structural and dynamical changes at the molecular level during protein adsorption. Using steered molecular dynamics simulations, we show that the PSBMA brush, due to its high local density, exhibits the greatest resistance to protein insertion. Cα root-mean-square deviation and interaction pattern analyses further reveal that the PSBMA brush also induces the most significant destabilization of lysozyme, while the PDMAEMA brush enhances protein stability through ion-mediated interactions. The PNOMA brush, while requiring the lowest force for protein adsorption, induces greater protein destabilization than the PDMAEMA brush, primarily due to electrostatic repulsion caused by a short carbon spacer length. Hydration analysis reveals that both the PSBMA brush and the lysozyme interacting with it exhibit the most rapid dehydration, attributed to the brush’s high local chain density, which results in the greatest lysozyme destabilization and the highest adsorption force. These findings highlight the dual role of antifouling polymer brushes: resisting protein adsorption and modulating protein structural dynamics. In conclusion, this study provides valuable insights for the rational design of next-generation antifouling materials and offers a framework for realistic model development in complex multicomponent systems.

Adsorption↗

Mobility analysis, simulation, and scale model testing for the design of wheeled planetary rovers

The use of computer based techniques to model and simulate wheeled rovers on rough natural terrains is considered. Physical models of a prototype vehicle can be used to test the correlation of the simulations in scaled testing. The computer approaches include a quasi-static planar or two dimensional analysis and design tool based on the traction necessary for the vehicle to have imminent mobility. The computer program modeled a six by six wheel drive vehicle of original kinematic configuration, called the Rocker Bogie. The Rocker Bogie was optimized using the quasi-static software with respect to its articulation parameters prior to fabrication of a prototype. In another approach used, the dynamics of the Rocker Bogie vehicle in 3-D space was modeled on an engineering workstation using commercial software. The model included the complex and nonlinear interaction of the tire and terrain. The results of the investigation yielded numerical and graphical results of the rover traversing rough terrain on the earth, moon, and Mars. In addition, animations of the rover excursions were also generated. A prototype vehicle was then used in a series of testbed and field experiments. Correspondence was then established between the computer models and the physical model. The results indicated the utility of the quasi-static tool for configurational design, as well as the predictive ability of the 3-D simulation to model the dynamic behavior of the vehicle over short traverses.

Lindemann, Randel A.↗

Vertical Profiles of Latent Heat Release and Their Retrieval for TOGA COARE Convective Systems Using a Cloud Resolving Model, SSM/I, and Ship-borne Radar Data

Latent heating profiles associated with three (TOGA COARE) Tropical Ocean and Global Atmosphere Coupled Ocean Atmosphere Response Experiment active convective episodes (December 10-17 1992; December 19-27 1992; and February 9-13 1993) are examined using the Goddard Cumulus Ensemble (GCE) Model and retrieved by using the Goddard Convective and Stratiform Heating (CSH) algorithm . The following sources of rainfall information are input into the CSH algorithm: Special Sensor Microwave Imager (SSM/1), Radar and the GCE model. Diagnostically determined latent heating profiles calculated using 6 hourly soundings are used for validation. The GCE model simulated rainfall and latent heating profiles are in excellent agreement with those estimated by soundings. In addition, the typical convective and stratiform heating structures (or shapes) are well captured by the GCE model. Radar measured rainfall is smaller than that both estimated by the GCE model and SSM/I in all three different COARE IFA periods. SSM/I derived rainfall is more than the GCE model simulated for the December 19-27 and February 9-13 periods, but is in excellent agreement with the GCE model for the December 10-17 period. The GCE model estimated stratiform amount is about 50% for December 19-27, 42% for December 11-17 and 56% for the February 9-13 case. These results are consistent with large-scale analyses. The accurate estimates of stratiform amount is needed for good latent heating retrieval. A higher (lower) percentage of stratiform rain can imply a maximum heating rate at a higher (lower) altitude. The GCE model always simulates more stratiform rain (10 to 20%) than the radar for all three convective episodes. SSM/I derived stratiform amount is about 37% for December 19-27, 48% for December 11-17 and 41% for the February 9-13 case. Temporal variability of CSH algorithm retrieved latent heating profiles using either GCE model simulated or radar estimated rainfall and stratiform amount is in good agreement with that diagnostically determined for all three periods. However, less rainfall and a smaller stratiform percentage estimated by radar resulted in a weaker (underestimated) latent heating profile and a lower maximum latent heating level compared to those determined diagnostically. Rainfall information from SSM/I can not retrieve individual convective events due to poor temporal sampling. Nevertheless, this study suggests that a good 4r, rainfall retrieval from SSM/I for a convective event always leads to a good latent heating retrieval. Sensitivity testing has been performed and the results indicate that the SSM/I derived time averaged stratiform amount may be underestimated for December 19-27. Time averaged heating profiles derived from SSM/I, however, are not in bad agreement with those derived by soundings for the December 10-17 convective period. The heating retrievals may be more accurate for longer time scales provided there is no bias in the sampling.

Tao, Wei-Kuo↗

An Integrated Approach to Modeling Solar Electric Propulsion Vehicles During Long Duration, Near-Earth Orbit Transfers

Recent NASA interest in utilizing solar electronic propulsion (SEP) technology to transfer payloads, e.g. from low-Earth orbit (LEO) to higher energy geostationary-Earth orbit (GEO) or to Earth escape, has necessitated the development of high fidelity SEP vehicle models and simulations. These models and simulations need to be capable of capturing vehicle dynamics and sub-system interactions experienced during the transfer trajectories which are typically accomplished with continuous-burn (potentially interrupted by solar eclipse), long duration "spiral out" maneuvers taking several months or more to complete. This paper presents details of an integrated simulation approach achieved by combining a high fidelity vehicle simulation code with a detailed solar array model. The combined simulation tool gives researchers the functionality to study the integrated effects of various vehicle sub-systems (e.g. vehicle guidance, navigation and control (GN&C), electric propulsion system (EP)) with time varying power production. Results from a simulation model of a vehicle with a 50 kW class SEP system using the integrated tool are presented and compared to the results from another simulation model employing a 50 kW end-of-life (EOL) fixed power level assumption. These models simulate a vehicle under three degree of freedom dynamics (i.e. translational dynamics only) and include the effects of a targeting guidance algorithm (providing a "near optimal" transfer) during a LEO to near Earth escape (C (sub 3) = −2.0 km (sup 2) / sec (sup −2) spiral trajectory. The presented results include the impact of the fully integrated, time-varying solar array model (e.g. cumulative array degradation from traversing the Van Allen belts, impact of solar eclipses on the vehicle and the related temperature responses in the solar arrays due to operating in the Earth's thermal environment, high fidelity array power module, etc.); these are used to assess the impact on vehicle performance (i.e. propellant consumption) and transit times.

vehicle↗