Abstract: Determining viable hydraulic fracture scenarios. At least some of the illustrative embodiments include: executing a fracture planning program, and determining a set of schedules from the fracture planning program, each schedule comprising a volume of fracture fluid, amount of proppant, and flow rate of the fracture fluid; providing each schedule of the set of schedules to a stress analysis program, executing the stress analysis program, and determining a set of indications from the stress analysis program, each indication indicative of whether a respect schedule exceeds engineering limits of a tubing string; and providing at least some of the schedules to a fracture simulation program, executing the fracture simulation program, and determining a set of fracture geometries from the fracture simulation program, each fracture geometry corresponding to a respective schedule.

Abstract: A computational “toolbox” can be used for the a priori design of optimized fluidic components. These components include a channel under low-Reynolds number, pressure-driven flow, with an arrangement of grooves cut into the top and/or bottom to generate a tailored cross-channel flow. An advection map for each feature (i.e., groove of a particular shape and orientation) predicts the lateral transport of fluid within the channel due to that feature. Applying the advection maps in sequence generates a representation of the outflow distribution for complex designs that combine one or more features. Therefore, the effect of the complex three-dimensional flow field can be predicted without solving the governing flow equations through the composite geometry, and the resulting distribution of fluids in the channel is used to evaluate how well a component performs a specified task. The toolbox is applied to determine optimal combinations of features for specified mixer sizes and mixing metrics.

Abstract: A method of simulating fluid flow in a three-dimensional volume with a moving boundary is presented. The moving boundary is provided within a computational fluid dynamics (CFD) mesh. The CFD mesh includes mesh points defining three-dimensional mesh cells. Locations of the mesh points at the moving boundary are revised predicated on known motion characteristics thereof such that the CFD mesh undergoes distortion near the moving boundary. The distortion is minimized using an elastic deformation model that does not involve dimensional physical parameters. As a result, a revised CFD mesh is generated. A fluid flow simulation is then computed for the three-dimensional volume using the revised CFD mesh.

Abstract: A simple and efficient method for simulating dispersed bubble flow is provided. Instead of modeling the complex hydrodynamics of numerous small bubbles explicitly, the method approximates the average motion of these bubbles using a continuum multiphase solver. The subgrid interactions among bubbles are computed using a new stochastic solver. Using the proposed scheme, complex scenes with millions of bubbles can be simulated efficiently.

Abstract: One aspect of the present invention provides a method for determining orientation of fibers in a fluid having polymer chains, characterized in that the determining of the orientation of the fibers is performed by taking into consideration an interaction between the fibers and the fluid, wherein the interaction between the fibers and the fluid comprises changes in configuration of the polymer chain to cause the entanglement or adsorption between the fibers and the polymer chains. Another aspect of the present invention provides a method for determining orientation of fibers in a fluid having polymer chains, the fibers in the fluid including a transitional movement and a rotatory movement, the method being characterized in that the determining of the orientation of the fibers is performed by taking into consideration a steric barrier effect on a rotary movement of the fibers.

Abstract: A computerized simulation method for evaluating a tire performance is disclosed. Firstly, a simulation of a rolling tire is made by contacting a tire model with a road surface model, and chronological data about node points of the tire model which node points appear in the outer surface of the tire model during rotating are stored. Then, a domain of a fluid is defined, and a simulation of the fluid is made by the use of the chronological data, wherein a small gap is formed between the tire model and the road surface model, and the fluid domain is also defined in the gap.

Abstract: A heat transfer apparatus and related methods are provided. The heat transfer apparatus and related methods more evenly distribute fluid flow in two-phase heat exchange systems by restricting fluid flow to one or more tube.

Abstract: A computerized method of identifying a deployment of a plurality of hydraulic sensors in a liquid distribution network, optionally according to leakage detection performance.

Abstract: A method and system of determining a volume of a vehicle component. The system includes a computing device having a host application and analytical tool. The host application is configured to present an interface to a user and allow the user to modify a component model. The component model includes a three-dimensional representation of at least a portion of a vehicle component and has a virtual housing at least partially defining an internal cavity. The analytical tool is configured to generate a solid body representing the internal cavity and calculate a volume of the solid body. The analytical tool is configured to determine whether the internal cavity is bounded by the virtual housing, counting a number of external and internal faces traversed by the fluid during the simulated fluid flow, and determining a shortest path between the selected external and internal faces.

