Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geometric modelling”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Effects of support and promoter on Ru catalyst activity in microwave-assisted ammonia synthesis

Microwave-assisted ammonia synthesis is a promising alternative to the energy-intensive Haber-Bosch process, specially at small- and medium-scale with renewable H 2 as resource. Here, we report that Cs promoted Ru/CeO 2 catalyst exhibits considerable activity at 533 K and ambient pressure. In this work, the combined theoretical and experimental approaches are adopted to optimize the electronic and geometric structures of Ru on the catalysts. Both DFT modeling work and structural characterization show that the strong interaction between Ru and CeO 2 results in the formation of highly dispersed Ru particles favoring ammonia synthesis. The higher electron donating ability of CeO 2 and lower electronegativity of Cs promoter result in higher electron density on Ru reducing the N≡N dissociation barrier. Finally, the work demonstrates the potential of microwave-assisted catalytic process in activating stable molecules for ammonia synthesis.

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

Molecular Dynamical and Quantum Mechanical Exploration of the Site-Specific Dynamics of Cy3 Dimers Internally Linked to dsDNA

Performing spectroscopic measurements on biomolecules labeled with fluorescent probes is a powerful approach to locating the molecular behavior and dynamics of large systems at specific sites within their local environments. The indocarbocyanine dye Cy3 has emerged as one of the most commonly used chromophores. The incorporation of Cy3 dimers into DNA enhances experimental resolution owing to the spectral characteristics influenced by the geometric orientation of excitonically coupled monomeric units. Various theoretical models and simulations have been utilized to aid in the interpretation of the experimental spectra. In this study, we employ all-atom molecular dynamics simulations to study the structural dynamics of Cy3 dimers internally linked to the dsDNA backbone. We used quantum mechanical calculations to derive insights from both the linear absorption spectra and the circular dichroism data. Furthermore, we explore potential limitations within a commonly used force field for cyanine dyes. The molecular dynamics simulations suggest the presence of four possible Cy3 dimeric populations. The spectral simulations on the four populations show one of them to agree better with the experimental signatures, suggesting it to be the dominant population. Furthermore, the relative orientation of Cy3 in this population compares very well with previous predictions from the Holstein–Frenkel Hamiltonian model.

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

Oregon Rivers Characterization Tool

DHS’s Regional Resiliency Assessment Program (RRAP) undertook the Oregon State Transportation Systems RRAP project between 2018-2021 to assess the seismic vulnerabilities of the state’s transportation system to a Cascadia Subduction Zone (CSZ) earthquake (CISA 2021). The project primarily focused on the effects of a CSZ earthquake on the mobility of emergency response efforts into the region. In the event that bridges located at crucial river crossings fail during a CSZ earthquake, rivers will become natural barriers affecting the ability to support response and recovery efforts. The characteristics of rivers (e.g., width, depth, and flow) are critical data to inform planning for the recovery of transportation systems across river barriers and determine reopening times following a CSZ earthquake. As part of the transportation system assessment, an Argonne National Laboratory (Argonne) project team developed a River Characterization Tool (RCT) to quantify the physical characteristics of the river system in Oregon, which can be used to inform a system-level analysis of transportation recovery and reopening times. This report mainly discusses approaches, methods, development, and applications of the RCT. The use of the river data generated from this tool for a system-level analysis of transportation recovery is discussed in the Oregon Transportation Systems RRAP Resiliency Assessment report (CISA 2021). The primary goal of the RCT is to generate statewide dataset of widths, depths, and flows for all rivers across Oregon to assess the effects of river barriers on emergency response and recovery efforts following a CSZ earthquake. The RCT leverages remote sensing data that multispectral sensors collect from satellites, employs a deep learning model to characterize river predictions, and computes river channel geometric dimensions with hydraulic relationships and river predictions. The entire process of river characterization is automated and implemented through the RCT.

58 GEOSCIENCES↗
DOE OSTI2022

Enhancing CO 2 Mineralization Rate and Extent of Iron and Steel Slag via Grinding

