Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Generalized additive models”

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

Regularization by denoising diffusion models for solving inverse PDE problems with application to full waveform inversion

Partial differential equation (PDE)-governed inverse problems are fundamental across various scientific and engineering applications; yet they face significant challenges due to nonlinearity, ill-posedness, and sensitivity to noise. Here, we introduce a computational framework, regularization by denoising using diffusion models for partial differential equations (RED-DiffEq), by integrating physics-driven inversion and data-driven learning. RED-DiffEq leverages pretrained diffusion models as a regularization mechanism for PDE-governed inverse problems. We apply RED-DiffEq to solve the full waveform inversion problem in geophysics, a challenging seismic imaging technique that seeks to reconstruct high-resolution subsurface velocity models from seismic measurement data. Our method shows enhanced accuracy and robustness compared to benchmark methods. Additionally, it exhibits strong generalization and domain decomposition capacity, enabling the inversion of more complex velocity models with larger domains than those used in training the diffusion model. Our framework can also be directly applied to diverse PDE-governed inverse problems.

Shan, Siming [Yale University, New Haven, CT (Unit

Models, In situ, and Remote sensing of Aerosols (MIRA): Formation of an International Working Group

There is a natural partitioning of scientific interest amongst three specialties of aerosol research: modeling, in situ measurements, and remote sensing. The community sees enhanced measurement capabilities when these groups interact, and this strengthens the overall scientific impact on climate and air quality. The Models, In situ, and Remote sensing of Aerosols (MIRA) working group connects members of the different aerosol communities through collaborative projects. What is MIRA? MIRA is forum that fosters international collaborations amongst the aerosol specialties. MIRA is also a collection of interdisciplinary projects with clear goals that are pursued by small working groups. Finally, MIRA projects are generally characterized by requests for additional scientific data (both observational and modeled). Why? The purpose of MIRA is to contextualize both observations and model results through the encouragement of holistic projects and collaborations. How does MIRA differ from other projects? MIRA focuses on interdisciplinarity to improve measurements and their utility, so MIRA complements the activities of other groups. For example, ensemble model runs of AeroCom could be used in a MIRA project with greater robustness than a similar effort that uses single-model analyses. We first present a description of MIRA and how your project can become a part of the MIRA community. This is followed by brief descriptions of some of the current MIRA projects, which include the Mapping of Aerosol lidar ratios for CALIPSO (MAC), the Tables of Aerosol Optics (TAO), and the Harmonization of aerosol Assimilation Models and Retrievals (HAMR). Finally, we discuss the immediate science goals and organization of MIRA.

Gregory Schuster

Models, In situ, and Remote sensing of Aerosols (MIRA)

There is a natural partitioning of scientific interest amongst three specialties of aerosol research: modeling, in situ measurements, and remote sensing. The community sees enhanced measurement capabilities when these groups interact, and this strengthens the overall scientific impact on climate and air quality. The Models, In situ, and Remote sensing of Aerosols (MIRA) working group connects members of the different aerosol communities through collaborative projects. What is MIRA? MIRA is a forum that fosters international collaborations amongst the aerosol specialties. MIRA is also a collection of interdisciplinary projects with clear goals that are pursued by small working groups. Finally, MIRA projects are generally characterized by requests for additional scientific data (both observational and modeled). Why? The purpose of MIRA is to contextualize both observations and model results through the encouragement of holistic projects and collaborations. How does MIRA differ from other projects? MIRA focuses on interdisciplinarity to improve measurements and their utility, so MIRA complements the activities of other groups. For example, ensemble model runs of AeroCom could be used in a MIRA project with greater robustness than a similar effort that uses single-model analyses. We present a description of MIRA and brief descriptions of some of the current MIRA topics, which include Satellite-Assisted Particulate Matter (SAPM), Mapping of Aerosol lidar ratios for CALIPSO (MAC), Tables of Aerosol Optics (TAO), and the Harmonization of aerosol Assimilation Models and Retrievals (HAMR). Finally, we discuss the immediate science goals and organization of MIRA. See https://science.larc.nasa.gov/mira-wg/ for more information.

Travis Toth

Program For Simulating Dynamics Of Aerospace Vehicles

SORT (Simulation and Optimization of Rocket Trajectories) is general-purpose three-degree-of-freedom with three axis static moment balance simulation of flight dynamics of arbitrary aerospace vehicle. Modular structure facilitates application to variety of trajectory-analysis problems. Contains math model of aerodynamics completely generalized. Computes both longitudinal and lateral forces and moments. In addition to fore-body coefficients, computes longitudinal base effect aerodynamic forces and moments. Simplified ballistic-coefficient model also available for analysis of ballistic entry. Written using ANSI FORTRAN 77.

Berning, M. J.

Refinement of a limit cycle oscillator model of the effects of light on the human circadian pacemaker

In 1990, Kronauer proposed a mathematical model of the effects of light on the human circadian pacemaker. Although this model predicted many general features of the response of the human circadian pacemaker to light exposure, additional data now available enable us to refine the original model. We first refined the original model by incorporating the results of a dose response curve to light into the model's predicted relationship between light intensity and the strength of the drive onto the pacemaker. Data from three bright light phase resetting experiments were then used to refine the amplitude recovery characteristics of the model. Finally, the model was tested and further refined using data from an extensive phase resetting experiment in which a 3-cycle bright light stimulus was presented against a background of dim light. In order to describe the results of the four resetting experiments, the following major refinements to the original model were necessary: (i) the relationship between light intensity (I) and drive onto the pacemaker was reduced from I1/3 to I0.23 for light levels between 150 and 10,000 lux; (ii) the van der Pol oscillator from the original model was replaced with a higher-order limit cycle oscillator so that amplitude recovery is slower near the singularity and faster near the limit cycle; (iii) a direct effect of light on circadian period (tau x) was incorporated into the model such that as I increases, tau x decreases, which is in accordance with "Aschoff's rule". This refined model generates the following testable predictions: it should be difficult to enhance normal circadian amplitude via bright light; near the critical point of a type 0 phase response curve (PRC) the slope should be steeper than it is in a type 1 PRC; and circadian period measured during forced desynchrony should be directly affected by ambient light intensity.

NASA Discipline Regulatory Physiology

Representation of the Auroral and Polar Ionosphere in the International Reference Ionosphere (IRI)

This issue of Advances in Space Research presents a selection of papers that document the progress in developing and improving the International Reference Ionosphere (IRI), a widely used standard for the parameters that describe the Earths ionosphere. The core set of papers was presented during the 2010 General Assembly of the Committee on Space Research in Bremen, Germany in a session that focused on the representation of the auroral and polar ionosphere in the IRI model. In addition, papers were solicited and submitted from the scientific community in a general call for appropriate papers.

IRI

An Empirical Temperature Variance Source Model in Heated Jets

An acoustic analogy approach is implemented that models the sources of jet noise in heated jets. The equivalent sources of turbulent mixing noise are recognized as the differences between the fluctuating and Favre-averaged Reynolds stresses and enthalpy fluxes. While in a conventional acoustic analogy only Reynolds stress components are scrutinized for their noise generation properties, it is now accepted that a comprehensive source model should include the additional entropy source term. Following Goldstein s generalized acoustic analogy, the set of Euler equations are divided into two sets of equations that govern a non-radiating base flow plus its residual components. When the base flow is considered as a locally parallel mean flow, the residual equations may be rearranged to form an inhomogeneous third-order wave equation. A general solution is written subsequently using a Green s function method while all non-linear terms are treated as the equivalent sources of aerodynamic sound and are modeled accordingly. In a previous study, a specialized Reynolds-averaged Navier-Stokes (RANS) solver was implemented to compute the variance of thermal fluctuations that determine the enthalpy flux source strength. The main objective here is to present an empirical model capable of providing a reasonable estimate of the stagnation temperature variance in a jet. Such a model is parameterized as a function of the mean stagnation temperature gradient in the jet, and is evaluated using commonly available RANS solvers. The ensuing thermal source distribution is compared with measurements as well as computational result from a dedicated RANS solver that employs an enthalpy variance and dissipation rate model. Turbulent mixing noise predictions are presented for a wide range of jet temperature ratios from 1.0 to 3.20.

Khavaran, Abbas

Parameter sensitivity in the dynamics of rotor-bearing systems

When designing a rotor system it is frequently desirable to have at hand a set of design sensitivity coefficients which quantitatively predict a change in specific system characteristics to changes in design parameters. This paper presents eigenvalue sensitivity coefficients for the damped natural frequencies of whirl of general linear rotor system modelled by finite element discretization. In addition, a simple and direct method for calculation of the damped critical speeds is presented, which utilizes the eigenvalue sensitivity with respect to the spin speed. It is shown that the combination of design parameter and spin speed whirl frequency sensitivity coefficients may be used to also evaluate the damped critical speed sensitivity coefficients.

Rajan, M.

Nonuniform sampling techniques for antenna applications

A two-dimensional sampling technique, which can employ irregularly spaced samples (amplitude and phase) in order to generate the complete far-field patterns is presented. The technique implements a matrix inversion algorithm, which depends only on the nonuniform sampled data point locations and with no dependence on the actual field values at these points. A powerful simulation algorithm is presented to allow a real-life simulation of many reflector/feed configurations and to determine the usefulness of the nonuniform sampling technique for the copolar and cross-polar patterns. Additionally, an overlapped window concept and a generalized error simulation model are discussed to identify the stability of the technique for recovering the field data among the nonuniform sampled data. Numerical results are tailored for the pattern reconstruction of a 20-m offset reflector antenna operating at L-band. This reflector is planned to be used in a proposed measurement concept of large antenna aboard the Space Shuttle, whereby it would be almost impractical to accurately control the movement of the Shuttle with respect to the RF source in prescribed directions in order to generate uniform sampled points. Also, application of the nonuniform sampling technique to patterns obtained using near-field measured data is demonstrated. Finally, results of an actual far-field measurement are presented for the construction of patterns of a reflector antenna from a set of nonuniformly distributed measured amplitude and phase data.

Rahmat-Samii, Yahya

Modeling the Diffuse Cloud-Top Optical Emissions from Ground and Cloud Flashes

A number of studies have indicated that the diffuse cloud-top optical emissions from intra-cloud (IC) lightning are brighter than that from normal negative cloud-to-ground (CG) lightning, and hence would be easier to detect from a space-based sensor. The primary reason provided to substantiate this claim has been that the IC is at a higher altitude within the cloud and therefore is less obscured by the cloud multiple scattering medium. CGs at lower altitudes embedded deep within the cloud are more obscured, so CG detection is thought to be more difficult. However, other authors claim that because the CG source current (and hence luminosity) is typically substantially larger than IC currents, the greater CG source luminosity is large enough to overcome the effects of multiple scattering. These investigators suggest that the diffuse cloud top emissions from CGs are brighter than from ICs, and hence are easier to detect from space. Still other investigators claim that the detection efficiency of CGs and ICs is about the same because modern detector sensitivity is good enough to "see" either flash type no matter which produces a brighter cloud top emission. To better assess which of these opinions should be accepted, we introduce an extension of a Boltzmann lightning radiative transfer model previously developed. It considers characteristics of the cloud (geometry, dimensions, scattering properties) and specific lightning channel properties (length, geometry, location, current, optical wave front propagation speed/direction). As such, it represents the most detailed modeling effort to date. At least in the few cases studied thus far, it was found that IC flashes appear brighter at cloud top than the lower altitude negative ground flashes, but additional model runs are to be examined before finalizing our general conclusions.

Solakiewicz, Richard

How Should Life Support Be Modeled and Simulated?

Why do most space life support research groups build and investigate large models for systems simulation? The need for them seems accepted, but are we asking the right questions and solving the real problems? The modeling results leave many questions unanswered. How then should space life support be modeled and simulated? Life support system research and development uses modeling and simulation to study dynamic behavior as part of systems engineering and analysis. It is used to size material flows and buffers and plan contingent operations. A DoD sponsored study used the systems engineering approach to define a set of best practices for modeling and simulation. These best practices describe a systems engineering process of developing and validating requirements, defining and analyzing the model concept, and designing and testing the model. Other general principles for modeling and simulation are presented. Some specific additional advice includes performing a static analysis before developing a dynamic simulation, applying the mass and energy conservation laws, modeling on the appropriate system level, using simplified subsystem representations, designing the model to solve a specific problem, and testing the model on several different problems. Modeling and simulation is necessary in life support design but many problems are outside its scope.

Jones, Harry W.

Multiphysics Modeling for Dimensional Analysis of a Self-Heated Molten Regolith Electrolysis Reactor for Oxygen and Metals Production on the Moon and Mars

The technology of direct electrolysis of molten lunar regolith to produce oxygen and molten metal alloys has progressed greatly in the last few years. The development of long-lasting inert anodes and cathode designs as well as techniques for the removal of molten products from the reactor has been demonstrated. The containment of chemically aggressive oxide and metal melts is very difficult at the operating temperatures ca 1600 C. Containing the molten oxides in a regolith shell can solve this technical issue and can be achieved by designing a self-heating reactor in which the electrolytic currents generate enough Joule heat to create a molten bath. In a first phase, a thermal analysis model was built to study the formation of a melt of lunar basaltic regolith irradiated by a focused solar beam This mode of heating was selected because it relies on radiative heat transfer, which is the dominant mode of transfer of energy in melts at 1600 C. Knowing and setting the Gaussian-type heat flux from the concentrated solar beam and the phase and temperature dependent thermal properties, the model predicts the dimensions and temperature profile of the melt. A validation of the model is presented in this paper through the experimental formation of a spherical cap melt realized by others. The Orbitec/PSI experimental setup uses an 3.6-cm diameter concentrated solar beam to create a hemispheric melt in a bed of lunar regolith simulant contained in a large pot. Upon cooling, the dimensions of the vitrified melt are measured to validate the thermal model. In a second phase, the model is augmented by multiphysics components to compute the passage of electrical currents between electrodes inserted in the molten regolith. The current through the melt generates Joule heating due to the high resistivity of the medium and this energy is transferred into the melt by conduction, convection and primarily by radiation. The model faces challenges in two major areas, the change of phase as temperature increases, and the dominance of radiative heat flux as heat transfer mechanism within the melt the change of phase concerns the regolith itself which is present in states ranging from a fine grain regolith with low thermal conductivity and low density to a vitrified melt with much higher thermal conductivity, and higher density. As the regolith is heated, it starts to soften around 1300 C the melt iS very viscous and evolving gas bubbles out in thick, lava-like fashion. By 1600 C the regolith is completely melted and the viscosity is low The second challenge resides in the proper modeling of the radiative heat flux requiring the addition of the computing-demanding radiative-heat-transfer function to the general heat transfer equation. The model Includes temperature-dependent properties (density, thermal conductivity, heat capacity, and viscosity, and absorption coefficients) and solves the radiative heat flux equation assuming gray (fine grains) and semi-transparent (melt) media and using an absorption coefficient spectral found in the literature for terrestrial minerals similar in composition to those of lunar regolith simulant

Dominguez, Jesus A.

Flux-Pinned Dynamics Model Parameterization and Sensitivity Study

Flux-pinned interfaces for spacecraft are an action-at-a-distance technology that can maintain a passively stable equilibrium between two spacecraft in close-proximity using the physics of magnetic flux pinning. Although flux pinning dynamics have been studied from a material-science perspective and at an interface level, there is a need to better understand the sensitivities and implications of system-level designs on the flux-pinned interface dynamics, especially in designs with multiple magnets and superconductors. These interfaces have highly nonlinear, coupled dynamics that are influenced by physical parameters including but not limited to strength of magnetic field sources, field-cooled position, and superconductor geometry. This paper addresses that gap by codifying parametric terms into an improved dynamics model, which can then be used to simulate the interaction of a multiple-superconductor-multiple-magnet interface. A standard starting point for modeling flux pinning dynamics is Kordyuk’s frozen image model, which defines a geometric mapping between magnetic field sources and their corresponding magnetic point source “images inside the volume of the superconductor.” The frozen image model successfully approximates the characteristics of flux pinning dynamics, but could provide more precise position and orientation predictions with the addition of various physical parameter refinements. The sensitivity of the general flux-pinned dynamics model is studied by varying the physical parameters and simulating the systems level dynamics. A predictive dynamics model is crucial to the maturation of this technology so it can be utilized in spacecraft systems, and this work represents a critical step in the development of that model.

Peck, Mason

Using Power Diagrams to Build Optimal Unstructured Meshes for C-Grid Models

The Model for Prediction Across Scales (MPAS) for Ocean (-O), Sea-Ice (-SI) and Land-Ice (-LI), in addition to the Coastal Ocean Marine Prediction Across Scales (COMPAS) are two novel general circulation models designed to resolve coupled ocean-ice dynamics over variable spatial scales using non-uniform unstructured grids. Both models are based on a conservative mimetic finite-difference/volume formulation (TRiSK), in which staggered momentum, vorticity and mass-based degrees- of-freedom are distributed over an orthogonal 'primal-dual' mesh.

Engwirda, Darren

OC7 phase I: Toward practical sea-state-dependent modeling of hydrodynamic viscous drag and damping

Here, this article presents a collaborative research campaign under the OC7 project on refining the engineering modeling approach for hydrodynamic viscous drag and damping on floating wind platforms, focusing on the adjustment of hydrodynamic drag and damping coefficients for different sea states. The participant simulation results show significant improvements over the previous OC6 project in predicting the low-frequency resonance motion under nonoperational conditions. The improvements are mainly due to enhanced modeling, including the adoption of wave stretching, and directly tuning the coefficients to measured platform motion in waves instead of free decay. For accurate predictions of mean- and slow-drift motion, the better performing models use a decreasing column splash zone drag coefficient and increasing surge damping/drag with increasing wave height. The model tuning for heave and pitch resonance shows less consistency. Generally, both quadratic drag and additional heave or pitch damping are needed for accurate predictions. Alternatively, a quadratic drag formulation with velocity filtering for the rectangular pontoons leads to improved predictions without additional damping. This model is also potentially more predictive, requiring minimal adjustment to its parameters for different conditions.

17 WIND ENERGY

A mathematical model of a single main rotor helicopter for piloted simulation

A mathematical model, suitable for piloted simulation of the flying qualities of helicopters, is a nonlinear, total force and moment model of a single main rotor helicopter. The model has ten degrees of freedom: six rigid body, three rotor flapping, and the rotor rotational degrees of freedom. The rotor model assumes rigid blades with rotor forces and moments radially integrated and summed about the azimuth. The fuselage aerodynamic model uses a detailed representation over a nominal angle of attack and sideslip range of + or - 15 deg., as well as a simplified curve fit at large angles of attack or sideslip. Stabilizing surface aerodynamics are modeled with a lift curve slope between stall limits and a general curve fit for large angles of attack. A generalized stability and control augmentation system is described. Additional computer subroutines provide options for a simplified engine/governor model, atmospheric turbulence, and a linearized six degree of freedom dynamic model for stability and control analysis.

Talbot, P. D.

On the stochastic dissemination of faults in an admissible network

The dynamic distribution of faults in a general type network is discussed. The starting point is a uniquely branched network in which each pair of nodes is connected by a single branch. Mathematical expressions for the uniquely branched network transition matrix are derived to show that sufficient stationarity exists to ensure the validity of the use of the Markov Chain model to analyze networks. In addition the conditions for the use of Semi-Markov models are discussed. General mathematical expressions are derived in an examination of branch redundancy techniques commonly used to increase reliability.

Kyrala, A.

Modeling the pressure-dilatation correlation

It is generally accepted that pressure dilatation, which is an additional compressibility term in turbulence transport equations, may be important for high speed flows. Recent direct simulations of homogeneous shear turbulence have given concrete evidence that the pressure dilatation is important insofar that it contributes to the reduced growth of turbulent kinetic energy due to compressibility effects. The problem of modeling pressure dilatation is addressed. A component of the pressure dilatation is isolated which exhibits temporal oscillations and, using direct numerical simulations of homogeneous shear turbulence and isotropic turbulence, show that it has a negligible contribution to the evolution of turbulent kinetic energy. Then, an analysis for the case of homogeneous turbulence is performed to obtain a model for the nonoscillatory pressure dilatation. This model algebraically relates the pressure dilatation to quantities traditionally obtained in incompressible turbulence closures. The model is validated by direct comparison with the pressure dilatation data obtained from the simulations.

Sarkar, S.