SEARCH · Search NASA
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.
Amesh generates discretized grids for numerical modeling of flow and transport problems in which the formulation is based on the integral finite difference method (IFDM). For example, the output of Amesh can be used directly as input to the TOUGH2 or TOUGH3 numerical simulators. Amesh can generate 1D, 2D, or 3D numerical grids for a given set of locations, i.e., the centers of each discrete sub-domain. In the 2D areal plane Voronoi tessellation is used.
APBuilder is a software tool designed to generate atmospheric profiles for use with AC2Dr by downloading weather model data and transforming it into 1D or 2D binary profiles. These binary profiles serve as required inputs for AC2Dr, an LLNL-developed, open source, two-dimensional numerical solver for the acoustic wave equation based on the finite difference method.
Seismic-H2 is an integrated software package for geological hydrogen reservoir simulation, optimization, and leakage monitoring. The package includes multiple components: (1) code used for modeling seismic wave propagation in 3D heterogeneous elastic media based on finite-difference method to support detection of geological hydrogen storage reservoir leakage; (2) 3D reservoir simulations of leaks from an underground reservoir and 3D simulations of saline aquifers and depleted gas reservoirs; (3) seismic monitoring costs of passive and active seismic monitoring required for UHS; (4) rock physics calculations and interpolations for converting the reservoir simulations from part (2) into the elastic media models in part (1); (5) pre-processing seismic data; and lastly (6), a GUI interface that combines these different components.
SAND2025-14430O Hydrogen Diffusion Through Stainless Steel is a tool that employs a finite difference method algorithm to solve Fick's law and other mass transport models pertinent to the diffusion of hydrogen through stainless steel. Additionally, the software uses particle swarm optimization to determine physical parameters based on experimental data. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.
The complex Ginzburg-Landau equation with additive noise is a stochastic partial differential equation that describes a remarkably wide range of physical systems which include coupled non-linear oscillators subject to external noise near a Hopf bifurcation instability and spontaneous structure formation in non-equilibrium systems, e.g., in cyclically competing populations or oscillatory chemical reactions. Here, we employ a finite-difference method to numerically solve the noisy complex Ginzburg-Landau equation on a two-dimensional domain with the goal to investigate its non-equilibrium dynamics when the system is quenched into the “defocusing spiral quadrant”. We observe slow coarsening dynamics as oppositely charged topological defects annihilate each other, and characterize the ensuing aging scaling behavior. We conclude that the physical aging features in this system are governed by non-universal aging scaling exponents. We also investigate systems with control parameters residing in the “focusing quadrant”, and identify slow aging kinetics in that regime as well. Finally, we provide heuristic criteria for the existence of slow coarsening dynamics and physical aging behavior in the complex Ginzburg-Landau equation.
The seismic simulation code SW4 was adapted to run on modern supercomputers. SW4 solves the elastic wave equation in second order formulation. It uses a fourth order accurate finite difference method with explicit time stepping. The approximation conserves energy perfectly, and is stable under a CFL restriction on the time step.
This is the final scientific report for the award DE-FE0023305, entitled “Quantitative Characterization of Impacts of Coupled Geomechanics and Flow on Safe and Permanent Geological Storage of CO 2 in Fractured Aquifers.” The work has been divided into six tasks. In Task 2, we characterized rock properties, which are important when developing a quantitative approach for understanding and predicting geomechanical effects on large-scale CO 2 injection and long-term storage in the subsurface. Rock properties of interest for this characterization include porosity, permeability, elastic constants, strength, and heat capacity. We measured rock properties for three different rock types: concrete, sandstone and shale. These properties were obtained from the acoustic measurement, permeability and porosity measurement, Brazilian test, the uniaxial compressive test, and heat capacity measurement. Then, using the Brazilian technique to fracture samples, we measured their permeability under brine injection for varying effective stresses, followed by similar measurements of permeability under sc-CO 2 injection. Permeability to brine and sc-CO 2 decreased as the effective stress increased. The apparent permeability to sc-CO 2 was an order of magnitude less than that for brine, a consequence of CO 2 being the non-wetting fluid. In Task 3, we developed understandings and correlations for CO 2 injection pressure induced fracturing. We initially fractured four eight inch cubes of concrete with brine to establish a standard test procedure. Then, additional concrete samples were fractured using brine. Confining stresses were exerted on the samples and brine was injected at a constant rate into a borehole that was previously drilled into the concrete cube. The borehole pressure was measured and fracture initiation was identified as a peak in the pressure profile. That was followed by fracture propagation, which was identified as a plateau in the pressure profile, and finally the fracture reached the sample surface. The sample surfaces were photographed before and after fracturing, which enabled determination of where the fracture reached the sample surface. Afterwards, aqueous food dye solution was injected into the fractured concrete samples to color the fracture planes. The borehole was filled with food dye solution in advance, a gas pulse was transmitted to the borehole, and the borehole was then shut in to allow the food dye to be driven into existing fractures. Then, the concrete samples were broken down by high-pressure gas injection to reveal the geometry and morphology of the hydraulic fractures. In addition, acoustic measurements were conducted on multiple locations on each face of the concrete samples before and after injection for comparison. This comparison was used to confirm changes in the internal structure of the concrete sample, namely the formation of fractures. There were six concrete samples fractured with brine. Afterwards, concrete samples were fracture using sc-CO 2 . The main change in experimental procedure was the inclusion of temperature control, since in the field, we mostly encounter temperature and pressure conditions above the supercritical point of CO 2 . The concrete samples were pre-heated in an oven to elevate their temperature, and thermal tape was used to heat the CO 2 to a temperature above its supercritical temperature. Two concrete samples were fractured using liquid CO 2 , four concrete samples were fractured using gaseous CO 2 , and twenty two concrete samples were fractured using sc-CO 2 . For the concrete samples fractured using sc-CO 2 , three had pre-existing hydraulic fractures on the surface near the wellbore, and two were water saturated (as opposed to containing only air). Two concrete samples were composite samples, with a high permeability ball in the center, in order to represent CO 2 injection into a high permeability zone surrounded by a low permeability sealing formation. Later, we fractured five shale samples, obtained from the Niobrara shale outcrop, with sc-CO 2 . The injection rate was constant, confining stresses were applied, the shale samples were pre-heated in an oven to the desired temperature, and the sample was then fractured. An analysis of the failure type was done for many of the concrete samples, based on the predicted break down pressure for tensile and shear failure, and it was found that most samples fractured due to shear failure. Also, the presence of induced fractures has little to no effect on the fracture initiation pressure, but significantly changes the fracture morphology. The CO 2 fracturing in water saturated samples behaved much similar to brine induced fracturing. There was a significant difference of breakdown pressure for injection of CO 2 and brine, where the breakdown pressure for CO 2 induced fracturing is generally around the minimum horizontal stress and that for brine is much higher. Finally, the fracture orientation for most samples with large confining stress differential was perpendicular to the minimum horizontal stress. When the stress differential, especially the difference between the two directions with the smallest stresses was small, the fracture orientation was determined not only by the minimum stress direction, but also the magnitude of the confining stress, the tensile strength of the rock, and the breakdown pressure. In Task 4, we modified our coupled flow-geomechanical models to model fracture growth and propagation in storage formations and caprocks. These flow-geomechanical models are TOUGH2-CSM and TOUGH2-FLAC. The TOUGH2-CSM fluid and heat flow formulation is based on the TOUGH2 formulation of mass and energy conservation equations that govern fluid and heat flow in general multiphase, multicomponent, multi-porosity systems. The TOUGH2-CSM simplified geomechanical formulation is based on the linear theory of elasticity applied to multi-porosity non-isothermal (thermo-multi-poroelastic) media. We previously derived, from the fundamental relations of the linear theory of elasticity, an equation relating mean stress, pore pressures, temperatures, and body force (the Mean Stress equation) that we added to the TOUGH2-CSM fluid and heat flow equations. We extended the TOUGH2-CSM simplified geomechanical formulation to calculate the entire stress tensor by deriving equations for stress tensor components from derivatives of the Cartesian thermo-multi-poroelastic Navier equation components. In addition, we derived equations for stress tensor components in rz-coordinates. We solved these geomechanical equations using the integral finite difference method. This method utilizes momentum fluxes obtained by “factoring” a divergence operator out of each geomechanical equation. In Task 5, we did a literature survey to determine suitable correlations for fracture initiation, growth, and propagation. After surveying numerous papers, we modified TOUGH2-CSM to model stress dependent fracture initiation and growth using the Mohr-Coulomb criterion for shear failure in faults and caprock and the condition of negative minimum effective stress for tensile failure. Tensile failure results in a fracture whose width depends on the difference between fracture pressure and minimum stress and whose propagation is based on the difference between the stress intensity factor at the fracture tip and the rock toughness. The TOUGH-FLAC model is based on the FLAC3D ubiquitous joint model, extended for the modeling of permeability changes induced by changes in effective normal stress across fractures as well permeability changes caused by shearing along fractures. In Task 6, we validated our coupled flow-geomechanical models using analytical solutions and problems from the literature. We validated the TOUGH2-CSM modifications using analytical solutions to the displacement from a uniform load on a semi-infinite elastic medium and the two-dimensional Mandel-Cryer effect. Those were followed by a sample problem to test our stress tensor calculations, which was for injection into a single-phase reservoir with constant properties. Comparisons of our simulator to published results were done on for the depletion of a single-phase reservoir with stress dependent porosity, two examples from the In Salah gas project, an axisymmetric baserock–reservoir–caprock system with a normal faulting stress regime that had thermally induced shear stresses, simulations of CO 2 leakage through caprock in a two-dimensional reservoir in Cartesian coordinates, and predictions of caprock failure from CO 2 injection into an axisymmetric reservoir. Finally, we simulated CO 2 injection pressure-induced fracturing from one of our laboratory studies, Sample 39, and obtained a fracture profile. We also validated the TOUGH2-FLAC model for fluid driven fracture growth against solutions based on the KGD model and a case with an inclined fracture that was loaded from the boundaries to achieve wing-crack propagation. The strain softening tensile behavior and softening of modulus considering a simple damage approach was verified by a simulation representing hydraulic fracturing stress measurement around a vertical well. Finally, the model was validated against deep fracture zone opening and surface uplift at In Salah with inverses analysis used to improve the match of simulation with field data. In Task 7, we developed a scheme based on inverse modeling that can be used to identify caprock leakage. Injection of fluid into a reservoir results in a time varying pressure profile that depends on the properties of the rock and the fluid. The presence of significant caprock leakage would effect this pressure profile and we used inverse modeling, namely the Levenberg-Marquardt method, to determine the leakage location from the pressure profile. A two-dimensional and a three-dimensional example based on published data were presented. The published simulation was run and its results were considered as “field” data. In the two-dimensional example, a simulation of CO 2 injection into a two-dimensional aquifer-caprock system, the caprock had a vertical fault through which leakage could occur. The location of this leakage was obtained from our inverse modeling scheme. In the three-dimensional example, the simulation domain had four geological layers, one of which was the injection zone that contained a horizontal injection well. We simulated fractures as high permeability gridblocks, introduced a fracture that spanned the caprock depth, and ran a simulation to obtain “field” data. Then, using our inverse modeling scheme, we obtained the location of that fracture.
Accurate modeling of gas hydrate reservoir productivity and geomechanical risks associated with subsurface dissociation of natural gas hydrates (NGH) requires the determination of model parameters through physical testing on natural hydrate-bearing sediments (HBS). This involves investigating the hydro-mechanical behavior of undisturbed hydrate samples from nature under in situ conditions using pressure core characterization and analysis, which provides a unique opportunity for research. By employing state-of-the-art micro computed tomography imagery on cryogenically preserved, hydrate-bearing sediment samples, we can determine hydrate saturation as well as permeability with and without the presence of hydrates in the sediment. Furthermore, utilizing a machine learning based image segmentation technique, it is possible to extract pore space and grain information. Subsections of the entire image volume were used to determine anisotropic permeabilities using a finite-difference method Stokes solver (FDMSS). Additionally, permeability measurements on whole pressure and temperature preserved hydrate-bearing core were analyzed by utilizing the National Energy Technology Laboratory’s (NETL) Pressure Core Characterization and X-ray CT Visualization Tool (PCXT) to manipulate, cut, and analyze pressure preserved sediment. Permeabilities were measured under a broad range of vertical stress states to simulate expected pressure changes during production scenarios, and the results show that permeabilities derived from images are in agreement with those from traditional core derived experiments. The collected stress-dependent permeability, permeability anisotropy, and corresponding gas hydrate saturations provide valuable input into numerical simulations of reservoir productivity. These properties have been proven to be key parameters determining a long-term reservoir response under depressurization.
The objective of this technical report is to describe the Courant–Friedrichs–Levy (CFL) limit for stability, as it applies to the solution of the Hamilton–Jacobi–Bellman equation with explicit finite difference methods.
This dataset consists of soil temperature time series that were used to estimate soil thermal diffusivity and its uncertainty trough the probabilistic modelling approach developed and presented in the article named "Probabilistic estimation of depth-resolved profiles of soil thermal diffusivity from temperature time series" and published in Earth Surface Dynamics. There are two compressed (.zip) files that contains synthetic (Synthetic_soiltemp_Teller.zip) and field (Field_soiltemp_Teller.zip) data. There is one MATLAB file that requires MATLAB to execute but any text editor can open it. The Synthetic_soiltemp_Teller.zip file includes 5 comma-delimited data files (.csv) each of which contains soil temperature time series generated through forward modeling (i.e., heat-conduction process in a heterogeneous medium using an explicit finite difference method) to mimic various types of temperature gradients, trend and fluctuations. These synthetic soil temperatures were used to investigate the impact of different environmental conditions on the uncertainty of thermal diffusivity estimates. The Field_soiltemp_Teller.zip file contains 28 comma-delimited data files (.csv) out of which (a) 27 files includes soil temperature time series recorded from 27 temperature probes located in a site along Teller Road about 40 km northwest of Nome, Alaska (64.72°N, 165.94°W), (b) one includes the name and coordinates of the 27 probes. These field soil temperatures were used to infer soil thermal diffusivity at numerous locations and depths in a discontinuous permafrost environment, and to evaluate the links between the estimated soil thermal diffusivity values and soil physical properties. The comma-delimited data files (.csv) of the synthetic and field soil temperature time series includes date and time (UTC) in the first column and soil temperature from 5 cm below the ground surface to 1.05 m depth (with 5 or 10 cm spacing between sensors) in the other columns. The measurements were acquired every 15 minutes. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).
The mechanism of perpendicular momentum input of lower hybrid waves and its influence on plasma rotation are studied. Discussion for parallel momentum input of lower hybrid waves is presented for comparison. It is found out that both toroidal and poloidal projections of perpendicular momentum input of lower hybrid waves are stronger than those of parallel momentum input. The perpendicular momentum input of lower hybrid waves therefore plays a dominant role in forcing the changes of rotation velocity observed during lower hybrid current drive. Lower hybrid waves convert perpendicular momentum carried by the waves into the momentum of dc electromagnetic field by inducing a resonant-electron flow across flux surfaces therefore charge separation and a radial dc electric field. The dc field releases its momentum into plasma through the Lorentz force acting on the radial return current driven by the radial electric field. Plasma is spun up by the Lorentz force. An improved quasilinear theory with gyro-phase dependent distribution function is developed to calculate the radial flux of resonant electrons. Rotations are determined by a set of fluid equations for bulk electrons and ions, which are solved numerically by applying a finite-difference method. Analytical expressions for toroidal and poloidal rotations are derived using the same hydrodynamic model.
A diffusion-based and a discrete ordinates method are used to simulate a neutron noise experiment in the research reactor AKR-2 at the Technical University in Dresden, Germany. The AKR-2 reactor provides an interesting case for the comparison between the two methods because it is characterized by large heterogeneities and regions with low macroscopic neutron cross-sections. For the calculations, the same spatial discretization and the same set of two-energy macroscopic neutron cross-sections with isotropic scattering are used. Significant discrepancies between the diffusion-based and discrete ordinates methods are found in regions of the systems where the diffusion approximation is expected to be inaccurate in reproducing characteristics of the static neutron flux and neutron noise. (authors)
We present computational studies of quasi three-dimensional nanowell (NW) and nanopost (NP) plasmonic crystals for applications in surface enhanced Raman spectroscopy (SERS). The NW and NP plasmonic crystals are metal coated arrays of cylindrical voids or posts, respectively, in a dielectric substrate characterized by a well/post diameter ( D ), relief depth ( RD ), periodicity ( P ), and metal thickness ( MT ). Each plasmonic crystal is modeled using the three-dimensional finite-difference time-domain (FDTD) method with periodic boundary conditions in the x- and y-directions applied to a computational unit cell to simulate the effect of a periodic array.
Explore the source record for details and available documents.
Abstract We present a general numerical solution method for control problems with state variables defined by a linear PDE over a finite set of binary or continuous control variables. We show empirically that a naive approach that applies a numerical discretization scheme to the PDEs to derive constraints for a mixed-integer linear program (MILP) leads to systems that are too large to be solved with state-of-the-art solvers for MILPs, especially if we desire an accurate approximation of the state variables. Our framework comprises two techniques to mitigate the rise of computation times with increasing discretization level: First, the linear system is solved for a basis of the control space in a preprocessing step. Second, certain constraints are just imposed on demand via the IBM ILOG CPLEX feature of a lazy constraint callback. These techniques are compared with an approach where the relations obtained by the discretization of the continuous constraints are directly included in the MILP. We demonstrate our approach on two examples: modeling of the spread of wildfire and the mitigation of water contamination. In both examples the computational results demonstrate that the solution time is significantly reduced by our methods. In particular, the dependence of the computation time on the size of the spatial discretization of the PDE is significantly reduced.
The standard particle-in-cell (PIC) method employs explicit finite-difference (FD) methods (e.g. the leap-frog scheme) for both spatial and temporal integrations. Here, we employ a pseudospectral method for solving the Poisson equation and a fully implicit time integration to achieve exact energy conservation. The advantage of a pseudospectral field solver is its spectral accuracy in solving field solutions. Earlier studies of implicit time integration of PIC FD equations can enforce exact energy exchange between field and particles, resulting in exact energy-conserving schemes. Here, we prove that the exact energy conservation property can be carried over to the pseudospectral scheme. Simultaneously, we provide a solution to ensure a pseudospectral charge continuity equation. We demonstrate the new scheme in a 2D electrostatic PIC code. In conclusion, theoretical results are confirmed via numerical examples.
Laser powder bed fusion (LPBF) additive manufacturing makes near-net-shaped parts with reduced material cost and time, rising as a promising technology to fabricate Ti-6Al-4V, a widely used titanium alloy in aerospace and medical industries. However, LPBF Ti-6Al-4V parts produced with 67° rotation between layers, a scan strategy commonly used to reduce microstructure and property inhomogeneity, have varying grain morphologies and weak crystallographic textures that change depending on processing parameters. Here, this study predicts LPBF Ti-6Al-4V solidification at three energy levels using a finite difference-Monte Carlo method and validates the simulations with large-area electron backscatter diffraction (EBSD) scans. The developed model accurately shows that a <001> texture forms at low energy and a <111> texture occurs at higher energies parallel to the build direction but with a lower strength than the textures observed from EBSD. A validated and well-established method of combining spatial correlation and general spherical harmonics representation of texture is developed to calculate a difference score between simulations and experiments. The quantitative comparison enables effective fine-tuning of nucleation density (N 0 ) input, which shows a nonlinear relationship with increasing energy level. Future improvements in texture prediction code and a more comprehensive study of N 0 with different energy levels will further advance the optimization of LPBF Ti-6Al-4V components. These developments contribute a novel understanding of crystallographic texture formation in LPBF Ti-6Al-4V, the development of robust model validation and calibration pipeline methodologies, and provide a platform for mechanical property prediction and process parameter optimization.
Not Available