Search NASA⌕ Search

SEARCH · Search NASA

Results for “Semi-analytical method”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Optimization of Ocean Color Algorithms: Application to Satellite Data Merging

The objective of the program is to develop and validate a procedure for ocean color data merging, which is one of the major goals of the SIMBIOS project. As part of the SIMBIOS Program, we have developed a merging method for ocean color data. Conversely to other methods our approach does not combine end-products like the subsurface chlorophyll concentration (chl) from different sensors to generate a unified product. Instead, our procedure uses the normalized water-leaving radiances L((sub wN)(lambda)) from single or multiple sensors and uses them in the inversion of a semi-analytical ocean color model that allows the retrieval of several ocean color variables simultaneously. Beside ensuring simultaneity and consistency of the retrievals (all products are derived from a single algorithm), this model-based approach has various benefits over techniques that blend end-products (e.g. chlorophyll): 1) It works with single or multiple data sources regardless of their specific bands; 2) It exploits band redundancies and band differences; 3) It accounts for uncertainties in the L((sub wN)(lambda)) data; 4) It provides uncertainty estimates for the retrieved variables.

Ritorena, Stephane↗
NASA NTRS1999

Electric organ discharges and electric images during electrolocation

Weakly electric fish use active electrolocation - the generation and detection of electric currents - to explore their surroundings. Although electrosensory systems include some of the most extensively understood circuits in the vertebrate central nervous system, relatively little is known quantitatively about how fish electrolocate objects. We believe a prerequisite to understanding electrolocation and its underlying neural substrates is to quantify and visualize the peripheral electrosensory information measured by the electroreceptors. We have therefore focused on reconstructing both the electric organ discharges (EODs) and the electric images resulting from nearby objects and the fish's exploratory behaviors. Here, we review results from a combination of techniques, including field measurements, numerical and semi-analytical simulations, and video imaging of behaviors. EOD maps are presented and interpreted for six gymnotiform species. They reveal diverse electric field patterns that have significant implications for both the electrosensory and electromotor systems. Our simulations generated predictions of the electric images from nearby objects as well as sequences of electric images during exploratory behaviors. These methods are leading to the identification of image features and computational algorithms that could reliably encode electrosensory information and may help guide electrophysiological experiments exploring the neural basis of electrolocation.

Review, Tutorial↗
NASA NTRS1997

CGRO Guest Investigator Program

The following are highlights from the research supported by this grant: (1) Theory of gamma-ray blazars: We studied the theory of gamma-ray blazars, being among the first investigators to propose that the GeV emission arises from Comptonization of diffuse radiation surrounding the jet, rather than from the synchrotron-self-Compton mechanism. In related work, we uncovered possible connections between the mechanisms of gamma-ray blazars and those of intraday radio variability, and have conducted a general study of the role of Compton radiation drag on the dynamics of relativistic jets. (2) A Nonlinear Monte Carlo code for gamma-ray spectrum formation: We developed, tested, and applied the first Nonlinear Monte Carlo (NLMC) code for simulating gamma-ray production and transfer under much more general (and realistic) conditions than are accessible with other techniques. The present version of the code is designed to simulate conditions thought to be present in active galactic nuclei and certain types of X-ray binaries, and includes the physics needed to model thermal and nonthermal electron-positron pair cascades. Unlike traditional Monte-Carlo techniques, our method can accurately handle highly non-linear systems in which the radiation and particle backgrounds must be determined self-consistently and in which the particle energies span many orders of magnitude. Unlike models based on kinetic equations, our code can handle arbitrary source geometries and relativistic kinematic effects In its first important application following testing, we showed that popular semi-analytic accretion disk corona models for Seyfert spectra are seriously in error, and demonstrated how the spectra can be simulated if the disk is sparsely covered by localized 'flares'.

Begelman, Mitchell C.↗
DOE OSTI2025

A semi-analytic estimate for the effective sound speed counterterm in the EFTofLSS

The Effective Field Theory of Large Scale Structure (EFTofLSS) has found tremendous success as a perturbative framework for the evolution of large scale structure, and it is now routinely used to compare theoretical predictions against cosmological observations. The model for the total matter field includes one nuisance parameter at 1-loop order, the effective sound speed, which can be extracted by matching the EFT to full N-body simulations. In this work we first leverage the Layzer-Irvine cosmic energy equation to show that the equation of state can be exactly computed with knowledge of the fully nonlinear power spectrum. When augmented with separate universe methods, we show one can estimate the effective sound speed. This estimate is in good agreement with simulation results, with errors at the few tens of percent level. Here, we apply our method to investigate the cosmology dependence of the effective sound speed and to shed light on what cosmic structures shape its value.

