Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation optimization”

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 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Implementation of a Mesh refinement algorithm into the quasi-static PIC code QuickPIC

Plasma-based acceleration (PBA) has emerged as a promising candidate for the accelerator technology used to build a future linear collider and/or an advanced light source. In PBA, a trailing or witness particle beam is accelerated in the plasma wave wakefield (WF) created by a laser or particle beam driver. The WF is often nonlinear and involves the crossing of plasma particle trajectories in real space and thus particle-in-cell methods are used. The distance over which the drive beam evolves is several orders of magnitude larger than the wake wavelength. This large disparity in length scales is amenable to the quasi-static approach. Three-dimensional (3D), quasi-static (QS), particle-in-cell (PIC) codes, e.g., QuickPIC, have been shown to provide high fidelity simulation capability with 2-4 orders of magnitude speedup over 3D fully explicit PIC codes. In PBA, the witness beam needs to be matched to the focusing forces of the WF to reduce the emittance growth. In some linear collider designs, the matched spot size of the witness beam can be 2 to 3 orders of magnitude smaller than the spot size (and wavelength) of the wakefield. Such an additional disparity in length scales is ideal for mesh refinement where the WF within the witness beam is described on a finer mesh than the rest of the WF. A mesh refinement scheme is described that has been implemented into the 3D QS PIC code, QuickPIC. Very fine (high) resolution is used in a small spatial region that includes the witness beam and progressively coarser resolutions in the rest of the simulation domain. A fast multigrid Poisson solver has been implemented for the field solve on the refined meshes and a Fast Fourier Transform (FFT) based Poisson solver is used for the coarse mesh. The code has been parallelized with both MPI and OpenMP, and the parallel scalability has also been improved by using pipelining. A preliminary adaptive mesh refinement technique is described to optimize the computational time for simulations with an evolving witness beam size. Several test problems are used to verify that the mesh refinement algorithm provides accurate results. Additionally, the results are benchmarked against highly resolved simulations exhibiting near-azimuthal symmetry, performed using QPAD—a novel hybrid QS PIC code that uses a PIC description in the coordinates (r, ct – z) and a gridless description in the azimuthal angle, Φ.

Linear collider↗
DOE OSTI2017

Development of the IBSAL-SimMOpt Method for the Optimization of Quality in a Corn Stover Supply Chain

Variability on the physical characteristics of feedstock has a relevant effect on the reactor’s reliability and operating cost. Most of the models developed to optimize biomass supply chains have failed to quantify the effect of biomass quality and preprocessing operations required to meet biomass specifications on overall cost and performance. The Integrated Biomass Supply Analysis and Logistics (IBSAL) model estimates the harvesting, collection, transportation, and storage cost while considering the stochastic behavior of the field-to-biorefinery supply chain. This paper proposes an IBSAL-SimMOpt (Simulation-based Multi-Objective Optimization) method for optimizing the biomass quality and costs associated with the efforts needed to meet conversion technology specifications. The method is developed in two phases. For the first phase, a SimMOpt tool that interacts with the extended IBSAL is developed. For the second phase, the baseline IBSAL model is extended so that the cost for meeting and/or penalization for failing in meeting specifications are considered. The IBSAL-SimMOpt method is designed to optimize quality characteristics of biomass, cost related to activities intended to improve the quality of feedstock, and the penalization cost. A case study based on 1916 farms in Ontario, Canada is considered for testing the proposed method. Analysis of the results demonstrates that this method is able to find a high-quality set of non-dominated solutions.

09 BIOMASS FUELS↗
NASA NTRS1989

Computer simulation of a pilot in V/STOL aircraft control loops

