Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coupled process modeling”

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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

A physics-based crystallographic modeling framework for describing the thermal creep behavior of Fe-Cr alloys

In this work, a physics-based thermal creep model is developed based on the understanding of the microstructure in Fe-Cr alloys. This model is associated with a transition state theory based framework that considers the distribution of internal stresses at sub-material point level. The thermally activated dislocation glide and climb mechanisms are coupled in the obstacle-bypass processes for both dislocation and precipitate-type barriers. A kinetic law is proposed to track the dislocation densities evolution in the subgrain interior and in the cell wall. The predicted results show that this model, embedded in the visco-plastic self-consistent (VPSC) framework, captures well the creep behaviors for primary and steady-state stages under various loading conditions. We also discuss the roles of the mechanisms involved.

36 MATERIALS SCIENCE↗
NASA NTRS2001

Volcanoes, Polar Clouds and Arctic Ozone

Satellite observations and model calculations show 5 to 10% local column ozone loss in some tropical and mid latitude locations, following El Chichon and Mount Pinatubo eruptions. The rapid deepening of the Antarctic ozone hole in the early 1980s has also been partially attributed to chemistry on volcanic aerosols from a number of large eruptions. Here the effects of volcanoes on Arctic polar processes are explored. Large polar stratospheric cloud particles that cause denitrification cannot form in a volcanically perturbed environment. Denitrification can increase Arctic ozone loss by up to 30% in a future colder climate. However, we show that enhanced chemical processing on volcanic aerosols can increase Arctic ozone loss in a cold year by about 60% independent of denitrification. A coupled chemistry-microphysics model is used to show that widespread distribution of volcanic aerosols in 2000 could have caused severe springtime ozone depletion in the Arctic stratosphere. While, volcanic aerosols can strongly affect the current Arctic column ozone abundance in a cold year, denitrification effects on ozone can only become important in a much colder lower stratosphere.

Tabazadeh, Azadeh↗
DOE OSTI2014

Semi-analytic models for the CANDELS survey: comparison of predictions for intrinsic galaxy properties

We compare the predictions of three independently developed semi-analytic galaxy formation models (SAMs) that are being used to aid in the interpretation of results from the CANDELS survey. These models are each applied to the same set of halo merger trees extracted from the 'Bolshoi' high-resolution cosmological N-body simulation and are carefully tuned to match the local galaxy stellar mass function using the powerful method of Bayesian Inference coupled with Markov Chain Monte Carlo or by hand. The comparisons reveal that in spite of the significantly different parameterizations for star formation and feedback processes, the three models yield qualitatively similar predictions for the assembly histories of galaxy stellar mass and star formation over cosmic time. Comparing SAM predictions with existing estimates of the stellar mass function from z = 0-8, we show that the SAMs generally require strong outflows to suppress star formation in low-mass halos to match the present-day stellar mass function, as is the present common wisdom. However, all of the models considered produce predictions for the star formation rates (SFRs) and metallicities of low-mass galaxies that are inconsistent with existing data. The predictions for metallicity-stellar mass relations and their evolution clearly diverge between the models. We suggest that large differences in the metallicity relations and small differences in the stellar mass assembly histories of model galaxies stem from different assumptions for the outflow mass-loading factor produced by feedback. Importantly, while more accurate observational measurements for stellar mass, SFR and metallicity of galaxies at 1 < z < 5 will discriminate between models, the discrepancies between the constrained models and existing data of these observables have already revealed challenging problems in understanding star formation and its feedback in galaxy formation. The three sets of models are being used to construct catalogs of mock galaxies on light cones that have the same geometry as the CANDELS survey, which should be particularly useful for quantifying the biases and uncertainties on measurements and inferences from the real observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Road Map for Development of Crystal-Tolerant High Level Waste Glasses

