Abstract: A computationally efficient general method of modeling or simulating matrix acidizing treatment when flow is not axisymmetric involves determining streamlines in the general flow field using complex potential theory to solve for the flow along the streamlines. The flow over a time step is used to model the propagation of the acid front and the creation and extension of wormholes.

Abstract: Various embodiments of systems and methods are described in which system simulation techniques are combined with graphical programming techniques in a common environment. For example, various embodiments of the methods comprise displaying a graphical data flow diagram connected to a system diagram, e.g., where the graphical data flow diagram and the system diagram are displayed together in a compound diagram. In the displayed compound diagram, the graphical data flow diagram may be connected to the system diagram, e.g., by a line or wire. In one embodiment the wire may visually indicate that the graphical data flow diagram is executable to produce a value that is provided as an input signal for a simulation performed based on the system diagram. In another embodiment the wire may visually indicate that the graphical data flow diagram receives an output value from the system simulation as input.

Abstract: Disclosed are techniques for representing and modeling one or more systems in which each system corresponds to an application mode. This may be done for one or more geometries using local and/or non-local couplings. For each application mode, physical quantities are modeled and may be defined using a graphical user interface. Physical properties may be used to model the physical quantities of each system. The physical properties may be defined in terms of numerical values or constants, and mathematical expressions that may include numerical values, space coordinates, time coordinates, and actual physical quantities. Physical quantities and any associated variables may apply to some or all of a geometric domain, and may also be disabled in other parts of a geometrical domain. Partial differential equations describe the physical quantities. One or more application modes may be combined using an automated technique into a combined system of partial differential equations as a multiphysics model.

Abstract: A method and system for optimal design of a water distribution network is provided. Three levels of optimization are available as options for selection by the user and these include least cost optimization, maximum benefit optimization and cost-benefit tradeoff optimization. The optimization models, in accordance with the present invention, include solutions generated by a competent genetic algorithm that can take into account multiple objective functions. Pareto-optimal solutions are produced for the whole range of a budget for the water distribution system. Tradeoff solutions allow engineers to apply engineering judgement to choose the true optimal solution under the specific application being considered.

Abstract: Disclosed is a method of designing an aerosol spray dispensing system. One or more to-be-determined attributes of the spray dispenser system are identified, including properties of a fluid to be dispensed, characteristics of the dispensing mechanism, and performance characteristic values. The identified attributes are determined using equations taking into account relationships among the properties of the fluid to be dispensed, the characteristics of the dispensing mechanism and the performance characteristic values. Once determined, the attributes are used to design the spray dispensing system.

Abstract: Embodiments of the present invention provide a finite difference method, a finite element projection, and a non-finite-element implementation of the finite element projection for solving the incompressible continuity equation on quadrilateral grids to obtain pressure and fluid velocity. An exemplary apparatus for simulating and analyzing the flow of an incompressible fluid is also provided. More particularly, embodiments of the present invention provide a non-finite-element-method (non-FEM) to implement the finite element projection to obtain a linear system whose coefficient matrix has an average bandwidth of five, in comparison to the finite difference scheme using conventional finite element projection and implementation with a bandwidth of nine. A two-dimensional-five-point-stencil-based non-FEM discretization is formulated to enforce the incompressible continuity equation.

Abstract: Method for treating a moving surface during a transient condition in a computational fluid dynamics system multi-cell computer model, the model having a fixed underlying geometric grid, includes identifying cells in the grid adjacent the moving surface at the beginning of a transient time step and identifying beginning adjacent cell boundaries formed in the grid by the moving surface. The method also includes computing effective thermophysical value fluxes through the beginning cell boundaries in accordance with the sense of the movement of the moving surface relative to the beginning adjacent cells. The method further includes solving transport equations for all the cells to determine cell thermophysical values at the end of the time step, using the computed effective fluxes and reforming the grid to comply with the moving surface at the end of the time step. Apparatus includes a digital computer programmed to carry out the above method.

Abstract: A method for quantifying uncertainties related to continuous and discrete (qualitative) parameters descriptive of a medium such as an underground zone and/or for managing the selection of a scenario from a series of possible scenarios relative to this medium, by construction of experiment designs and results analysis suited to the experiment designs constructed.