Roughly 10% of the CO 2 emissions from iron and steel making are attributable to the direct release of CO 2 from the thermal decomposition of carbonates to produce flux, mainly CaO, used for impurity removal. Notably, these direct emissions remain even if carbon-based steelmaking is replaced by hydrogen-based steelmaking. After removing impurities from the molten metal, this flux becomes the solid waste product called ‘slag’, a primarily Ca-silicate material. The transformation of slag back into carbonates is thermodynamically spontaneous with negative ΔG in the ambient environment, meaning that ~10% of the CO 2 emissions from iron and steel making could be negated if equipment and methods were developed to support CO 2 mineralization. However, the rate of CO 2 mineralization using slag is slowed by several environmental, geometric, and processing factors. We leverage an experimentally verified model of CO 2 mineralization to determine how to efficiently accelerate the process. Increasing the crystallinity of slag, increasing the relative humidity, and reducing the grain size of slag particles provide the greatest increase in CO 2 mineralization rate at the lowest energy penalty. Increasing the concentration of CO 2 and the temperature provide only modest increases in the CO 2 mineralization rate while incurring a substantial energy penalty. For steelmaking slags, CO 2 mineralization represents low-hanging fruit as the current reuse pathways are low value. For ironmaking slag, replacing the production of amorphous slag for the cement industry with the production of crystalline slag for CO 2 mineralization becomes financially preferable when a carbon price/tax exceeds 67.40 USD/t-CO 2 .

42 ENGINEERING↗
DOE OSTI2023

Design optimization of an ethanol heavy-duty engine using design of experiments and bayesian optimization

Diesel-fueled engines still hold a large market share in the medium and heavy-duty transportation sector. However, the increase in fossil fuel prices and the strict emission regulations are leading engine manufacturers to seek cleaner alternatives without a compromise in performance. Alcohol-based fuels, such as ethanol, offer a promising alternative to diesel fuel in meeting regulatory demands. Ethanol provides cleaner combustion and lower levels of soot due to its chemical properties, in particular its lower level of carbon content. In addition, the stoichiometric operating conditions of alcohol fueled engines enable the mitigation of NOx emissions in aftertreatment stage. With the promise of retrofitting diesel engines to run on ethanol to reduce emissions, the thermal efficiency of these engines remains the primary optimization target. In order to find the optimal ethanol-fueled engine design that maximizes the thermal efficiency, a large design space needs to be investigated using engineering tools. In this study, previous research by the authors on optimizing the design of a single-cylinder ethanol-fueled engine was extended to explore the design space for a heavy-duty multi-cylinder engine configuration. A heavy-duty engine setup with multiple operating conditions at different engine speeds and loads were considered. A design optimization analysis was performed to identify the potential designs that maximize the indicated thermal efficiency in an ethanol-fueled compression ignition engine. First, a computational fluid dynamics (CFD) model of the engine was validated using experimental data for four drive cycle points. Using a design of experiments (DoE) approach and a parameterized piston bowl geometry, the model was then exercised to explore the relationship among geometric features of the piston bowl and spray targeting angle and indicated thermal efficiency across all tested operating conditions. After evaluating 165~candidate designs, a piston bowl geometry was identified that yielded an increase between 1.3 to 2.2 percentage points in indicated thermal efficiency for all tested conditions, while satisfying the operational design constraints for peak pressure and maximum pressure rise rate. The increased performance was attributed to enhanced mixing that led to the formation of a more homogeneous distribution of in-cylinder temperature and equivalence ratio, higher combustion temperatures, and shorter combustion duration. Finally, a Bayesian optimization (BOpt) analysis was employed to find the optimal piston bowl geometry with a fixed spray injector angle for one of the operating conditions. Using BOpt, a piston candidate was identified that resulted in a 1.9~percentage point increase in thermal efficiency from the baseline design, yet only required 65\% of the design samples investigated using the DoE approach.

Tekgul, Bulut↗
DOE OSTI2025

Fluoride-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model Updates