The U.S. Department of Energy (DOE) is building a Tank Waste Treatment and Immobilization Plant (WTP) at the Hanford Site in Washington to remediate 55 million gallons of radioactive waste that is being temporarily stored in 177 underground tanks. Efforts are being made to increase the loading of Hanford tank wastes in glass while meeting melter lifetime expectancies and process, regulatory, and product quality requirements. This road map guides the research and development for formulation and processing of crystal-tolerant glasses, identifying near- and long-term activities that need to be completed over the period from 2014 to 2019. The primary objective is to maximize waste loading for Hanford waste glasses without jeopardizing melter operation by crystal accumulation in the melter or melter discharge riser. The potential applicability to the Savannah River Site (SRS) Defense Waste Processing Facility (DWPF) will also be addressed in this road map. The planned research described in this road map is motivated by the potential for substantial economic benefits (significant reductions in glass volumes) that will be realized if the current constraints (T 1% for WTP and T L for DWPF) are approached in an appropriate and technically defensible manner for defense waste and current melter designs. The basis of this alternative approach is an empirical model predicting the crystal accumulation in the WTP glass discharge riser and melter bottom as a function of glass composition, time, and temperature. When coupled with an associated operating limit (e.g., the maximum tolerable thickness of an accumulated layer of crystals), this model could then be integrated into the process control algorithms to formulate crystal-tolerant high-level waste (HLW) glasses targeting high waste loadings while still meeting process related limits and melter lifetime expectancies. The modeling effort will be an iterative process, where model form and a broader range of conditions, e.g., glass composition and temperature, will evolve as additional data on crystal accumulation are gathered. Model validation steps will be included to guide the development process and ensure the value of the effort (i.e., increased waste loading and waste throughput). A summary of the stages of the road map for developing the crystal-tolerant glass approach, their estimated durations, and deliverables is provided.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2011

Search for a Light Higgs Boson at BaBar

We search for evidence of a light Higgs boson (A{sup 0}) in the radiative decays of the narrow {Upsilon}(3S) resonance: {Upsilon}(3S) {yields} {gamma}A{sup 0}, where A{sup 0} {yields} invisible or A{sup 0} {yields} {mu}{sup +}{mu}{sup -}. Such an object appears in extensions of the Standard Model, where a light CP-odd Higgs boson naturally couples strongly to b-quarks. We find no evidence for such processes in a sample of 122 x 106 {Upsilon}(3S) decays collected by the BABAR collaboration at the PEP-II B-factory, and set 90% C.L. upper limits on the product of the branching fractions {Beta}({Upsilon}(3S) {yields} {gamma}A{sup 0}) x {Beta}(A{sup 0} {yields} invisible) at (0.7-31) x 10{sup -6} in the mass range mA{sup 0} {le} 7.8 GeV, and on the product {Beta}({Upsilon}(3S) {yields} {gamma}A{sup 0}) x {Beta}(A{sup 0} {yields} {mu}{sup +}{mu}{sup -}) at (0.25-5.2) x 10{sup -6} in the mass range 0.212 {le} mA{sup 0} {le} 9.3GeV. We also set a limit on the dimuon branching fraction of the recently discovered {eta}{sub b} meson {Beta}({eta}{sub b} {yields} {mu}{sup +}{mu}{sup -}) < 0.8% at 90% C.L. The results are preliminary.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

A rigorous process modeling methodology for biomass fast pyrolysis with an entrained-flow reactor

A biomass fast pyrolysis model was developed for implementation in equation-oriented modeling software. Based on a previous framework of coupled 1-dimensional mass, momentum, and heat balance equations, this model includes updated reaction kinetics to provide a more detailed representation of biomass components, intermediates, and products. A recently published derivation of thermodynamic properties for the species in this model has allowed the energy balance around the pyrolysis reactor to be rigorously redefined. With these improvements, the optimum pyrolysis temperature for bio-oil production predicted by the model is increased by up to 50°C, bringing it in line with experimental data and increasing the overall agreement. More importantly, the reactor energy balance is strictly enforced. The resulting model can be used for the design and optimization of biomass fast pyrolysis processes, and comparisons with other options for biomass utilization.

09 BIOMASS FUELS↗
DOE OSTI2021

Development of Time Lapse VSP Integration Workflow: A Case Study at Farnsworth CO2-EOR Project