Cosmological perturbation theory in GR and beyond↗
NASA NTRS1999

Entry Guidance for the Reusable Launch Vehicle

The X-33 Advanced Technology Demonstrator is a half-scale prototype developed to test the key technologies needed for a full-scale single-stage reusable launch vehicle (RLV). The X-33 is a suborbital vehicle that will be launched vertically, and land horizontally. The goals of this research were to develop an alternate entry guidance scheme for the X-33 in parallel to the actual X-33 entry guidance algorithms, provide comparative and complementary study, and identify potential new ways to improve entry guidance performance. Toward these goals, the nominal entry trajectory is defined by a piecewise linear drag-acceleration-versus-energy profile, which is in turn obtained by the solution of a semi-analytical parameter optimization problem. The closed-loop guidance is accomplished by tracking the nominal drag profile with primarily bank-angle modulation on-board. The bank-angle is commanded by a single full-envelope nonlinear trajectory control law. Near the end of the entry flight, the guidance logic is switched to heading control in order to meet strict conditions at the terminal area energy management interface. Two methods, one on ground-track control and the other on heading control, were proposed and examined for this phase of entry guidance where lateral control is emphasized. Trajectory dispersion studies were performed to evaluate the effectiveness of the entry guidance algorithms against a number of uncertainties including those in propulsion system, atmospheric properties, winds, aerodynamics, and propellant loading. Finally, a new trajectory-regulation method is introduced at the end as a promising precision entry guidance method. The guidance principle is very different and preliminary application in X-33 entry guidance simulation showed high precision that is difficult to achieve by existing methods.

Lu, Ping↗
NASA NTRS2008

Black Hole Mergers as Probes of Structure Formation

Observations of gravitational waves from massive black hole (MBH) mergers can provide us with important clues about the era of structure formation in the early universe. Previous research in this field has been limited to calculating merger rates of MBHs using different models where many assumptions are made about the specific values of physical parameters of the mergers, resulting in merger rate estimates that span 5 to 6 orders of magnitude. We develop a semi-analytical, phenomenological model that includes plausible combinations of several physical parameters involved in the mergers. which we then turn around to determine how well LISA observations will be able to enhance our understanding of the universe during the critical z approximately equal to 5-30 structure formation era. We do this by generating synthetic LISA observable data (masses, redshifts, merger rates), which are then analyzed using a Markov Chain Monte Carlo (MCMC) method. This allows us to constrain the physical parameters of the mergers.

Alicea-Munoz, Emily↗
NASA NTRS2008

Black Hole Mergers as Probes of Structure Formation

Intense structure formation and reionization occur at high redshift, yet there is currently little observational information about this very important epoch. Observations of gravitational waves from massive black hole (MBH) mergers can provide us with important clues about the formation of structures in the early universe. Past efforts have been limited to calculating merger rates using different models in which many assumptions are made about the specific values of physical parameters of the mergers, resulting in merger rate estimates that span a very wide range (0.1 - 104 mergers/year). Here we develop a semi-analytical, phenomenological model of MBH mergers that includes plausible combinations of several physical parameters, which we then turn around to determine how well observations with the Laser Interferometer Space Antenna (LISA) will be able to enhance our understanding of the universe during the critical z 5 - 30 structure formation era. We do this by generating synthetic LISA observable data (total BH mass, BH mass ratio, redshift, merger rates), which are then analyzed using a Markov Chain Monte Carlo method. This allows us to constrain the physical parameters of the mergers. We find that our methodology works well at estimating merger parameters, consistently giving results within 1- of the input parameter values. We also discover that the number of merger events is a key discriminant among models. This helps our method be robust against observational uncertainties. Our approach, which at this stage constitutes a proof of principle, can be readily extended to physical models and to more general problems in cosmology and gravitational wave astrophysics.

Alicea-Munoz, E.↗
NASA NTRS2008

Using Black Hole Mergers to Explore Structure Formation

Observations of gravitational waves from massive black hole mergers will open a new window into the era of structure formation in the early universe. Past efforts have concentrated on calculating merger rates using different physical assumptions, resulting in merger rate estimates that span a wide range (0.1 - 1 0A4 mergers/year). We develop a semi-analytical, phenomenological model of massive black hole mergers that includes plausible combinations of several physical parameters, which we then turn around to determine how well observations with the Laser Interferometer Space Antenna (LISA) will be able to enhance our understanding of the universe during the critical z approx. 5 - 30 epoch. Our approach involves generating synthetic LISA observable data (total BH masses, BH mass ratios, redshifts, merger rates), which are then analyzed using a Markov Chain Monte Carlo method, thus finding constraints for the physical parameters of the mergers. We find that our method works well at estimating merger parameters and that the number of merger events is a key discriminant among models, therefore making our method robust against observational uncertainties. Our approach can also be extended to more physically-driven models and more general problems in cosmology.