Abstract: A method for modeling deformation in subsurface strata, including defining physical boundaries for a geomechanical system. The method also includes acquiring one or more mechanical properties of the subsurface strata within the physical boundaries, and acquiring one or more thermal properties of the subsurface strata within the physical boundaries. The method also includes creating a computer-implemented finite element analysis program representing the geomechanical system and defining a plurality of nodes representing points in space, with each node being populated with at least one of each of the mechanical properties and the thermal properties. The program solves for in situ stress at selected nodes within the mesh.

Abstract: A method for predicting dynamic behavior of an aircraft structure, the aircraft including at least one rotary device including at least one rotor guided in rotation by at least one fixed-bracket bearing containing a shock-absorbing fluid film. The method includes using a global digital model of the aircraft structure including a basic digital model of each rotary device, generating a non-linear digital model of each shock-absorbing film (for example by developing a cavitation-free model and a model with cavitation), integrating the model of each shock-absorbing film with the global model, applying interference to at least one rotor of a rotary device, and calculating the vibration frequency of at least one rotor of each rotary device and the corresponding frequency of the vibrations generated in a critical portion of the aircraft structure.

Abstract: Exemplary embodiments of the present invention relate to a numerical analysis data providing system. The system can include a server side computer which may provide numerical analysis data to a user side computer connected to a network. The server side computer can record a material name and a property item with corresponding to any one or more types of property values from a mechanical property value, a thermal physical property value, and an electromagnetic property value as for plural materials. Such server side computer can extract one or more types of material property data from the mechanical property value, the thermal physical property value, and the electromagnetic property value corresponding to the material name and the property item recorded by a material property data storage arrangement based on the material name and the property item transmitted from the user side computer. In addition, the server side computer can transmit the extracted data to the user side computer.

Abstract: The present method aims at preventing breakdown by selecting conditions for preventing the occurrence of voids via analysis. According to the present method, a database related to the time variation of amount of generated gas or number of moles when a solid member is heated is experimentally formulated in advance, and the amount of generated gas from the member and the specific heat ratio of gases are entered for analysis in order to predict via analysis the generation of voids within the resin material when the solid member is heated, to thereby select conditions for preventing voids.

Abstract: Methods and systems are provided for generating a development plan for a hydrocarbon asset. A high-fidelity computer model of a hydrocarbon asset is created. A low-fidelity computer model of the hydrocarbon asset is created. The low-fidelity computer model is iterated on to an interim solution. A comparison is generated of the interim solution to a solution obtained from a simulation of the high-fidelity computer model at the variables of the interim solution. The low-fidelity computer model is calibrated based, at least in part, on the comparison. The development plan for the hydrocarbon asset is generated based, at least in part, on a result from the calibrated low-fidelity computer model. The low-fidelity computer model is a mixed-integer nonlinear programming problem with complementarity.

Abstract: Fluid-flow simulation over a computer-generated surface is generated using inviscid and viscous simulations. A fluid-flow mesh of fluid cells is obtained. At least one inviscid fluid property for the fluid cells is determined using an inviscid fluid simulation that does not simulate fluid viscous effects. A set of intersecting fluid cells that intersects the surface are identified. A surface mesh polygon of the surface mesh is identified for each intersecting fluid cell. At least one boundary-layer fluid property for each identified surface mesh polygon is determined using the at least one inviscid fluid property of the corresponding intersecting fluid cell and a boundary-layer simulation that simulates fluid viscous effects.

Abstract: The present invention incorporates a device for measuring fluid movement conditions at a distance, such as LIDAR, into an underlying device for extracting useful work from a moving fluid, such as a wind turbine, for the purpose reducing the cost of designing and manufacturing the underlying device.

Abstract: A prediction of whether a point on a computer-generated surface is adjacent to laminar or turbulent flow is made using a transition prediction technique. A plurality of boundary-layer properties at the point are obtained from a steady-state solution of a fluid flow in a region adjacent to the point. A plurality of instability modes are obtained, each defined by one or more mode parameters. A vector of regressor weights is obtained for the known instability growth rates in a training dataset. For each instability mode in the plurality of instability modes, a covariance vector is determined, which is the covariance of a predicted local growth rate with the known instability growth rates. Each covariance vector is used with the vector of regressor weights to determine a predicted local growth rate at the point. Based on the predicted local growth rates, an n-factor envelope at the point is determined.