Abstract: A computer implemented method for simulating a final pattern of a droplet of a fluid having a plurality of fluid properties is disclosed. The method includes representing a position of a moving fluid boundary at a first point in time with a plurality of fluid markers, wherein the moving fluid boundary provides a boundary for a single-phase fluid. The moving fluid boundary separates a simulation space into a fluid space and a vacuum space, and the single-phase fluid inhabits the fluid space. The method further includes evaluating a plurality of indicator function defined on a quadrilateral grid for a velocity field, a pressure, and the plurality of fluid markers. The quadrilateral grid is formed by a plurality of points inside the fluid space and the vacuum space.

Abstract: A numerical model related to fluttering of an insect, when an equivalent model of actual structure of a wing of the insect is moved in the air in accordance with a model of fluttering motion of the wing of the insect is calculated by fluid-structure interactive analysis, in which behavior of the wing and behavior of the air are given as numerical values, including interaction therebetween. Thereafter, a method of controlling a fluttering robot, wing shape and the like are determined by modifying numerical models of fluttering of an insect prepared by fluid-structure interactive analysis, in accordance with sensitivity analysis. Accordingly, a method of preparing numerical models of wing and air considering the behavior of the wing of the insect in the air is provided and, in addition, a method of manufacturing a fluttering robot utilizing the numerical model prepared by this method of preparing numerical model can be provided.

Abstract: A method is disclosed for executing a physics simulation in a system comprising a computational platform, a main application stored in the computational platform, a secondary application stored in the computational platform, and a smoothed particle hydrodynamics (SPH) application programming interface (API) implemented in the computational platform. The method defines a SPH call in the SPH API, and by operation of the main application, invokes a software routine using the SPH call. Additionally, by operation of the secondary application, a state of the physics simulation is updated in response to the software routine.

Abstract: A method and its implementation as a program system are disclosed herein for generating an air flow report based upon a cellular model of the air flow of a hard disk drive including a parameterized component approximated by a component parameter list and a map of the partitioned regions of the hard disk drive. The air flow report is generated based upon the cellular model and a partition/region map. The partition/region map lumps the simulation domain into a small number of regions and then calculates the fluxes across boundaries of the regions. These fluxes accumulate the results of many individual cells, averaging out small variations caused by rounding and/or the convergence properties of the specific cellular approach used. A simulation figure of merit is calculated from the air flow report that further refines the accuracy, effectively removing even more noise.

Abstract: Disclosed herein are methods, systems, and devices for determining parameters of an earth formation. Pressure transient data from a formation test can be recorded and analyzed using an analytical model including one or more correction factors derived from an assumption that an induced flow within the formation is hemispherical. Regression analysis of the refined analytical model and the pressure transient data results in accurate earth formation parameters.

Abstract: A method is disclosed for using a three-dimensional seismic image of a subsurface earth volume to construct a geologic model specifying the spatially-varying grain size distribution, porosity, and permeability throughout the volume. The method applies to earth volumes composed of water-lain clastic sedimentary deposits and involves, in one embodiment, (a) identifying the outline forms of geologic bodies in geologic data; (b) using the outline forms of the geologic bodies to determine the spatially-varying grain size distribution within the bodies, guided by assumptions about the nature and behavior of the paleoflow that deposited the bodies; (c) determining rock properties such as, porosity and permeability within the geologic bodies based on grain-size distribution, mineralogy and burial history information.

Abstract: The present invention is a method and system for simulating fluid flow in a cavity, having relevance to the modeling of viscous flows within thin cavities of complex shapes in which heat exchange with the cavity walls may be a governing factor, as in injection molding of plastic parts. Its automated discretization scheme, in which the model cavity is partitioned into macrocells, each macrocell having substantial contact areas with the model cavity walls, eliminates the need for time consuming and expensive manual model configuration required by some modeling methods. Because of the simplicity of the discretization, models based on the invention are less demanding on computer resources then conventional finite element methods, and execute more quickly. A key innovation of the method is a function characterizing the shape of a macrocell that appears in a coefficient in equations governing the flow.