This work presents the latest improvements to, and investigations performed with, the Fluoride-Cooled High-Temperature Pebble-Bed Reactor reference plant models for the United States Nuclear Regulatory Commission. These models, developed with the Comprehensive Reactor Analysis Bundle, or BlueCRAB, serve as the foundation for the future development of detailed design evaluation models based on license applications. BlueCRAB is the code suite proposed for non-light-water reactor systems safety analyses, and it incorporates various simulation tools developed by the Nuclear Energy Advance Modeling and Simulation program, including the Griffin code for reactor physics, the Pronghorn and SAM codes for core thermal fluids, the BISON code for solid conduction and fuel performance, and the SAM code for system analysis. The primary objective of this work is to assess the level of readiness of BlueCRAB for modeling fluoride-cooled high-temperature pebble-bed reactors. We first developed numerical models in BlueCRAB that include the key physics for this technology, ensuring an adequate level of fidelity for modeling the core performance during accident scenarios. This was followed by simulation of transient scenarios, two loss-of-forced-cooling events (one protected and one unprotected), and two control rod withdrawal events (one delayed and one prompt supercritical reactivity insertion). The analysis includes comparisons between the 2-D thermal fluid porous media models in Pronghorn and SAM, comparisons between coupled Pronghorn-Griffin and coupled SAM-Griffin models for two loss-of-forced cooling events and one control rod withdrawal event, and comparisons between SAM single-solve and domain-overlapping approaches for multi-scale thermal fluid coupling. In addition, we performed comparisons between 3-D, 2-D, and 0-D neutronic models for the two control rod withdrawal scenarios with Pronghorn-Griffin. The results show that the BlueCRAB models led to physically intuitive solutions for the scenarios examined. The changes in the various scalar and vector fields such as the neutron flux, power, temperatures, densities, pressures, and velocities are all within the expected ranges, and their distributions can be explained from the system response of the transients and the geometric and material variations. Several comparisons suggest that the porous media models in Pronghorn and SAM can lead to similar solutions, even though they are based on different methodologies. The simulations demonstrate that there are differences between the various levels of fidelity, and it is advisable to have flexible tools that can cover the breadth and depth of needs that may arise in future technical evaluations. We believe that BlueCRAB’s capabilities represent a significant asset for confirmatory analyses aimed at resolving important safety questions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Renormalization of the Standard Model Effective Field Theory from geometry

S-matrix elements are invariant under field redefinitions of the Lagrangian. They are determined by geometric quantities such as the curvature of the field-space manifold of scalar and gauge fields. We present a formalism where scalar and gauge fields are treated together, with a metric on the combined space of both types of fields. Scalar and gauge scattering amplitudes are given by the Riemann curvature Rijkl of this combined space, with indices i, j, k, l chosen to be scalar or gauge indices depending on the type of external particle. One-loop divergences can also be computed in terms of geometric invariants of the combined space, which greatly simplifies the computation of renormalization group equations. We apply our formalism to the Standard Model Effective Field Theory (SMEFT), and compute the renormalization group equations for even-parity bosonic operators to mass dimension eight.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Machine learning of 27Al NMR electric field gradient tensors for crystalline structures from DFT

NMR crystallography has emerged as a promising technique for the determination and refinement of atomic coordinates in crystal structures. The crystal structure of compounds containing quadrupolar nuclei, such as 27Al, can be improved by directly comparing solid-state NMR measurements to DFT computations of the electric field gradient (EFG) tensor. The non-negligible computational cost of these first-principles calculations limits the applicability of this method to all but the most well-defined structures. We developed a fast, low-cost machine learning model to predict EFG parameters based on local structural motifs and elemental parameters. We computed 8081 EFG tensors from 1681 27Al crystalline solids using DFT and benchmarked them against 105 experimentally measured 27Al sites. Surprisingly, simple local geometric features dominate the predictive performance of the resulting random-forest model, yielding an R2 value of 0.98 and an RMSE of 0.61 MHz for CQ, the quadrupolar coupling constant. This model accuracy should enable pre-refining future structural assignments before finally validating with first-principles calculations. Such a catalogue of 27Al NMR tensors can serve as a tool for researchers assigning complex NMR spectra influenced by the nuclear electric quadrupole interaction.

Sun, He↗
DOE OSTI2022

Manifold Learning-Based Polynomial Chaos Expansions for High-Dimensional Surrogate Models