The objective was to develop a computerized adaptive pilot model for the computer model of the research aircraft, the Harrier II AV-8B V/STOL with special emphasis on propulsion control. In fact, two versions of the adaptive pilot are given. The first, simply called the Adaptive Control Model (ACM) of a pilot includes a parameter estimation algorithm for the parameters of the aircraft and an adaption scheme based on the root locus of the poles of the pilot controlled aircraft. The second, called the Optimal Control Model of the pilot (OCM), includes an adaption algorithm and an optimal control algorithm. These computer simulations were developed as a part of the ongoing research program in pilot model simulation supported by NASA Lewis from April 1, 1985 to August 30, 1986 under NASA Grant NAG 3-606 and from September 1, 1986 through November 30, 1988 under NASA Grant NAG 3-729. Once installed, these pilot models permitted the computer simulation of the pilot model to close all of the control loops normally closed by a pilot actually manipulating the control variables. The current version of this has permitted a baseline comparison of various qualitative and quantitative performance indices for propulsion control, the control loops and the work load on the pilot. Actual data for an aircraft flown by a human pilot furnished by NASA was compared to the outputs furnished by the computerized pilot and found to be favorable.

Vogt, William G.↗
NASA NTRS2005

Predicting Flows of Rarefied Gases

DSMC Analysis Code (DAC) is a flexible, highly automated, easy-to-use computer program for predicting flows of rarefied gases -- especially flows of upper-atmospheric, propulsion, and vented gases impinging on spacecraft surfaces. DAC implements the direct simulation Monte Carlo (DSMC) method, which is widely recognized as standard for simulating flows at densities so low that the continuum-based equations of computational fluid dynamics are invalid. DAC enables users to model complex surface shapes and boundary conditions quickly and easily. The discretization of a flow field into computational grids is automated, thereby relieving the user of a traditionally time-consuming task while ensuring (1) appropriate refinement of grids throughout the computational domain, (2) determination of optimal settings for temporal discretization and other simulation parameters, and (3) satisfaction of the fundamental constraints of the method. In so doing, DAC ensures an accurate and efficient simulation. In addition, DAC can utilize parallel processing to reduce computation time. The domain decomposition needed for parallel processing is completely automated, and the software employs a dynamic load-balancing mechanism to ensure optimal parallel efficiency throughout the simulation.

LeBeau, Gerald J.↗
DOE OSTI2019

Simulation and design of a compact GaAs based tunable dual-wavelength diode laser system

We present our design of a compact, integrated and tunable dual-wavelength diode laser system emitting around 785 nm, which is of interest for several applications like Raman spectroscopy and the generation of THz radiation. To achieve a more compact device compared to previous GaAs based designs two etch depths are realized, leading to shallowly etched ridge waveguides in regions were optical gain is applied and deeply etched waveguides used to enable compact integrated waveguide components. The device parameters are optimized using a numerically efficient simulation tool for passive waveguides. Subsequently, the entire laser system is further analyzed applying a sophisticated traveling-wave equation based model for active devices giving access to internal intensity and carrier density distributions. It is shown that active laser simulations are crucial to deduce critical and performance limiting design aspects not accessible via an all-passive simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Simulation study for particle identification with the dRICH of the ePIC experiment at the EIC

The dual-radiator Imaging Cherenkov detector (dRICH), employing an aerogel and a gas radiator, is a key component of the forward particle identification system for the ePIC experiment at the Electron-Ion Collider (EIC). This study evaluates the dRICH performance using Geant4 simulations in the context of the global ePIC simulation stack, focusing on the optimization of the aerogel radiator and the impact of sensor noise. We compare two aerogel configurations: the initial design (n = 1.019) and the current default (n = 1.026). The latter, characterized by improved optical properties and a higher refractive index, demonstrates enhanced π/K separation at high momenta, effectively extending the operational overlap with the C 2 F 6 gas radiator (n = 1.00076 at 25 °C, as implemented in the simulation software). Additionally, the study investigates the impact of Silicon Photomultiplier (SiPM) dark noise, showing that a 300 kHz noise rate per 3mm x 3mm channel leads to a moderate reduction (approximately 1.5 GeV/c) in the 3σ separation threshold. These results validate the current dRICH design and quantify the purity levels achievable for both radiators under expected experimental conditions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