Abstract: A system for simulating a wellbore used for cuttings re-injection, that includes functionality to obtain as input to the system at least one wellbore design parameter for the wellbore, at least one operating parameter for the cuttings re-injection, and a slurry design for a slurry to be injected into the wellbore, functionality to segment the wellbore into a plurality of elements, wherein each element includes a plurality of nodes, and functionality to perform a simulation at a current time interval. The functionality to perform the simulation includes functionality to update a solid accumulation at a bottom of the wellbore at the current time interval and functionality to perform for each of the plurality of nodes, until the wellbore reaches a steady-state condition for the current time interval, at least one calculation using the at least one wellbore design parameter, the at least one operating parameter, and/or the slurry design.

Abstract: The method provides a new fluid simulation technique that significantly reduces the non-physical dissipation of velocity using particles and derivative information. In solving the conventional Navier-Stokes equations, the method replace the advection part with a particle simulation. When swapping between the grid-based and particle-based simulators, the physical quantities such as the level set and velocity must be converted. A novel dissipation-suppressing conversion procedure that utilizes the derivative information stored in the particles as well as in the grid points is developed. Through several experiments, the proposed technique can reproduce the detailed movements of high-Reynolds-number fluids, such as droplets/bubbles, thin water sheets, and whirlpools. The increased accuracy in the advection, which forms the basis of the proposed technique, can also be used to produce better results in larger scale fluid simulations.

Abstract: The method of automatic optimization is applied to a natural gas transport network in the steady state comprising at one and the same time a set of passive works such as pipelines or resistances, and a set of active works comprising regulating valves, isolating valves, compression stations, storage or supply devices, consumption devices, elements for bypassing the compression stations and elements for bypassing the regulating valves, the passive works and the active works being linked together by junctions. The optimization method comprises the determination of values for continuous variables. Intervals of values for the continuous variables and sets of values for the discrete variables are chosen as initial state of the optimization.

Abstract: In some embodiments, a method is provided. A sinusoidal signal is generated that is representative of a wave at an average surface of a liquid. A distance between the average surface of the liquid and a bottom of the liquid is determined. A characteristic of the sinusoidal signal is adjusted as a function of the distance.

Abstract: A method of designing a magnetorheological (MR) fluid energy absorbing damper is provided that uses hydromechanical analysis with lumped parameters to allow a determination as to whether a potential damper design will provide predetermined characteristics, such as a desired dynamic force range and maximum piston velocity, with a selected MR fluid and yield stress and preferably meeting predetermined geometric limitations.

Abstract: Simulating a physical process includes storing, in a computer-accessible memory, state vectors for voxels, where the state vectors correspond to a model and include entries that correspond to particular momentum states of possible momentum states at a voxel. Interaction operations are performed on the state vectors. The interaction operations model interactions between elements of different momentum states according to the model. Move operations performed on the state vectors reflect movement of elements to new voxels according to the model. The model is adapted to simulate a high-Knudsen number flow that has a Knudsen number greater than 0.1.

Abstract: The invention relates to a computer-aided method for predicting particle uptake by a surface of an object moving in a flow current, comprising the following steps: preparing a discrete model of said object and the surrounding flow, including a projection area (21) for projecting said particles formed by a mesh of cells (23) with an identical size, and a mesh of cells (13) of said surface (11); simulating the projection of particles with a different size from said projection area (21) on the surface (11), both the density Di of projected particles and their size distribution being identical for the cell (23); obtaining for each cell (13) of said surface (11) the local uptake efficiency parameter ? by dividing the density Df of particles impacting against said cell (13) by the density Di of particles projected from the projection area (21).

Abstract: A method of measuring flow sections of a turbomachine nozzle sector relative to two reference vanes, wherein: a) a three-dimensional numerical model of the sector is produced by digitizing; b) numerical models are provided for the two reference vanes; c) the reference vanes are reset at the ends of the numerical model of the sector in their assembly relative positions in the nozzle; and d) the flow sections of the sector are determined from sections of the numerical models of the sector and of the reference vanes; the nozzle sector having contact surfaces, and being positioned relative to adjacent nozzle sectors by putting said contact surfaces into abutment, the numerical models of the sector and of the reference vanes including models of the contact surfaces involved in putting them into their assembly relative position, and step c) being performed by putting the contact surfaces into correspondence.