Abstract This study aims to develop a 4D Vertical Seismic Profile (VSP) integration workflow to improve the prediction of subsurface stress changes. The selected study site is a 5-spot pattern within the ongoing CO2-EOR operations at the Farnsworth Field Unit FWU in Ochiltree County, Texas. The specific pattern has undergone extensive geological and geomechanical characterization through the acquisition of 3D seismic data, geophysical well logs, and core. This workflow constrains a numerical hydromechanical model by applying a penalty function formed between "modeled" versus "observed" time-lapse compressional and shear seismic velocity changes. Analyses of geophysical logs and ultra-sonic measurements on core exhibit measurable sensitivities to changes in both fluid saturation and mean effective stress. These data are used to develop a site-specific rock physics model and stress-velocity relationship, which inform the numerical models used to generate the "modeled" portion of the penalty function. The "observed" portion of the penalty function is provided by a novel elastic full-waveform inversion of the available 3D baseline and three monitor surveys to produce high-quality estimates of time-lapse compressional and shear seismic velocity changes. The modeling workflow accounts sequentially for fluid substitution and stress impacts. Hydrodynamic and geomechanical properties of the 3D coupled numerical model are estimated through geostatistical integration of well log and core data with 3D seismic inversion products. Changes in seismic velocities due to fluid substitution are computed using the Biot-Gassmann workflow and site-specific rock physics. Stress impacts on time-lapse seismic velocity changes are modeled from the effective stress output of the hydromechanical model and are initially based on the velocity versus effective stress relationship extracted from core mechanical testing. Based on the principle of superposition of seismic wavefields, seismic velocity changes attributed to fluid substitution and that due to changes in mean effective stress are treated as linearly additive. The modeled results are upscaled using Backus averaging to reconcile scale discrepancies between the modeled and measured datasets to formulate the penalty function. This manuscript presents the forward modeling process and concludes that for the base case, the seismic velocity changes due to mean effective stress dominates over the seismic velocity changes attributed to fluid substitution because of the extensive range of the pressure perturbations. Successful minimization of this penalty function calibrates the coupled hydrodynamic geomechanical numerical model and affirms the suitability of acoustic time-lapse measurements such as 4D-VSP for geomechanical calibration.

02 PETROLEUM↗
DOE OSTI2012

The role of asymmetric excitation in self-organized nanostructure formation upon femtosecond laser ablation

Surface pattering upon multi-pulse femtosecond laser ablation is modeled by a non-linear-dynamic erosion/smoothing model, similar to structure formation during ion sputtering. The model is adopted to account for the influence of laser polarization on nanostructure features. Based on a nonlinear equation of the Kuramoto-Siavshinsky type, it is shown that the directional anisotropy in the pattern formation may result from a spatial anisotropy of the initial excitation/energy-coupling process, such as resonant coupling to surface plasmons/polaritons, or electron diffusion properties. Also, anisotropy of elasto-dynamic surface and diffusion properties may be involved. A comparison of numeric simulations based on the model with corresponding experi-mental results gives a very good agreement.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A self-consistent and spatially dependent model of the multiband emission of pulsar wind nebulae

A self-consistent and spatially dependent model is presented to investigate the multiband emission of pulsar wind nebulae (PWNe). In this model, a spherically symmetric system is assumed and the dynamical evolution of the PWN is included. The processes of convection, diffusion, adiabatic loss, radiative loss, and photon–photon pair production are taken into account in the electron’s evolution equation, and the processes of synchrotron radiation, inverse Compton scattering, synchrotron self-absorption, and pair production are included for the photon’s evolution equation. Both coupled equations are simultaneously solved. The model is applied to explain observed results of the PWN in MSH 15–52. Our results show that the spectral energy distributions (SEDs) of both electrons and photons are all a function of distance. The observed photon SED of MSH 15–52 can be well reproduced in this model. With the parameters obtained by fitting the observed SED, the spatial variations of photon index and surface brightness observed in the X-ray band can also be well reproduced. Moreover, it can be derived that the present-day diffusion coefficient of MSH 15–52 at the termination shock is κ{sub 0}=6.6×10{sup 24} cm{sup 2} s{sup −1}, the spatial average has a value of κ-bar =1.4×10{sup 25} cm{sup 2} s{sup −1}, and the present-day magnetic field at the termination shock has a value of B{sub 0}=26.6 μG and the spatial averaged magnetic field is B-bar =14.9 μG. The spatial changes of the spectral index and surface brightness at different bands are predicted.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Multi-fidelity modeling to predict the rheological properties of a suspension of fibers using neural networks and Gaussian processes