An iterative forward analysis technique to determine the equation of state of dynamically compressed materials

Here, we developed an iterative forward analysis (IFA) technique with the ability to use hydrocode simulations as a fitting function for analysis of dynamic compression experiments. The IFA method optimizes over parameterized quantities in the hydrocode simulations, breaking the degeneracy of contributions to the measured material response. Velocity profiles from synthetic data generated using a hydrocode simulation are analyzed as a first-order validation of the technique. We also analyze multiple magnetically driven ramp compression experiments on copper and compare with more conventional techniques. Excellent agreement is obtained in both cases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Mixing with applications to inertial-confinement-fusion implosions

Approximate one-dimensional (1D) as well as 2D and 3D simulations are playing an important supporting role in the design and analysis of future experiments at National Ignition Facility. This work is mainly concerned with 1D simulations, used extensively in design and optimization. We couple a 1D buoyancy-drag mix model for the mixing zone edges with a 1D inertial confinement fusion simulation code. This analysis predicts that National Ignition Campaign (NIC) designs are located close to a performance cliff, so modeling errors, design features (fill tube and tent) and additional, unmodeled instabilities could lead to significant levels of mix. The performance cliff we identify is associated with multimode plastic ablator (CH) mix into the hot-spot deuterium and tritium (DT). The buoyancy-drag mix model is mode number independent and selects implicitly a range of maximum growth modes. Our main conclusion is that single effect instabilities are predicted not to lead to hot-spot mix, while combined mode mixing effects are predicted to affect hot-spot thermodynamics and possibly hot-spot mix. Combined with the stagnation Rayleigh-Taylor instability, we find the potential for mix effects in combination with the ice-to-gas DT boundary, numerical effects of Eulerian species CH concentration diffusion, and ablation-driven instabilities. With the help of a convenient package of plasma transport parameters developed here, we give an approximate determination of these quantities in the regime relevant to the NIC experiments, while ruling out a variety of mix possibilities. Plasma transport parameters affect the 1D buoyancy-drag mix model primarily through its phenomenological drag coefficient as well as the 1D hydro model to which the buoyancy-drag equation is coupled.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1988

Cooperative control of two arms in the transport of an inertial load in zero gravity

In designing a robot control system for dual arm configurations, the control engineer is faced with two challenges: to derive the equations of motion for a given situation, and to meet certain desired control requirements (for instance, minimum energy). The former may involve closed kinematic chains, such as the case when the two arms are grasping a common object. The latter usually involves nonlinear optimization. These issues are considered in the context of transporting an inertial load using two planar three-link arms. A generalized 'reduction transformation' is applied to the dynamics to remove the singularity in the system equations. A suboptimal minimum energy method is presented to reduce a difficult 12-state, six-control nonlinear optimization to two independent, nonconflicting suboptimizations. A simulation example is provided to illustrate the degree of energy reduction possible using the optimal arm torque distribution that was developed.

Carignan, Craig R.↗
NASA NTRS1996

Advanced Models for Aeroelastic Analysis of Propulsion Systems

This report describes an integrated, multidisciplinary simulation capability for aeroelastic analysis and optimization of advanced propulsion systems. This research is intended to improve engine development, acquisition, and maintenance costs. One of the proposed simulations is aeroelasticity of blades, cowls, and struts in an ultra-high bypass fan. These ducted fans are expected to have significant performance, fuel, and noise improvements over existing engines. An interface program was written to use modal information from COBSTAN and NASTRAN blade models in aeroelastic analysis with a single rotation ducted fan aerodynamic code.

Keith, Theo G., Jr.↗
NASA NTRSDate not supplied

Estimation of Timing Performance with Pulsed Laser Links for Small Satellite PNT