Abstract: A method for simulating digital watercolor image includes: receiving a background texture on a virtual canvas where a watercolor image will be painted; receiving parameters for fluid simulation; and creating brushstrokes by a mouse pointer's movement on the virtual canvas. Further, the method includes converting coordinates of the brushstrokes on the virtual canvas to fit into simulation grids; calculating movement of colors and water through the fluid simulation by using the parameters and the brushstrokes; and simulating the watercolor image based on the brushstrokes by using the calculated movement of the colors and water.

Abstract: A method for automatic feedback control and/or monitoring of an oil and/or gas production system, or part(s) thereof, which includes one or more processing unit(s) and/or one or more flow line(s) and/or one or more controls. The method includes measuring or estimating value(s) of at least one level or pressure or flow rate or load indicator associated with at least one of the processing unit(s) and/or flow line(s) and at least one of calculating setting(s) for the controls based on at least one controlled variable in form of the measured or estimated level or pressure or flow rate or load indicator, including compensating for nonlinearities of the control means or estimating at least one of the oil-, gas-, water-, or liquid flow rate(s) into and/or out of at least one of the processing units and/or flow line(s).

Abstract: Apparatus and method are provided for facilitating simulation of one or more operating characteristics of an electronics rack. The apparatus includes a rack frame, one or more air-moving devices associated with the rack frame, and an adjustable heat source associated with the rack frame. The one or more air-moving devices establish airflow through the rack frame from an air inlet side to an air outlet side thereof, wherein the established airflow through the rack frame is related to airflow through the electronics rack to be simulated. The adjustable heat source heats air passing through the rack frame, with heated air exhausting from the air outlet side of the rack frame simulating heated air exhausting from the electronics rack.

Abstract: A method is provided for facilitating installation of one or more electronics racks within a data center. The method includes: placing a plurality of thermal simulators in the data center in a data center layout to establish a thermally simulated data center, each thermal simulator simulating at least one of airflow intake or heated airflow exhaust of a respective electronics rack of a plurality of electronics racks to be installed in the data center; monitoring temperature within the thermally simulated data center at multiple locations, and verifying the data center layout if measured temperatures are within respective acceptable temperature ranges for the data center when containing the plurality of electronics racks; and establishing the plurality of electronics racks within the data center using the verified data center layout, the establishing including at least one of reconfiguring or replacing each thermal simulator with a respective electronics rack.

Abstract: Multi-Element Structured Arrays (MESAs) present an effective approach to utilizing high pressure drop media that previously were cost prohibitive for use in a traditional filter element. The operational velocity is significantly slowed by the incorporation of numerous elements; thus, eliminating the Darcian flow resistances of the media. The MESAs are readily fitted with two or more replaceable “filter elements” that are able to employ commercially available media, microfibrous entrapped catalyst/sorbent media, unique multilayers and more.

Abstract: Method and apparatus are disclosed for implementing a geometric multi-cell system dynamics model having an embedded grid, the embedded grid having cells with a finer grid size relative to the grid size of other cells. The apparatus includes a digital computer having a computational fluid dynamics model program stored therein, the program having software for iteratively solving transport equations for thermophysical values for the embedded cells and the other cells, and for solving residual equations for the values for each cell using a multi-grid computation method. The program also has other software for computationally manipulating the embedded grid cells to provide composite cells of the same grid size as the other cells and having averaged residual thermophysical values, for allowing the multi-grid computation method to operate on the embedded grid cells.

Abstract: A method and apparatus for analyzing steady state and transient flow in a complex fluid network, modeling phase changes, compressibility, mixture thermodynamics, external body forces such as gravity and centrifugal force and conjugate heat transfer. In some embodiments, a graphical user interface provides for the interactive development of a fluid network simulation having nodes and branches. In some embodiments, mass, energy, and specific conservation equations are solved at the nodes, and momentum conservation equations are solved in the branches. In some embodiments, contained herein are data objects for computing thermodynamic and thermophysical properties for fluids. In some embodiments, the systems of equations describing the fluid network are solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods.