In this work we introduce a manifold learning-based method for uncertainty quantification (UQ) in systems describing complex spatiotemporal processes. Our first objective is to identify the embedding of a set of high-dimensional data representing quantities of interest of the computational or analytical model. For this purpose, we employ Grassmannian diffusion maps, a two-step nonlinear dimension reduction technique which allows us to reduce the dimensionality of the data and identify meaningful geometric descriptions in a parsimonious and inexpensive manner. Polynomial chaos expansion is then used to construct a mapping between the stochastic input parameters and the diffusion coordinates of the reduced space. An adaptive clustering technique is proposed to identify an optimal number of clusters of points in the latent space. The similarity of points allows us to construct a number of geometric harmonic emulators which are finally utilized as a set of inexpensive pretrained models to perform an inverse map of realizations of latent features to the ambient space and thus perform accurate out-of-sample predictions. Thus, the proposed method acts as an encoder-decoder system which is able to automatically handle very high-dimensional data while simultaneously operating successfully in the small-data regime. The method is demonstrated on two benchmark problems and on a system of advection-diffusion-reaction equations which model a first-order chemical reaction between two species. In all test cases, the proposed method is able to achieve highly accurate approximations which ultimately lead to the significant acceleration of UQ tasks.

42 ENGINEERING↗
DOE OSTI2024

Deep-learning based artificial intelligence tool for melt pools and defect segmentation

Accelerating fabrication of additively manufactured components with precise microstructures is important for quality and qualification of built parts, as well as for a fundamental understanding of process improvement. Accomplishing this requires fast and robust characterization of melt pool geometries and structural defects in images. This paper proposes a pragmatic approach based on implementation of deep learning models and self-consistent workflow that enable systematic segmentation of defects and melt pools in optical images. Deep learning is based on an image-to-image translation–conditional generative adversarial neural network architecture. An artificial intelligence (AI) tool based on this deep learning model enables fast and incrementally more accurate predictions of the prevalent geometric features, including melt pool boundaries and printing-induced structural defects. We present statistical analysis of geometric features that is enabled by the AI tool, showing strong spatial correlation of defects and the melt pool boundaries. The correlations of widths and heights of melt pools with dataset processing parameters show the highest sensitivity to thermal influences resulting from laser passes in adjacent and subsequent layer passes. The presented models and tools are demonstrated on the aluminum alloy and datasets produced with different sets of processing parameters. However, they have universal quality and could easily be adapted to different material compositions. The method can be easily generalized to microstructural characterizations other than optical microscopy.

additive manufacturing↗
DOE OSTI2020

Aeroelastic oscillations of a pitching flexible wing with structural geometric nonlinearities: Theory and numerical simulation

In this paper, we focus on the derivation of an analytical model of the aeroelastic dynamics of an elastically mounted flexible wing. This equations of motion obtained serve to help understand the behaviour of the aeroelastic wind tunnel setup in question, which consists of a rectangular wing with a uniform NACA 0012 airfoil profile, whose base is free to rotate rigidly about a longitudinal axis. Of particular interest are the structural geometric nonlinearities primarily introduced by the coupling between the rigid body pitch degree-of-freedom and the continuous system. A coupled system of partial differential equations (PDEs) coupled with an ordinary differential equation (ODE) describing axial-bending-bending-torsion-pitch motion is derived using Hamilton's principle. A finite dimensional approximation of the system of coupled differential equations is obtained using the Galerkin method, leading to a system of coupled nonlinear ODEs. Subsequently, these nonlinear ODEs are solved numerically using Houbolt's method. The results that are obtained are verified by comparison with the results obtained by direct integration of the equations of motion using a finite difference scheme. Adopting a linear unsteady aerodynamic model, it is observed that the system undergoes coalescence flutter due to coupling between the rigid body pitch rotation dominated mode and the first flapwise bending dominated mode. Finally, the behaviour of the limit cycle oscillations is primarily influenced by the structural geometric nonlinear terms in the coupled system of PDEs and ODE.

42 ENGINEERING↗
DOE OSTI2024

Multilevel Monte Carlo methods for the Grad-Shafranov free boundary problem