Satellite constellations are growing in size and in utility for applications as diverse as providing reliable, low-latency internet service to rural areas orEarth science missions.This expansion of satellite constellations brings attention to the need for responsive and reliable satellite communication.Current research in the PrecisionSpace Systems Laboratory at the University of Florida focuses on using pulsed laser communications as the method of relaying time signals instead of the radio frequency (RF)signals.The upcoming CLICK (CubeSat Laser Infrared CrosslinK) missions (Massachusetts Institute of Technology, University of Florida, NASA Ames Research Center) will demonstrate a spacecraft-to-spacecraft laser communication link and time-transfer.For pre-flight technology assessment of missions such as CLICK, a simulation of a spacecraft constellation in Earth orbit has been created to estimate clock synchronization and precision orbit determination based on measured instrumentation performance.We have developed a novel MATLAB-based numerical simulator to model spacecraft-to-spacecraft laser time-transfer and estimate the offset between the spacecraft clocks over time.This simulation includes timing errors associated with laser pulse generation and detection, as well as non-Gaussian clock drift models.The two on-board clocks modeled are a cesium-based Chip-Scale Atomic Clock (CSAC) and a rubidium-based Miniature Atomic Clock (MAC), both produced by Microchip.An example case of two spacecraft in a circular, low-Earth orbit receiving GPS position, GPS timing, and laser pulse time-of-flight measurements is simulated.The positions and velocities of the two spacecraft at a reference epoch and the constant clock model coefficients are estimated.Polynomial models of different orders were used as clock models. The effect of clock model order on the root-mean-square (RMS) of the clock error is apparent in the case using GPS and lasing measurements, showing that the clock model improves with increasing clock model order. Results compare the estimated clock model of a mission operation that only uses GPS measurements and a mission operation that uses both GPS and laser pulse time-of-flight measurements between spacecraft referenced to their on-board CSACs or on-board MACs. Including lasing measurements reduces the RMS clock model error by approximately 80% of the RMS of the cases with only GPS measurements.This simulation tool can be used to optimize the lasing operations schedule based on mission timing performance objectives. During this presentation, I will discuss the details of the simulation and clock model estimation, as well as discuss the results.I will compare the results of incorporating lasing timing measurements to using only GPS measurements and the outcomes of using different orders of polynomials.

Timing Performance↗
DOE OSTI2022

Digital Twin Development for Real Time Optimization

In the real-time optimization operation process of the integrated energy system, digital twin (DT) can support operators to find optimal dispatching strategy. In this research, the DT framework interfacing (external) system-level simulation model, RAVEN, and Python-based optimization framework has been developed. Simple demonstration of finding optimal input values minimizing an objective function with the DT framework has been presented.

25 ENERGY STORAGE↗
DOE OSTI2026

Adaptive Computing and Multi-Fidelity Learning

We describe our ongoing research in adaptive computing. Our goal is to use a combination of low- and high-fidelity simulation models to enable computationally efficient optimization and uncertainty quantification. We develop optimization formulations that take into account the compute resources currently available, which act as a constraint with regards to the fidelity level simulation we can run while maximizing information gain. We will discuss a few application examples that can benefit from this approach, especially when considering challenges arising in scaling up experiments and simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A Hands-on Course on Intensified Membrane Processes for Sustainable Water Purification

The lack in process intensification (PI) training may be limiting processes to achieve sustainability goals. Furthermore, this project developed a hands-on course on intensified membrane processes to address knowledge gaps in PI. Bench-scale activities introduced concepts of membrane separation, operating engineering-scale modules provided a near-industrial experience, and simulation software allowed modeling and optimization of a process. Student evaluations showed that the combination of modeling, experimental, and peer-instruction work helped them understand the application of theoretical concepts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Initial conceptual demonstration of control co-design for WEC optimization