Abstract: Water surface and other effects are efficiently simulated to provide real time or near real time imaging on low-capacity computer graphics computation platforms. Water and other surfaces are modeled using multiple independent layers can be dynamically adjusted in response to real time events. The number of layers used in a given area can be adjusted to reduce computational loading as needed. Different algorithms can be employed on different layers to give different effects. The multiple layer modeling is preferably converted to polygons using an adaptive polygon mesh generator based on camera location and direction in the 3D world to provide automatic level of detailing and generating a minimal number of polygons. The visual effects of water droplets and other coatings on see-through surfaces can be modeled and provided using indirect texturing.

Abstract: A system and method for simulating the flow of a viscoelastic fluid through a channel. The simulation including a interface between a first fluid and a second fluid. The simulation including the formation of a droplet. The simulation includes solving equations governing the viscoelastic flow of the first fluid through the channel, including viscoelastic stress equations that include a normalized relaxation time greater than or equal to 5. The calculations simulate the flow of the first fluid through the channel. The simulation is stable over a period time in which a droplet is formed. The simulation including a level set function that describes the position of the interface between the first and second fluids, and the evolution of the level set function over time describes the shape and position of the interface.

Abstract: Heat trace system design is disclosed. In some embodiments, pipe system data including data associated with pipe system components is received, and in response to receiving a selection of a pipe system component, data associated with a set of one or more pipe system components including the selected pipe system component is automatically extracted from the received pipe system data. One or more heat trace system components for the set of pipe system components are determined based at least in part on the extracted data.

Abstract: A method and system for characterization of fault conditions within a subterranean volume. In one embodiment, the system comprises means for mathematically modeling stress conditions, to predict breakout conditions along a borehole trajectory. The system further comprises means for sensing actual breakout conditions along the borehole . Predictive breakout data is compared with sensed breakout conditions to assess correlation between predictive data and actual data, verifying the accuracy of the stress model. The mathematical model may be revised to reflect the presence of an active fault plane in the volume, the presumed fault plane not being intersected by the borehole. The revised model is used to generate new predictive data. Revising the stress model and assessing correlation between predictive and actual breakout conditions is repeatable to achieve an optimally accurate stress model reflecting fault conditions proximal to but not necessarily penetrated by the borehole.

Abstract: A nonlinear technical product or process described by stochastic system output target values dependent on stochastic system input parameter values, thereby stating discrete technical system dependencies, is optimized by using a Response Surface Methods based on discrete technical system dependencies to generate at least one continuous auxiliary function for the real dependencies of the target values on the input parameter values. Next, an auxiliary function is used to generate at least one optimizing parameter for optimization by an objective function, thereby generating an additional discrete technical system dependence. The technical system is adaptively optimized by repeating the above, using the additional discrete technical system dependence, until the difference of successive optimized optimizing parameters is below a threshold. The final additional discrete technical system dependence is an optimal technical system operating point.

Abstract: A fully-parallelized, highly-efficient compositional implicit hydrocarbon reservoir simulator is provided. The simulator is capable of solving giant reservoir models, of the type frequently encountered in the Middle East and elsewhere in the world, with fast turnaround time. The simulator may be implemented in a variety of computer platforms ranging from shared-memory and distributed-memory supercomputers to commercial and self-made clusters of personal computers. The performance capabilities enable analysis of reservoir models in full detail, using both fine geological characterization and detailed individual definition of the hydrocarbon components present in the reservoir fluids.

Abstract: The present invention generally relates to microfluidics and more particularly to the design of customized microfluidic systems using a microfluidic computer aided design system. In one embodiment of the present invention a microfluidic circuit design method is provided. The method includes developing synthesizable computer code for a design. Next, a microfluidic circuit schematic, including a plurality of symbols for microfluidic components, is generated either interactively or using the synthesizable computer code. The microfluidic circuit schematic is then functionally simulated. The microfluidic components are placed and routed on a template to form a physical layout. Then the physical layout is physically simulated using dynamic simulation models of the microfluidic components; and the physical layout is written to a layout file.