The equilibrium configuration of a plasma in an axially symmetric reactor is described mathematically by a free boundary problem associated with the celebrated Grad-Shafranov equation. The presence of uncertainty in the model parameters introduces the need to quantify the variability in the predictions. This is often done by computing a large number of model solutions on a computational grid for an ensemble of parameter values and then obtaining estimates for the statistical properties of solutions. In this study, we explore the savings that can be obtained using multilevel Monte Carlo methods, which reduce costs by performing the bulk of the computations on a sequence of spatial grids that are coarser than the one that would typically be used for a simple Monte Carlo simulation. We examine this approach using both a set of uniformly refined grids and a set of adaptively refined grids guided by a discrete error estimator. Numerical experiments show that multilevel methods dramatically reduce the cost of simulation, with cost reductions typically on the order of 60 or more and possibly as large as 200. Furthermore, adaptive griding results in more accurate computation of geometric quantities such as x-points associated with the model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Deformed Fredkin model for the ν = 5 / 2 Moore-Read state on thin cylinders

We propose a frustration-free model for the Moore-Read quantum Hall state on sufficiently thin cylinders with circumferences ≲ 7 magnetic lengths. While the Moore-Read Hamiltonian involves complicated long-range interactions between triplets of electrons in a Landau level, our effective model is a simpler one-dimensional chain of qubits with deformed Fredkin gates. We show that the ground state of the Fredkin model has high overlap with the Moore-Read wave function and accurately reproduces the latter's entanglement properties. Moreover, we demonstrate that the model captures the dynamical response of the Moore-Read state to a geometric quench, induced by suddenly changing the anisotropy of the system. We elucidate the underlying mechanism of the quench dynamics and show that it coincides with the linearized bimetric field theory. The minimal model introduced here can be directly implemented as a first step towards quantum simulation of the Moore-Read state, as we demonstrate by deriving an efficient circuit approximation to the ground state and implementing it on an IBM quantum processor. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

304L Can Crush Validation Studies

Accurate prediction of ductile behavior of structural alloys up to and including failure is essential in component or system failure assessment, which is necessary for nuclear weapons alteration and life extensions programs of Sandia National Laboratories. Modeling such behavior requires computational capabilities to robustly capture strong nonlinearities (geometric and material), rate- dependent and temperature-dependent properties, and ductile failure mechanisms. This study's objective is to validate numerical simulations of a high-deformation crush of a stainless steel can. The process consists of identifying a suitable can geometry and loading conditions, conducting the laboratory testing, developing a high-quality Sierra/SM simulation, and then drawing comparisons between model and measurement to assess the fitness of the simulation in regards to material model (plasticity), finite element model construction, and failure model. Following previous material model calibration, a J 2 plasticity model with a microstructural BCJ failure model is employed to model the test specimen made of 304L stainless steel. Simulated results are verified and validated through mesh and mass-scaling convergence studies, parameter sensitivity studies, and a comparison to experimental data. The converged mesh and degree of mass-scaling are the mesh discretization with 140,372 elements, and a mass scaling with a target time increment of 1.0e-6 seconds and time step scale factor of 0.5, respectively. Results from the coupled thermal-mechanical explicit dynamic analysis are comparable to the experimental data. Simulated global force vs displacement (F/D) response predicts key points such as yield, ultimate, and kinks of the experimental F/D response. Furthermore, the final deformed shape of the can and field data predicted from the analysis are similar to that of the deformed can, as measured by 3D optical CMM scans and DIC data from the experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Stochastic parametric skeletal dosimetry model for humans: Pediatric and adult computational skeleton phantoms for internal bone marrow dosimetry