Abstract: A computer-implemented method of determining the pitch stability of an airfoil system, comprising using a computer to numerically integrate a differential equation of motion that includes terms describing PID controller action. In one model, the differential equation characterizes the time-dependent response of the airfoil's pitch angle, ?. The computer model calculates limit-cycles of the model, which represent the stability boundaries of the airfoil system. Once the stability boundary is known, feedback control can be implemented, by using, for example, a PID controller to control a feedback actuator. The method allows the PID controller gain constants, KI, Kp, and Kd, to be optimized. This permits operation closer to the stability boundaries, while preventing the physical apparatus from unintentionally crossing the stability boundaries. Operating closer to the stability boundaries permits greater power efficiencies to be extracted from the airfoil system.

Abstract: The disclosed methods, systems, and software are described to optimize fracture characteristics and simulate fluid flow rates in a well model. The well model, which includes at least one fracture intersecting a production well, is generated with static and dynamic data. Fluid flow in the well model is simulated to obtain simulated fluid flow rates between fractures and the well. Fracture properties, such as length, height and aperture, are then updated responsive to measured and simulated fluid flow rates. Multiple simulation runs and updating of the fracture properties can be performed until the simulated fluid flow rates converge with the measured fluid flow rates. Pressure data can be used to determine gridblock permeability, which in turn helps constrain the model, thus providing more reliable fracture properties. Uncertainty ranges of the fracture properties can also be calculated.

Abstract: Methods and systems for optimizing the design of aerodynamic surfaces are provided. They include in the design process of an aerodynamic surface of a body which moves through a flow field, the following: preparation of a CFD simulation of said body; resolution of the CFD calculation for specific conditions of the flow field and for obtaining flow lines of the flow on said surface; calculation of at least one geometrical variable of the surface, such as, for example, the normal curvature, according to the direction of said flow lines; and simultaneous visualization of the surface of said body and of at least one geometric variable. Additionally, at least one variable of the aerodynamic flow, such as, for example, a pressure gradient, may be calculated and visualized.

Abstract: Methods, systems, and computer readable media are provided for fast updating of oil and gas field production optimization using physical and proxy simulators. A base model of a reservoir, well, or a pipeline network is established in one or more physical simulators. A decision management system is used to define uncertain parameters for matching with observed data. A proxy model is used to fit the uncertain parameters to outputs of the physical simulators, determine sensitivities of the uncertain parameters, and compute correlations between the uncertain parameters and output data from the physical simulators. Parameters for which the sensitivities are below a threshold are eliminated. The decision management system validates parameters which are output from the proxy model in the simulators. The validated parameters are used to make production decisions.

Abstract: A method of modeling a subsurface hydrogeologic system within a study area can include applying an electric current directly to the subsurface hydrogeologic system and monitoring a magnetic field generated by the electric current in the subsurface hydrogeologic system. An electric current distribution can be identified within the study area by an inversion process based on the magnetic field. A model of the electric current distribution can be created.

Abstract: An improved method and computer system is provided, for analyzing the two- or three-dimensional movement of a fluid in an absorbent article that comprises fluid-swellable composite material, which comprises a fluid-swellable solid material, and that comprises void spaces in the fluid-swellable composite material, the fluid-swellable composite material being defined by a virtual two-dimensional or three-dimensional mesh, to determine the performance or property of at least one feature of the fluid-swellable composite material thereof, or the absorbent article, the method comprising: inputting one or more properties for the fluid-swellable (composite) material into a model and evaluating this to determine the performance of a feature of the absorbent article or fluid-swellable composite material thereof, the virtual model.

Abstract: Apparatus and method are provided for facilitating simulation of heated airflow exhaust of an electronics subsystem, electronics rack or row of electronics racks. The apparatus includes a thermal simulator, which includes an air-moving device and a fluid-to-air heat exchanger. The air-moving device establishes airflow from an air inlet to air outlet side of the thermal simulator tailored to correlate to heated airflow exhaust of the electronics subsystem, rack or row of racks being simulated. The fluid-to-air heat exchanger heats airflow through the thermal simulator, with temperature of airflow exhausting from the simulator being tailored to correlate to temperature of the heated airflow exhaust of the electronics subsystem, rack or row of racks being simulated. The apparatus further includes a fluid distribution apparatus, which includes a fluid distribution unit disposed separate from the fluid simulator and providing hot fluid to the fluid-to-air heat exchanger of the thermal simulator.

Abstract: A method is provided for simulating a physical process such as fluid flow in porous media by performing a fine-grid calculation of the process in a medium and re-using the fine grid solution in subsequent coarse-grid calculations. For fluid flow in subsurface formations, the method may be applied to optimize upscaled calculation grids formed from geologic models. The method decreases the cost of optimizing a grid to simulate a physical process that is mathematically described by the diffusion equation.