Unveiling the rheological properties of fiber suspensions is of paramount interest to many industrial applications. There are multiple factors, such as fiber aspect ratio and volume fraction, that play a significant role in altering the rheological behavior of suspensions. Three-dimensional (3D) numerical simulations of coupled differential equations of the suspension of fibers are computationally expensive and time-consuming. Machine learning algorithms can be trained on the available data and make predictions for the cases where no numerical data are available. However, some widely used machine learning surrogates, such as neural networks, require a relatively large training dataset to produce accurate predictions. Multi-fidelity models, which combine high-fidelity data from numerical simulations and less expensive lower fidelity data from resources such as simplified constitutive equations, can pave the way for more accurate predictions. Here, we focus on neural networks and the Gaussian processes with two levels of fidelity, i.e., high and low fidelity networks, to predict the steady-state rheological properties, and compare them to the single-fidelity network. High-fidelity data are obtained from direct numerical simulations based on an immersed boundary method to couple the fluid and solid motion. The low-fidelity data are produced by using constitutive equations. Multiple neural networks and the Gaussian process structures are used for the hyperparameter tuning purpose. Results indicate that with the best choice of hyperparameters, both the multi-fidelity Gaussian processes and neural networks are capable of making predictions with a high level of accuracy with neural networks demonstrating marginally better performance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A Kinetic Model for Predicting Trace Gas Uptake and Reaction

A model is developed to describe trace gas uptake and reaction with applications to aerosols and microdroplets. Gas uptake by the liquid is formulated as a coupled equilibria that links gas, surface, and bulk regions of the droplet or solution. Previously, this framework was used in explicit stochastic reaction-diffusion simulations to predict the reactive uptake kinetics of ozone with droplets containing aqueous aconitic acid, maleic acid, and sodium nitrite. With the use of prior data and simulation results, a new equation for the uptake coefficient is derived, which accounts for both surface and bulk reactions. Lambert W functions are used to obtain closed form solutions to the integrated rate laws for the multiphase kinetics; similar to previous expressions that describe Michaelis-Menten enzyme kinetics. Together these equations couple interface and bulk processes over a wide range of conditions and do not require many of the limiting assumptions needed to apply resistor model formulations to explain trace gas uptake and reaction.

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

Maximizing long-term biohydrogen production with Clostridium thermocellum for high solids conversion of lignocellulosic biomass

Biological hydrogen production from lignocellulosic biomass sustainably couples organic waste reduction with renewable energy generation. Efficient conversion is challenged by the structural complexity of lignocellulose and resulting recalcitrance to enzymatic degradation. Clostridium thermocellum natively breaks down biomass with highly effective hemi-/cellulases systems (i.e., cellulosomes) and generates hydrogen in anaerobic cultivation, creating a compelling platform for lignocellulosic biohydrogen production. Achieving commercially viable production rates requires balancing high biomass loading and throughput against uniform mixing conditions required for enzyme dispersion, pH and temperature control, and efficient hydrogen and metabolite removal in continuous operation. To address these barriers to process intensification, we implemented novel reactor and process designs for high-solids lignocellulosic biomass fermentations using the C. thermocellum KJC19-9 strain, genetically engineered for co-utilization of cellulose and hemicellulose sugars (i.e., xylose). Via computational fluid dynamics (CFD) modeling and experimental validation, we achieved a >50% improvement in biohydrogen production with an improved anchor-type impeller morphology, coupled to a threefold reduction in agitation rate. To further reduce rheological constraints and accumulation of toxic metabolites, we then transitioned the process to sequencing fed-batch operation. The resulting process generated 24.87 L H 2 L −1 from 160 g L −1 of deacetylated and mechanically refined (DMR)-pretreated corn stover biomass over 16 days while solubilizing >95% of influent cellulose and hemicellulose, setting a new performance benchmark for continuous production of biohydrogen from lignocellulose.

08 HYDROGEN↗
NASA NTRS1992

A fixed tuned broadband matching structure for submillimeter SIS receivers