Currently, computational phantoms that simulate skeletal tissues are used in active red bone marrow (AM) internal dosimetry. Up-to-date reference computational phantoms recommended by the ICRP are based on the analysis of CT-images of cadavers. Such phantoms have significant disadvantages. One disadvantage is that the assessment of uncertainty due to the population variability of skeleton dimensions and microstructure results from the limited availability of autopsy material. Another disadvantage is the simplified modelling of cortical layer and bone microarchitecture. A method of stochastic parametric skeletal dosimetry modelling of the bone structures – SPSD modelling – has been developed as an alternative to the ICRP reference phantoms. In the framework of this approach, skeletal phantom parameters are evaluated based on extensively reviewed results of published measurements of real bones. The SPSD approach allows for the assessment of both population-average values and their variability. SPSD-phantoms of the skeleton are modelled in voxel representation. They consist of smaller phantoms of the bone sites – segments – described by simple geometric shapes with uniform microarchitecture parameters. Such segmentation makes it possible to account for non-homogeneous skeletal microarchitecture and to model the bone structure with the required voxel resolution to elaborate suitable skeletal phantoms. The current study presents the parameters of the SPSD skeletal phantoms for the following age-groups: newborn, 1-year-old, 5-year-old, 10-year-old, 15-year-old (male and female), and adult (male and female). This skeletal phantom can be used for dosimetry as an alternative to available reference phantoms for bone-seeking radionuclides. The above-mentioned age- and sex-specific skeletal phantoms are comprised of 289 unique segments. The characteristics of the SPSD phantoms do not contradict published data and are in good agreement with the measurement results of real bones.

Science & Technology - Other Topics↗
DOE OSTI2021

Steepest-descent algorithm for simulating plasma-wave caustics via metaplectic geometrical optics

The design and optimization of radiofrequency-wave systems for fusion applications is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at wave cutoffs and caustics. To accurately model the wave behavior in these regions, more advanced and computationally expensive “full-wave” simulations are typically used, but this is not strictly necessary. A new generalized formulation called metaplectic geometrical optics (MGO) has been proposed that reinstates GO near caustics. The MGO framework yields an integral representation of the wavefield that must be evaluated numerically in general. We present an algorithm for computing these integrals using Gauss-Freud quadrature along the steepest-descent contours. Benchmarking is performed on the standard Airy problem, for which the exact solution is known analytically. Furthermore, the numerical MGO solution provided by the new algorithm agrees remarkably well with the exact solution and significantly improves on previously derived analytical approximations of the MGO integral.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Anomalous Hall crystals in rhombohedral multilayer graphene. II. General mechanism and a minimal model

Here we propose a minimal "three-patch model"for the anomalous Hall crystal (AHC), a topological electronic state that spontaneously breaks both time-reversal symmetry and continuous translation symmetry. The proposal for this state is inspired by the recently observed integer and fractional quantum Hall states in rhombohedral multilayer graphene at zero magnetic field. There, interaction effects appear to amplify the effects of a weak moiré potential, leading to the formation of stable, isolated Chern bands. It has been further shown that Chern bands are stabilized in mean-field calculations even without a moiré potential, enabling a realization of the AHC state. Our model is built on the dissection of the Brillouin zone into patches centered around high-symmetry points. Within this model, the wave functions at high-symmetry points fully determine the topology and energetics of the state. We extract two quantum geometrical phases of the noninteracting wave functions that control the stability of the topologically nontrivial AHC state. The model predicts that the AHC state wins over the topological trivial Wigner crystal in a wide range of parameters, and agrees very well with the results of full self-consistent Hartree-Fock calculations of the rhombohedral multilayer graphene Hamiltonian.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Optimal design of microseismic monitoring network: Synthetic study for the Kimberlina CO2 storage demonstration site

Microseismic monitoring can play a crucial role to ensure safe long-term geological carbon storage. For reliable long-term monitoring for CO2-injection-induced microseismic events, a surface seismic array is desirable in addition to borehole geophone sensors. Optimal design of surface seismic network is of great interest to achieve cost-effective monitoring. We develop a methodology to determine the optimal number of surface seismic stations with a geometrically satisfactory distribution for given monitoring regions. We design an optimal microseismic monitoring network based on widely-accepted guiding principles, and the relationship between the location accuracy of microseismic events and the total number of seismic stations. We determine the optimal number of seismic stations based on the trade-off curve of the event location accuracy vs. the total number of seismic stations. We apply our optimal design method to the Kimberlina carbon storage site in California. We use a synthetic Kimberlina model to show that approximately 20 surface seismic stations with a geometrically satisfactory distribution are preferred for the best trade-off between the cost and the event location accuracy.

58 GEOSCIENCES↗