Alicea-Munoz, E.↗
NASA NTRS2009

Using Black Hole Mergers to Explore Structure Formation

Observations of gravitational waves from massive black hole mergers will open a new window into the era of structure formation in the early universe. Past efforts have concentrated on calculating merger rates using different physical assumptions, resulting in merger rate estimates that span a wide range (0.1 - 10(exp 4) mergers/year). We develop a semi-analytical, phenomenological model of massive black hole mergers that includes plausible combinations of several physical parameters, which we then turn around to determine how well observations with the Laser Interferometer Space Antenna (LISA) will be able to enhance our understanding of the universe during the critical z approximately equal to 5-30 epoch. Our approach involves generating synthetic LISA observable data (total BH masses, BH mass ratios, redshifts, merger rates), which are then analyzed using a Markov Chain Monte Carlo method, thus finding constraints for the physical parameters of the mergers. We find that our method works well at estimating merger parameters and that the number of merger events is a key discriminant among models, therefore making our method robust against observational uncertainties. Our approach can also be extended to more physically-driven models and more general problems in cosmology. This work is supported in part by the Cooperative Education Program at NASA/GSFC.

Alicea-Munoz, E.↗
NASA NTRS1992

Elastic Analysis of Three-Dimensional Solids with Fiber Inclusions by Bem

A very efficient method of analysis for the elastic behavior of three-dimensional solids with a large number of embedded fiber inclusions has been developed. For specific applications to the analysis of cylindrical shaped fiber inclusions significant gain in efficiency can be achieved by defining these inclusions by a system of curvilinear line elements with a prescribed diameter and by assuming a variation in the traction and displacement field in the circumferential directions in terms of a trigonometric shape function together with a linear or quadratic variation along the longitudinal direction. The resulting integrals are then treated semi-analytically. A number of examples of the analysis which has been developed for composite elements are described.

Banerjee, P. K.↗
NASA NTRS2012

Launch Window Analysis for the Magnetospheric Multiscale Mission

The NASA Magnetospheric Multiscale (MMS) mission will fly four spinning spacecraft in formation in highly elliptical orbits to study the magnetosphere of the Earth. This paper describes the development of an MMS launch window tool that uses the orbitaveraged Variation of Parameter equations as the basis for a semi-analytic quantification of the dominant oblateness and lunisolar perturbation effects on the MMS orbit. This approach, coupled with a geometric interpretation of all of the MMS science and engineering constraints, allows a scan of 180(sup 2) = 32,400 different (RAAN, AOP) pairs to be carried out for a specified launch day in less than 10 s on a typical modern laptop. The resulting plot indicates the regions in (RAAN, AOP) space where each constraint is satisfied or violated: their intersection gives, in an easily interpreted graphical manner, the final solution space for the day considered. This tool, SWM76, is now used to provide launch conditions to the full fidelity (but far slower) MMS simulation code: very good agreement has been observed between the two methods.

MAGNETOSPHERIC↗
DOE OSTI2025

Semi-analytical covariance matrices for two-point correlation function for DESI 2024 data

We present an optimized way of producing the fast semi-analytical covariance matrices for the Legendre moments of the two-point correlation function, taking into account survey geometry and mimicking the non-Gaussian effects. We validate the approach on simulated (mock) catalogs for different galaxy types, representative of the Dark Energy Spectroscopic Instrument (DESI) Data Release 1, used in 2024 analyses. We find only a few percent differences between the mock sample covariance matrix and our results, which can be expected given the approximate nature of the mocks, although we do identify discrepancies between the shot-noise properties of the DESI fiber assignment algorithm and the faster approximation (emulator) used in the mocks. Importantly, we find a close agreement (≤ 8% relative differences) in the projected errorbars for distance scale parameters for the baryon acoustic oscillation measurements. This confirms our method as an attractive alternative to simulation-based covariance matrices, especially for non-standard models or galaxy sample selections, making it particularly relevant to the broad current and future analyses of DESI data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Myna: Connecting powder bed fusion build data to simulation tools for digital twin applications