Abstract While some engineering fields have benefited from systematic design optimization studies, wave energy converters have yet to successfully incorporate such analyses into practical engineering workflows. The current iterative approach to wave energy converter design leads to sub-optimal solutions. This short paper presents an open-source MATLAB toolbox for performing design optimization studies on wave energy converters where power take-off behavior and realistic constraints can be easily included. This tool incorporates an adaptable control co-design approach, in that a constrained optimal controller is used to simulate device dynamics and populate an arbitrary objective function of the user’s choosing. A brief explanation of the tool’s structure and underlying theory is presented. To demonstrate the capabilities of the tool, verify its functionality, and begin to explore some basic wave energy converter design relationships, three conceptual case studies are presented. In particular, the importance of considering (and constraining) the magnitudes of device motion and forces in design optimization is shown.

16 TIDAL AND WAVE POWER↗
DOE OSTI2019

Reduced energetic particle transport models enable comprehensive time-dependent tokamak simulations

Time-dependent integrated simulations through codes such as TRANSP are becoming an indispensable tool for the interpretation of existing experiments and predictions of optimized scenarios. For many practical cases, quantitative simulations need to include the effect of plasma instabilities on the evolution of a tokamak discharge. An example is the degradation in energetic particle (EP) confinement induced by instabilities, which in turn affects important source terms for heating, non-inductive current, and momentum in a simulation. The reduced-physics "kick model" provides phase-space resolved transport probability matrices to TRANSP that are used to account for enhanced EP transport by instabilities in addition to neoclassical transport. The model has recovered the measured Alfvén eigenmode (AE) spectrum on NSTX, NSTX-U and DIII-D, and has reproduced details of phase-space resolved fast ion diagnostic data measured on DIII-D for EP-driven modes and tearing modes. In general, the kick model has proven the potential of phase-space resolved EP simulations to unravel details of EP transport for detailed theory/experiment comparison and for scenario planning based on optimization of Neutral Beam (NB) injection parameters. In this work, the extension of the kick model to low-frequency instabilities such as Tearing Modes and fishbones, in addition to AEs, is assessed. The goal is to enable TRANSP simulations that retain the main effects of multiple types of instabilities through a common framework. Results from the NSTX/NSTX-U and DIII-D tokamaks show that the extension to "multi-mode" scenarios can expand the range of applicability of the model for more reliable, quantitative integrated simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

AC Resistance Reduction Using Orthogonal Air Gaps in High-Frequency Inductors: Preprint

This paper presents a relatively simple technique to reduce winding losses due to fringing fields in high-frequency inductors. In high-frequency power electronics, ac inductor winding losses are affected by skin and proximity effects, including uneven current distribution due to fringing magnetic fields around airgaps. It is well known how fringing effects can be mitigated using distributed airgaps, at the expense of non-standard core or winding geometry. In planar magnetics, the orthogonal-airgap approach proposed in this paper combines airgaps in core segments parallel with the windings with conventionally placed airgaps in segments perpendicular to the windings. The approach is developed using a 1D analytical framework and validated by 2D finite-element simulations. Analytical guidelines are presented to optimize the airgaps to achieve minimum ac resistance. As a case study, a planar inductor is designed for an 8~kW SiC-based buck converter operating at 250~kHz. It is shown how the orthogonal airgaps result in more than 45\% reduction in ac resistance and substantially reduced inductor losses compared to the design using standard airgaps. The results are verified by loss measurements on an experimental converter prototype.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

Numerical investigation on a new designed hollow-core photonic crystal fiber with large modal separation

We investigate the characteristics of optical modes in a new designed hollow-core photonic crystal fiber (PCF) with large modal separation. The large effective refractive index difference between the adjacent second-order vector modes can be maintained higher than 0.001 within the wavelength range from 1 to 2 μm. The maximum value can arrive 0.01199 between HE{sub 21} and TM{sub 01} mode at 2 μm. The fundamental mode of the proposed PCF has flat dispersion in a wide spectral range, which can be applied on wideband supercontinuum generation. The influence on the characteristics by the PCF structure has also been investigated. The increase of the size of the central air hole and the duty cycle affects the flatness of the dispersion curve to some extent. Our simulated results could provide an optimized fiber structure for the experimental study on stable transmission of vector modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