Abstract: Method for reducing a 3D joint inversion of at least two different types of geophysical data acquired by 3-D surveys (21) to an equivalent set of ID inversions. First, a 3D inversion is performed on each data type separately to the yield a 3-D model of a physical property corresponding to the data type (22). Next, a ID model of the physical property is extracted at selected (x,y) locations. A ID simulator (23) and the ID model of the physical property is then used at each of the selected locations to create a synthetic ID data set at each location (24). Finally, the ID synthetic data sets for each different type of geophysical data are jointly inverted at each of the selected locations, yielding improved values of the physical properties. Because the joint inversion is a ID inversion, the method is computationally advantageous, while recognizing the impact of 3-D effects.

Abstract: A system and method for accurately measuring supply gas consumed by a particular process control component within a process control system is disclosed. Enhanced measurement accuracy is derived from measuring the consumption of the process control component in a normal operating mode of the process control system. The amount of fluid expended by one process control component is separated by a fluid control system from the amount of supply gas expended in actuating other process control components. The amount of fluid expended by each component may be determined by measuring a decrease in a fluid within a vessel having a known quantity that independently supplies supply gas to each component during its operation.

Abstract: A method of improving a 2D or 3D multiblock-structured mesh of an object that moves through a fluid medium, such as an aircraft. The method is used in the design of the object in connection with a defined analytical scheme, the mesh having been generated with rectangular (2D) or hexahedral (3D) blocks under constraints that may cause the mesh to include irregular blocks comprising the following steps: finding irregular blocks in the mesh; degenerating the irregular blocks merging contiguous edges or sides in one edge or side so that triangular (2D) or pentahedral (3D) blocks are formed and extending the modifications all along the blocks that are propagated; and generating several rectangular (2D) or hexahedral (3D) regular blocks inside the triangular (2D) or pentahedral (3D) blocks. A system for carrying out the method.

Abstract: Techniques and apparatus inhibit, limit, or remove biofouling and certain inorganic accumulations, to increase the longevity of accurate in-situ oceanographic and other underwater measurements and transducing processes. The invention deters formation of an initial bacterial layer and other precipitation, without harming the environment. The invention integrates an ultrasonic source into a sensor or other device, or its supporting structures. The ultrasonic source vibrates one or more critical surfaces of the device at a frequency and amplitude that dislodge early accumulations, thus preventing the rest of the fouling sequence. The ultrasonic driver is activated for short periods and low duty cycles, and in some cases preferably while the device is not operating.

Abstract: A computer-aided method for predicting particle uptake by a surface of an object moving in a flow current. The method includes preparing a discrete model of the object and the surrounding flow, including a projection area for projecting the particles formed by a mesh of cells with an identical size, and a mesh of cells of the surface. The method also includes simulating the projection of particles with a different size from the projection area on the surface, both the density Di of projected particles and their size distribution being identical for the cell. The method further includes obtaining for each cell of the surface the local uptake efficiency parameter by dividing the density Df of particles impacting against the cell by the density Di of particles projected from the projection area.

Abstract: Methods of using computer based models for simulating the physical behavior of bodily fluids with absorbent articles.

Abstract: A computer-implemented method for simulating fluid flow using a lattice Boltzmann (LB) approach and for solving scalar transport equations is described herein. In addition to the lattice Boltzmann functions for fluid flow, a second set of distribution functions is introduced for transport scalars.

Abstract: Methods, systems, and apparatus, including computer programs encoded on a computer storage medium, for interactively simulating an injection mold model. A three-dimensional CAD model is identified representing an injection mold cavity. The mold cavity includes a location of at least one gate. A potential fill pattern is determined for injection of a material into the injection mold cavity. The determined fill pattern is based at least in part on the geometry and dimensions of the modeled mold cavity and the location of the at least one gate. A strip model of the CAD model is generated based at least in part on the determined fill pattern. The strip model is used to perform a strip analysis simulating injection of the material within the injection mold cavity.

Abstract: A method and system for assessment of virtual stent implantation in an aortic aneurysm is disclosed. A patient-specific 4D anatomical model of the aorta is generated from the 4D medical imaging data. A model representing mechanical properties of the aorta wall is adjusted to reflect changes due to aneurysm growth at a plurality of time stages. A stable deformation configuration of the aorta is generated for each time stages by performing fluid structure interaction (FSI) simulations using the patient-specific 4D anatomical model at each time stage based on the adjusted model representing the mechanical properties of the aorta wall at each time stage. Virtual stent implantation is performed for each stable deformation configuration of the aorta and FSI simulations are performed for each virtual stent implantation.