Additive manufacturing (AM), as a digital process, can generate a detailed digital thread linking a part’s design and manufacturing to its operational performance. As AM systems advance, an increasing amount of process data is stored in manufacturing databases. In principle, this data can be utilized by simulation-based digital twin approaches, such as real-time process control and asynchronous post-processing guidance. However, few tools currently exist for systematically integrating digital thread data with computational tools. Here, in this study, we propose a software package, called Myna, for connecting data from powder bed fusion processes to simulation tools. The utility of such a platform is demonstrated using build data from the Oak Ridge National Laboratory Manufacturing Demonstration Facility “Peregrine v2023-10” public dataset to automatically configure and run 54 semi-analytical 3DThesis melt pool simulations, 78 numerical Additive FOAM melt pool simulations, and 3 ExaCA microstructure simulations. The simulated, spatially registered microstructures are then compared directly with electron backscatter diffraction characterization of the corresponding as-built part locations. The resulting simulated microstructure showed variation as a function of process parameters, particularly stripe width; however, the experimental data had little variation between the microstructure texture and grain size resulting from different processing conditions. Analysis of the discrepancies suggest that it is possible a two-phase ferritic-austenitic solidification model is needed to accurately predict grain size and texture for certain stainless steel 316L feedstock compositions under powder bed fusion conditions, providing direction for future research. As illustrated here, due to the number and complexity of the simulations involved in AM process-structure–property predictions, automated methods to connect process data and simulations will remain necessary tools for testing hypotheses and implementing digital twin applications.