Abstract: Disclosed are techniques for representing and modeling one or more systems in which each system corresponds to an application mode. This may be done for one or more geometries using local and/or non-local couplings. For each application mode, physical quantities are modeled and may be defined using a graphical user interface. Physical properties may be used to model the physical quantities of each system. The physical properties may be defined in terms of numerical values or constants, and mathematical expressions that may include numerical values, space coordinates, time coordinates, and actual physical quantities. Physical quantities and any associated variables may apply to some or all of a geometric domain, and may also be disabled in other parts of a geometrical domain. Partial differential equations describe the physical quantities. One or more application modes may be combined using an automated technique into a combined system of partial differential equations as a multiphysics model.

Abstract: A computer-implemented system for designing a pressure-dosed drain field comprising: a screen that allows a user to specify length and diameter for pipes that comprise the forcemain and discharge assembly, specify the pumping liquid level and discharge elevation, select the number of laterals and enter a specified residual pressure, enter or select a coefficient to be used in calculating friction losses, and enter the diameter and length of each manifold segment; a minor loss screen that calculates minor losses based on fitting selections made by the user; a pump file; a tabulated results screen; a gauge screen; and system and pump curves. A computer-implemented method of using the above system to design a pressure-dosed drain field. The purpose of the system and method is to enable an engineer to design a pressure-dosed drain field that will meet regulatory requirements for the system flow differential.

Abstract: A method for numerical calculation related to fluids which is efficient and less liable to cause numerical diffusion and does not rely upon the functional form of an equation of state and hence can be used as a generalized numerical calculation method. A flux is divided into an advection part, a pressure-dependent part, and a dissipation part, according to physical properties thereof, to calculate a one-dimensional basic equation. A time evolution is executed on the advection part and the pressure-dependent part using an upwind difference type leap-frog difference method and a central difference type leap-frog difference method, respectively.

Abstract: The present invention discloses a method for performing a stochastic analysis of one or more hydrocarbon reservoir exploitation strategies taking into consideration one or more uncertain parameters. The method optimizes an objective function that considers the gain in value of a reservoir management goal attributable to these exploitation strategies. The methodology may be used to decide whether or not to implement a strategy. Alternatively, it may be used to decide which competing strategy will yield the maximum benefit. In another embodiment of the present invention, the value of information obtained from the installation of additional sensors or new measurements is also considered.

Abstract: A computer implemented method and system for the generation of software thermal profiles for applications executed on a set of processors in a simulated environment. Execution of a software program being run on a software simulator is detected and hardware operations for the software program being executed by the set of processors are analyzed to create analyzed information. Then, a thermal index is generated based on the analyzed information.

Abstract: A method for modeling fluid depletion in a reservoir is disclosed. A map is divided into cells. For each of the cells a value is stored that is based at least in part on a physical characteristic of the cell. At least one cell that contains a depletion location is identified along with a depletion amount corresponding to that location. An amount of walkers associated with the depletion location is determined. For each walker, a plurality of steps are calculated with each step to an adjacent cell. Each walker starts in the cell containing the depletion location associated with that walker. The visits of all the walkers are recorded by cell. The fluid depletion of each cell is then assessed based at least in part on the number of walker visits for each cell.

Abstract: A method, practiced by a Hydraulic Fracturing Simulator software, is adapted for simulating a hydraulic fracture in an Earth formation, the formation including a multilayered reservoir, a mesh adapted to overlay the fracture in the formation thereby defining a plurality of elements, the mesh and the fracture defining one or more partially active elements, the method comprising the step of: setting up an Influence Coefficient Matrix, the matrix having a plurality of numbers, the plurality of numbers of the Influence Coefficient Matrix being calculated and determined such that the method, practiced by the Hydraulic Fracturing Simulator software for simulating the hydraulic fracture, models, or take into account an existence of, the multilayered reservoir and the partially active elements.