Abstract: A client-server based system for building and executing flows (i.e., interconnected systems of algorithms). The client allows a user to build a flow specification and send the flow specification to the server. The server assembles the flow from the flow spec and executes the flow. A decision flow builder allows the user to build a flow targeted for the analysis/optimization of decisions (modeled by decision variables) regarding a plurality of assets in view of various underlying uncertainties (modeled by uncertainty variables). The user may specify a global objective (as a function of asset level statistics) as well and one or more constraints for the optimization. The flow may account for inter-asset correlations and inter-asset dependencies between uncertainty variables, and, inter-asset constraints between decision variables.

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: Fluid-flow simulation over a computer-generated aircraft surface is generated using inviscid and viscous simulations. A fluid-flow mesh of fluid cells is obtained. At least one inviscid fluid property for the fluid cells is determined using an inviscid fluid simulation that does not simulate fluid viscous effects. A set of intersecting fluid cells that intersects the aircraft surface are identified. One surface mesh polygon of the surface mesh is identified for each intersecting fluid cell. A boundary-layer prediction point for each identified surface mesh polygon is determined. At least one boundary-layer fluid property for each boundary-layer prediction point is determined using the at least one inviscid fluid property of the corresponding intersecting fluid cell and a boundary-layer simulation that simulates fluid viscous effects. At least one updated fluid property for at least one fluid cell is determined using the at least one boundary-layer fluid property and the inviscid fluid simulation.

Abstract: A method and system for multi-scale anatomical and functional modeling of coronary circulation is disclosed. A patient-specific anatomical model of coronary arteries and the heart is generated from medical image data of a patient. A multi-scale functional model of coronary circulation is generated based on the patient-specific anatomical model. Blood flow is simulated in at least one stenosis region of at least one coronary artery using the multi-scale function model of coronary circulation. Hemodynamic quantities, such as fractional flow reserve (FFR), are computed to determine a functional assessment of the stenosis, and virtual intervention simulations are performed using the multi-scale function model of coronary circulation for decision support and intervention planning.

Abstract: Methods and apparatus for simulating fluid motion using procedural shape growth. In a vector-based, fluid motion simulation technique, fluid location may be defined by groups of one or more polygons deposited on a digital canvas. Two or more polygons may overlap. The polygons may be semitransparent. To simulate fluid motion, vertices that specify the edges of each deposited polygon are independently moved. By moving the vertices, a polygon may grow, and fluid motion effects may be simulated, including but not limited to directional flow and blending effects. A randomization technique may be applied to the movement at each vertex to simulate the non-uniform spreading of fluids. Overlapped polygons may be blended with overlapping polygons to simulate the mixing of fluids. The technique may be applied, for example, in watercolor painting simulation, where groups of one or more polygons are deposited using brush strokes.

Abstract: A simulation system (101) includes an X-ray CT device (11) that obtains a tomographic image of a porous sample, and a simulation device (14) that simulates a mercury intrusion method by processing a laminated tomographic image of the sample. The simulation device (14) includes a modeling means which processes the laminated tomographic image of the sample, and which models an internal structure of the sample, a minimum-diameter obtaining means that obtains a minimum entrance diameter when mercury enters in a pore of the sample at a predetermined pressure based on a surface energy of the sample and a pressure, and a means that simulates a liquid entering in the interior of the pore from one surface of the sample based on a diameter of the pore of the modeled sample and the minimum entrance diameter.

Abstract: A system is disclosed including a three-dimensional object having a non-conforming region, and a photogrammetry device adapted to scan the three-dimensional object. The system further includes optical reference targets and a controller structured to perform functions of repairing the three-dimensional object. The controller commands the photogrammetry device to scan the three-dimensional object, and calculates a nominal surface location and contour for the three-dimensional object. The controller further commands the photogrammetry device to scan the non-conforming region of the three-dimensional object, and calculates a material removal tool path comprising a path adapted to remove material from the object located beyond the nominal surface location and contour. The controller generates a solid model of the damaged region of the object based on the nominal surface location and contour, and computes a material addition tool path according to the solid model.

Abstract: A computer-implemented method for simulating flow and acoustic interaction of a fluid with a porous medium includes simulating activity of a fluid in a first volume adjoining a second volume occupied by a porous medium, the activity of the fluid in the first volume being simulated so as to model movement of elements within the first volume and using a first model having a first set of parameters, simulating activity of the fluid in the second volume occupied by the porous medium, the activity in the second volume being simulated so as to model movement of elements within the second volume and using a second model having a second set of parameters and differing from the first model in a way that accounts for flow and acoustic properties of the porous medium, and simulating movement of elements between the first volume and the second volume at an interface between the first volume and the second volume.