Knapp, Gerald L. [Oak Ridge National Laboratory (O↗
DOE OSTI2026

Semi-Analytical Hierarchical Bayesian Inference of Nonlinear Model Structure in Stochastic Dynamics: Applied to Compartmental Models of Infectious Diseases

A Bayesian computational framework for parsimonious inference in stochastic nonlinear dynamical systems is presented. This framework enables the concurrent estimation of system states, time-varying parameters, time-invariant parameters, and the optimal sparsity structure of the model parameters. Because differential equation-based models are often simplified mechanistic or phenomenological representations, robust inference from noisy measurement data requires explicit treatment of model error and uncertainty. Model error and time-varying parameters can be represented as random processes, enabling inference while making minimal assumptions about the underlying sources of discrepancy and variability. Adopting stochastic differential equation representations affords the model significant flexibility, but can also render it susceptible to overfitting during statistical inversion, where the inferred model may track noise rather than the underlying signal. To alleviate the effects of overfitting and to enable the discovery of the optimal sparse representation of the time-invariant parameters, a Bayesian sparse learning algorithm is embedded within the framework. This sparse learning framework adopts an approximate hierarchical Bayesian setting defined by a series of semi-analytical expressions. The model structure inference framework is validated using a stochastic compartmental model for tracking and forecasting active cases of an infectious disease. Compartmental models describe population-level infectious disease dynamics through interactions among population fractions grouped by disease state. Mathematically, such models consist of a system of coupled ordinary differential equations. This example adopts an expressive compartmental model that includes multiple possible interactions between disease states, motivated by early uncertainty surrounding COVID-19 reinfection dynamics and their implications for long-term epidemic forecasting. The sparse learning exercise permits the inference of a priori unknown epidemiological dynamics from simulated public health data, discovering the nested compartmental model that optimizes the trade-off between average data-fit and model complexity. It is shown that inducing sparsity among the model parameters eliminates redundant interactions between compartments, equivalently revealing the optimal coupling structure between differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Probabilistic inference of the structure and orbit of Milky Way satellites with semi-analytic modelling

Semi-analytic modelling furnishes an efficient avenue for characterizing dark matter haloes associated with satellites of Milky Way-like systems, as it easily accounts for uncertainties arising from halo-to-halo variance, the orbital disruption of satellites, baryonic feedback, and the stellar-to-halo mass (SMHM) relation. We use the SatGen semi-analytic satellite generator, which incorporates both empirical models of the galaxy–halo connection as well as analytic prescriptions for the orbital evolution of these satellites after accretion onto a host to create large samples of Milky Way-like systems and their satellites. By selecting satellites in the sample that match observed properties of a particular dwarf galaxy, we can infer arbitrary properties of the satellite galaxy within the cold dark matter paradigm. For the Milky Way’s classical dwarfs, we provide inferred values (with associated uncertainties) for the maximum circular velocity v max and the radius r max at which it occurs, varying over two choices of baryonic feedback model and two prescriptions for the SMHM relation. While simple empirical scaling relations can recover the median inferred value for v max and r max , this approach provides realistic correlated uncertainties and aids interpretability. We also demonstrate how the internal properties of a satellite’s dark matter profile correlate with its orbit, and we show that it is difficult to reproduce observations of the Fornax dwarf without strong baryonic feedback. Furthermore, the technique developed in this work is flexible in its application of observational data and can leverage arbitrary information about the satellite galaxies to make inferences about their dark matter haloes and population statistics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1968

Internal low Reynolds number turbulent heat transfer

The results of a semi-analytical and experimental study of adiabatic tube flow and an analytical and experimental investigation of the thermal entrance region for gas flowing through electrically heated circular tubes are presented. Emphasis is placed on the low Reynolds number turbulent flow regime--defined as fully turbulent flow at bulk Reynolds numbers from 3,000 to about 15,000. Adiabatic air velocity and friction data and localized heat transfer measurements for air and helium, at low heating rates, are presented for this range. The adiabatic data were obtained in a 1.61 inch 1D tube for flow at bulk Reynolds numbers from 3,000 to 15,000. A continuous, Reynolds number-dependent, profile is developed from the data by using a modification of Reichardts' wall and middle law eddy diffusivity expressions. The velocity profile satisfies continuity. It is valid for all Reynolds numbers in excess of 3,000 for which the flow is fully turbulent and the Blasius friction factor expression is valid. The thermal entrance problem for a fully developed velocity profile is solved analytically by the method of Sparrow, Hallman, and Siegel. The solution is based on the profile developed from the velocity study. Tabular values of the eigenvalues and normalized Nusselt numbers for gases are presented for it range of Reynolds numbers from 3,000 to 50,000. The axial variation of Nusselt number is found to be correlated by [Nu/Nu_(∞)] = 1 + 0.8 (1+ 70,000 Re^(-3/2)) ((x/D)^(-1)) to within ± 5 per cent for x/D ≥2. The fully developed value agrees with the Dittus - Boelter correlation, Nu_(∞) = 0.021 Re^(0.8) Pr^(0.4) For the eigenvalues, λ^(2)_(n), and the associated constants, A_(n), correlations of the form λ^(2)_(n) =A_(1,n) Re^(-b_(1,n) + C_(1,n) Re^(-d_(1,n)) A_(n) = -A_(2,n) Re^(-b_(2,n) + C_(2,n) Re^(-d_(2,n)) are obtained. The coefficients and powers are presented in tabular form. Heat transfer data are presented, primarily for helium, for the low Reynolds number turbulent range. A one-quarter inch, resistively heated, vertical, circular tube was used for the study. The data cover an axial range from 1.2 to 96 diameters; wall-to-bulk temperature ratios vary from 1 to 1.4. In the low Reynolds number turbulent regime, these data clearly support the present analytical solution rather than the prediction obtained by applying the eddy diffusivity distribution used by Sparrow, Hallman, and Siegel.

Harold C. Reynolds, Jr.↗
NASA NTRS1990

A numerical method for phase-change problems

A highly accurate and efficient finite-difference method for phase-change problems with multiple moving boundaries of irregular shape is developed by employing a coordinate transformation that immobilizes moving boundaries and preserves the conservative forms of the original governing equations. The numerical method is first presented for one-dimensional phase-change problems (involving large density variation between phases, heat generation, and multiple moving boundaries) and then extended to solve two-dimensional problems (without change of densities between phases). Numerical solutions are obtained non-iteratively using an explicit treatment of the interfacial mass and energy balances and an implicit treatment of the temperature field equations. The accuracy and flexibility of the present numerical method are verified by solving some phase-change problems and comparing the results with existing analytical, semi-analytical and numerical solutions. Results indicate that one- and two-dimensional phase-change problems can be handled easily with excellent accuracies.

Kim, Charn-Jung↗
NASA NTRS2000

Efficient Simulation of Wing Modal Response: Application of 2nd Order Shape Sensitivities and Neural Networks

At the preliminary design stage of a wing structure, an efficient simulation, one needing little computation but yielding adequately accurate results for various response quantities, is essential in the search of optimal design in a vast design space. In the present paper, methods of using sensitivities up to 2nd order, and direct application of neural networks are explored. The example problem is how to decide the natural frequencies of a wing given the shape variables of the structure. It is shown that when sensitivities cannot be obtained analytically, the finite difference approach is usually more reliable than a semi-analytical approach provided an appropriate step size is used. The use of second order sensitivities is proved of being able to yield much better results than the case where only the first order sensitivities are used. When neural networks are trained to relate the wing natural frequencies to the shape variables, a negligible computation effort is needed to accurately determine the natural frequencies of a new design.

Kapania, Rakesh K.↗