Abstract: A system and method for drag reduction allows thrust output, fuel efficiency or both to be maximized. Specifically, a rear portion of a body or motor vehicle may be modified to increase thrust output, fuel efficiency or both by creating a stagnation area, a suction inlet and a convex cusp area formed on the rear portion of the motor vehicle. Increasing the concavity or camber or sharpness of the radius of the stagnation area results in greater local pressure coefficient, which results in greater thrust output. The size and shape of the suction inlet and the convex cusp area will also have an effect on thrust output and fuel efficiency.

Abstract: Disclosed are methods for designing vanes that, in use, are expected to be exposed to vibratory loading, in particular rotor or stator vanes for aero engines or turbomachinery. A quantitative characteristic (modeshape) is used which is a measure of correlation between the excitation force to which the vane is to be subjected in use and a vibrational mode of the vane, corresponding to a characteristic frequency of the excitation. The modeshape of the first design is determined. Then the first design is modified to give a second design, by one or more of leaning, sweeping or twisting the blade design, or by altering the local shape of the design, or by altering the material of the design. Then the modeshape of the second design is determined. The design modification gives rise to a reduction in the forced response levels.

Abstract: A simulation system [200] models and optimizes parameters for a pulsed liquid slug boring system employ an energetic fluid [7]. The simulation system [200] employs a fluid flow energy unit [251], an exhaust and retention energy unit [253] and a comminuting energy unit [255] to calculate energies of the system. Total energy unit [257] combines these energies. Fluid flow energy unit 251 receives fluid volume and calculates the fluid flow energy. Exhaust and retention energy unit 253 receives input from the exhaust energy volume unit [243] and mission duration unit [211] to determine the exhaust and retention energy. Comminuting energy unit 255 receives hole depth and hole diameter and specific energy of rock to determine the require comminuting energy. The simulation system [200] operates to determine optimum values of design parameters by searching for the minimum energy solution.

Abstract: An embodiment of the present invention may be a system or method for simulating the flow of a single-phase fluid flow. Markers represent a moving fluid boundary of the single-phase fluid at a first point in time. The moving fluid boundary separates a simulation space into a fluid space and a non-fluid space. The single-phase fluid inhabits the fluid space. A signed distance function is evaluated at points surrounding the moving fluid boundary based upon markers. The curvature of the moving fluid boundary based on the signed distance function is evaluated near the markers in the non-fluid space. The curvature is not evaluated at the moving fluid boundary. The velocity of the fluid is calculated based upon the curvature of the level set in the non-fluid space. Update the position of the moving fluid boundary at a second point in time based on the velocity of the fluid.

Abstract: A system for simulating a wellbore used for cuttings re-injection, that includes functionality to obtain as input to the system at least one wellbore design parameter for the wellbore, at least one operating parameter for the cuttings re-injection, and a slurry design for a slurry to be injected into the wellbore, functionality to segment the wellbore into a plurality of elements, wherein each element includes a plurality of nodes, and functionality to perform a simulation at a current time interval. The functionality to perform the simulation includes functionality to update a solid accumulation at a bottom of the wellbore at the current time interval and functionality to perform for each of the plurality of nodes, until the wellbore reaches a steady-state condition for the current time interval, at least one calculation using the at least one wellbore design parameter, the at least one operating parameter, and/or the slurry design.

Abstract: Two finite difference algorithms on rectangular grids are provided for viscoelastic ink ejection simulations. The first algorithm is first-order accurate in time, while the second algorithm is second-order accurate in time. Each of these two algorithms seamlessly couples or incorporates several features: (1) a coupled level set-projection method for incompressible immiscible two-phase fluid flows; (2) a higher-order Godunov type algorithm for the convection terms in the momentum equations and level set convection equation; (3) the central difference for viscosity, surface tension, and upper-convected derivative terms; (4) an equivalent circuit for inflow pressure (or inflow rate). In the first algorithm, a simple first-order upwind algorithm is designed for the convection term in the viscoelastic stress equations. But in the second algorithm, a second-order Godunov type upwind scheme is employed for the convection term in the viscoelastic stress equations.