Abstract: Methods, computer-readable mediums, and systems analyze hydrocarbon production data from a subsurface region to determine geologic time scale reservoir connectivity and production time scale reservoir connectivity for the subsurface region. Compartments, fluid properties, and fluid distribution are interpreted to determine geologic time scale reservoir connectivity and production time scale reservoir connectivity for the subsurface region. A reservoir connectivity model based on the geologic time scale and production time scale reservoir connectivity for the subsurface region is constructed, wherein the reservoir connectivity model includes a plurality of production scenarios each including reservoir compartments, connections, and connection properties for each scenario. Each of the production scenarios is tested and refined based on production data for the subsurface region.

Abstract: A finite difference level set projection algorithm on multi-staged quadrilateral grids for simulation of a split liquid film between two rollers. From one stage to the next, the number of meshes in the vertical direction changes by a select factor to provide a good balance between resolution and number of meshes. Having fewer meshes at the nip (the smallest gap between the rollers), the multi-staged quadrilateral grid allows a much bigger time step and requires much less CPU time to reach the steady state.

Abstract: A method is disclosed having application notably towards ranking earth models responsive to dynamic heterogeneity. A plurality of earth models representing a subsurface reservoir are provided. Streamline analysis for each of the plurality of earth models is conducted. Flow Capacity (F) vs. Storage Capacity (?) curves are constructed for each of the plurality of earth models based on the streamline analysis. Dynamic heterogeneity for each of the plurality of earth models is computed from the Flow Capacity (F) vs. Storage Capacity (?) curves constructed for each of the plurality of earth models. The plurality of earth models are ranked responsive to dynamic heterogeneity.

Abstract: A method is disclosed including simulating a hydraulic fracture in an Earth formation, the formation including a multilayered reservoir, wherein a mesh overlays the hydraulic fracture thereby defining a plurality of fracture elements. The method further includes calculating and determining an influence coefficient matrix having spatially related indices, wherein the spatially related indices relate influence coefficient matrix elements to corresponding fracture elements. The mesh overlays the hydraulic fracture in more than one layer of the multilayered reservoir. In one embodiment, the method is executed on a program storage device readable by a machine tangibly embodying a program of instructions executable to perform the method.

Abstract: A method for constructing a velocity vector field from a grid and a set of fluxes for each face of the grid cells. The cells are first subdivided and internal fluxes are calculated for each cell subject to the constraints of the flux for each cell and to achieve the minimum energy state for the each cell. The minimum energy state is computed efficiently using a divergence-free correction method without introducing a pressure variable. Then, the velocity vector field is constructed from the subcell fluxes using mixed finite element interpolation.

Abstract: A method of predicting earth stresses in response to changes in a hydrocarbon-bearing reservoir within a geomechanical system includes establishing physical boundaries for the geomechanical system, acquiring logging data from wells drilled, and acquiring seismic data for one or more rock layers. The well and seismic data are automatically converted into a three-dimensional digital representation of one or more rock layers within the geomechanical system, thereby creating data points defining a three-dimensional geological structure. The method also includes (a) applying the data points from the geological structure to derive a finite element-based geomechanical model, and (b) initializing a geostatic condition in the geomechanical model, and then running a geomechanics simulation in order to determine changes in earth stresses associated with changes in pore pressure or other reservoir characteristics within the one or more rock layers.

Abstract: Method and systems for modeling fractures in quasi-brittle materials are provided. An exemplary method included generating a model that incorporates a unified creep-plasticity (UCP) representation into a constitutive model for a ductile rock. The model may be used in a finite element analysis to model hydraulic fractures in the ductile rock.

Abstract: A computer system for modelling and controlling a hydrocarbon reservoir through management of fluid flow at individual wells. The computer system has program instructions which operate a computer model which uses oilfield production data to provide a model of future production. The model comprises an optimal regression model which represents injector and producer wells whose fluid flow characteristics are highly correlated with the fluid flow characteristics of the well of interest; the application of parsimonious information criterion techniques to identify well pairs that statistically contribute information to the optimal regression model; and a statistical reservoir model comprising the product of the optimal regression model and a significance matrix. The system is also provided with control means, responsive to the output of the computer model in order to control wells in the hydrocarbon reservoir.