The authors have designed, fabricated, and tested a quasi-optical submillimeter wave receiver with an Nb/AlO(x)/Nb tunnel junction. This design incorporates a hybrid antenna fed by a planar logarithmic spiral structure in order to couple to the radiation field from the telescope. The novel matching circuit requires several layers of photolithographic processing on top of the actual tunneling device. Computer modeling of the device correctly predicted the measured bandwidth and the characteristic frequencies to within 8 percent. A good match has been obtained from 200 to 475 GHz between the antenna and a relatively large area (1 sq mm) tunnel junction with omega(R(n))C about 2-4. Noise measurements at 318 GHz, 395 GHz, 426 GHz, and 492 GHz yielded uncorrected double sideband receiver noise temperatures of 200 K, 230 K, 220 K, and 500 K, respectively. Using the same optics, coupling efficiencies between the receiver and the Caltech Submillimeter Telescope were found to have values approaching those achieved by the best waveguide-horn-based receiver systems.

Buttgenbach, Thomas H.↗
DOE OSTI2016

Sensitivity Analysis of Passive Heat Removal Capability of FHR Based on Efficient Surrogate Model and GenFlow

Code Scaling, Applicability and Uncertainty Evaluation Methodology (CSAU) is used for quantifying reactor safety margins and works as an alternative method for licensing nuclear power plant. In the framework of CSAU, sensitivity and uncertainty analysis (SA and UA) should be performed. However, for large-scale complicated thermal hydraulics analysis simulated by code (e.g. RELAP5, CFD and subchannel analysis, the sensitivity and uncertainty analysis is really time-consuming. Furthermore, sometimes it is not feasible to perform SA and UA. Specifically, the computation time for Monte Carlo sampling or advanced LHS sampling can be neglected, but the computation time for thermal hydraulics code makes up main part of time during SA and UA process. In order to accelerate the speed of SA and UA, thermal hydraulics code can be improved based on advanced mathematic algorithm or using parallel method. However, the development for new mathematic method needs long time and strict verification and validation process is needed for a new code, which usually takes a couple of years or even longer time to have an approved code that can be used in nuclear industry. An alternative way to improve the speed of SA and UA is using surrogate model. The surrogate model is derived from the input as well as output data of the original code, assuming the code as a black box. The most common method is using response surface method (RSM) as a surrogate for the code, the input parameters are then sampled from their assumed distribution and the corresponding value of the response surface function is calculated. However, a couple of input parameters need hundreds or even thousands of original code calculations. The large number of runs is an obvious disadvantage of this method. Besides, the second disadvantage is an appropriate response surface function must be selected (e.g. quadratic, cubic, quartic, with or without cross product terms, or using a non-polynomial form) to avoid the overfit and underfit problem. This process is time consuming and expensive. Polynomial chaos expansion method is a reasonably fast and attractive alternative, stochastic model. The chaos expansion provides an efficient high-order method of including non-linear effects in stochastic analysis. Usually, the information about the model parameters and their probability distribution is limited and the dependence between them is not known. However, if arbitrary polynomial chaos is used, instead of assuming probability distribution for input parameters which may introduces additional subjectivity into the analysis, the available information can directly and most purely be used in stochastic analysis based on arbitrary data distribution. Thus, it is not necessary to know the exact probability density functions to perform SA and UA based on arbitrary polynomial chaos (aPC). The aPC has been applied for uncertainty quantification in many fields (e.g. in solid mechanics, fluid flows, thermal sciences), and become widely used in nuclear industry recently. Mk1 PB-FHR is a small modular fluoride-salt-cooled, high-temperature reactors (FHRs) that is a concept design by University of California, Berkeley (UCB). It has inherently safe features like natural-circulation-driven systems (e.g. direct reactor auxiliary cooling system (DRACS)) to passively remove decay heat from the reactor core during design basis and beyond design. It is necessary to perform SA and UA, since the licensing process of this new kind of reactor is under way and there are many design parameters that have not been fixed. Current nuclear commercial codes, such as RELAP5, have mainly been developed for light water reactors, which require the use of comprehensive two-phase flow models. This capability is not required to simulate single-phase flow phenomena in FHRs, whose baseline coolant, flibe, remains liquid up to 1430 deg. C at atmospheric pressure. The FHR advanced natural circulation analysis (FANCY) code has been developed specifically for FHR technology modeling and validation work has been performed by using CIET data in UCB. The code has been renamed General Flow (GenFlow) because it has been further developed to have the capability to model complicated pipework system based on object-oriented concept for user friendly purpose. Besides, it can be used to simulate both the force circulation as well as natural circulation. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Bridging Multiscale Processes in Earth System Models with Physics-Guided Hierarchical Machine Learning

Focal Area(s): The focal area of this whitepaper is, “Predictive modeling through the use of AI techniques and AI-derived model components; the use of AI and other tools to design a prediction system comprising of a hierarchy of models (e.g., AI driven model/component/parameterization selection).” The ideas and frameworks described herein are deemed site agnostic. As a use case, this group will initially focus on the coupling between land processes, surface/subsurface hydrological processes, coastal processes, and human activities in the U.S. Gulf states.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Intermittency in Precipitation: Duration, Frequency, Intensity, and Amounts Using Hourly Data

Abstract Intermittency is a core characteristic of precipitation, not well described by data and very poorly modeled. Detailed analyses are made of near-global gridded (about 1°) hourly or 3-hourly precipitation rates from two updated observational datasets [3-hourly TRMM 3B42, version 7, and hourly CMORPH, version 1.0, bias corrected (CRT)] and from special runs of CESM from January 1998 to December 2013 to obtain hourly values. The analyses explore the intermittency of precipitation: the frequency, intensity, duration, and amounts. A comparison is made for all products using several metrics with a focus on the duration of events, and a new metric is proposed based on the ratio of the frequency of precipitation at certain rates (0.2–2 mm h−1) for hourly versus 3-hourly versus daily data. For all seasons and rain rates, TRMM values are similar in pattern to CMORPH, but durations are about 80%–85%. It is mainly over land in the monsoons that CMORPH exceeds TRMM rain durations. Observed duration of precipitation events in CMORPH over oceans are 12–15 h in the tropics and subtropics, much less than the ~20 h for CESM. Hence, the observational results differ somewhat but both are considerably different from the model, which has too much precipitation overall, and it precipitates far too often at low rates and not enough for intense rates, with the divide about 1–2 mm h−1. There is a need to properly represent precipitation phenomena and processes either explicitly or implicitly (parameterized).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Transition times between the extremum points of the current–voltage characteristic of a resonant tunneling diode with hysteresis

A numerical solution to the problem of transient processes in a resonant tunneling diode featuring a current–voltage characteristic with hysteresis is found for the first time in the context of a coherent model (based on the coupled Schrödinger and Poisson equations) taking into account the Fermi distribution of electrons. The transitions from the high-current to the low-current state and vice versa, which result from the existence of hysteresis and are of great practical importance for ultrafast switches based on resonant tunneling diodes, are studied in detail. It is shown that the transition times for such processes initiated by the application of a small voltage can significantly exceed the characteristic time ℏ/Γ (where G is the width of the resonance level). It is established for the first time that the transition time can be reduced and made as short as the characteristic time ℏ/Γ by applying a sufficiently high voltage. For the parameters of the resonant-tunnelingdiode structure considered in this study, the required voltage is about 0.01 V.

36 MATERIALS SCIENCE↗
DOE OSTI2013

GPU accelerated fully space and time resolved numerical simulations of self-focusing laser beams in SBS-active media

A computer cluster equipped with Graphics Processing Units (GPUs) is used for simulating nonlinear optical wave packets undergoing Kerr self-focusing and stimulated Brillouin scattering in fused silica. We first recall the model equations in full (3+1) dimensions. These consist of two coupled nonlinear Schrödinger equations for counterpropagating optical beams closed with a source equation for light-induced acoustic waves seeded by thermal noise. Compared with simulations on a conventional cluster of Central Processing Units (CPUs), GPU-based computations allow us to use a significant (16 times) larger number of mesh points within similar computation times. Reciprocally, simulations employing the same number of mesh points are between 3 and 20 times faster on GPUs than on the same number of classical CPUs. Performance speedups close to 45 are reported for isolated functions evaluating, e.g., the optical nonlinearities. Since the field intensities may reach the ionization threshold of silica, the action of a defocusing electron plasma is also addressed.

97 MATHEMATICS AND COMPUTING↗