Crack surfaces are assumed to stick or slip according to a linear friction law. For an array of echelon cracks the ratio of maximum slip to array length significantly underestimates the difference between the driving stress and frictional resistance. The ratio of maximum slip to crack length slightly overestimates this difference. HampsonRussell offers a fully integrated suite of world-class geophysical interpretation tools for reservoir characterization.Known for its ease of use, HampsonRussell makes sophisticated geophysical techniques accessible to all levels of geophysicists.New features for reservoir characterization. New features in HampsonRussell version 10.4.1 delivers continual improvements in reservoir.
IcyDwarf calculates the coupled physical-chemical evolution of an icy dwarf planet or moon. As of version 16.3, the code calculates:
- The thermal evolution of an icy planetary body (moon or dwarf planet), with no chemistry, but with rock hydration, dehydration, hydrothermal circulation, core cracking, tidal heating, and porosity. The depth of cracking and a bulk water:rock ratio by mass in the rocky core are also computed.
- Whether cryovolcanism is possible by the exsolution of volatiles from cryolavas.
- Equilibrium fluid and rock chemistries resulting from water-rock interaction in subsurface oceans in contact with a rocky core, up to 200ºC and 1000 bar.
IcyDwarfPlot creates interactive displays of outputs from the following IcyDwarf functionalities:
- Thermal Evolution
- Core cracking
- Equilibrium fluid and rock compositions.
There is currently no display of cryovolcanism or geochemistry outputs from IcyDwarf.
The two codes can run independently of each other, so it's not necessary to install both if you only need one.
The installation steps outlined below are valid for Mac OS 10.9+. IcyDwarf and IcyDwarfPlot could also run on Windows and Linux, but compilation instructions are not set up and external I/O handling needs to be modified in the source code.
Install R
R is needed only for IcyDwarf, to run the geochemistry package CHNOSZ.Go to http://www.r-project.org and follow instructions.
Install CHNOSZ
CHNOSZ is needed only for IcyDwarf. Open R using either the installed application icon or in a terminal by typing
In R, type the command
Install Rcpp and RInside
Rcpp and RInside are libraries that allow R applications to be embedded in C or C++ codes. They are needed only for IcyDwarf. Go to http://cran.r-project.org/web/packages/Rcpp/index.html and http://cran.r-project.org/web/packages/RInside/index.html, or directly to http://dirk.eddelbuettel.com/code/rcpp to download the respective archives. On Mac, unzip the archives in /Library/Frameworks/R.framework/Resources/library/, so that Rcpp and RInside are two subfolders of library.
Install IPHREEQC
The IPHREEQC library is a module that allows the PHREEQC application to be embedded in C or C++ codes. It is needed only for IcyDwarf. Go to http://wwwbrr.cr.usgs.gov/projects/GWC_coupled/phreeqc to download IPHREEQC and follow the default installation instructions:
Install gcc
In Mac OS 10.8+, the default compiler clang has replaced the compiler gcc, which is needed to take advantage of the parallel computing capabilities of the PlanetSystem and WaterRock_ParamExplor routines of IcyDwarf. Go to http://hpc.sourceforge.net and follow the instructions there to download and install gcc.Once installed, you might need to break the symbolic link between the command gcc and clang by typing:
Install SDL2 (IcyDwarfPlot only)
SDL2 is a graphic library. Go to http://www.libsdl.org/projects. Download and install SDL2, SDL2_image, and SDL2_ttf. SDL2_mixer is not needed as the code doesn't play music for you yet.
Install IcyDwarf
Go to https://github.com/MarcNeveu/IcyDwarf. Click the green Clone or download button to the right of the page, then either:
- Download ZIP on the bottom right. Unzip IcyDwarf-master.zip. Rename the unzipped folder IcyDwarf-master to IcyDwarf.Move the renamed IcyDwarf folder to any folder you would like, we will call it Path_to_GitFolder here.
- if you are familiar with GitHub, you can clone the directory with your favorite tool (I use Git within the Eclipse developing environment).
All source files should be in:
- /Path_to_GitFolder/IcyDwarf/IcyDwarf and subfolders
- /Path_to_GitFolder/IcyDwarf/IcyDwarfPlot and subfolders.
Start the code
The executable files are:
- /Path_to_GitFolder/IcyDwarf/IcyDwarf/Release/IcyDwarf (no extension)
- /Path_to_GitFolder/IcyDwarf/IcyDwarfPlot/Release/IcyDwarfPlot (no extension)
Input files
The respective input files are located in the IcyDwarf/Input and IcyDwarfPlot/Input folders.
Output files
Thermal (± orbital) evolution code
For each file name, the initial character x is 0 for the first/only object and incremented by 1 for each additional object. Thermal and crack output files can be read and displayed by IcyDwarfPlot.
- xCrack_stresses.txt: Internal stresses accounted for by the core cracking subroutine (Neveu et al. 2015, JGR). There are n_zones rows (one per grid zone from the center to the surface) printed at each time interval. Columns list, respectively:
- grid zone radius (in km)
- pressure (in MPa)
- brittle strength (in MPa)
- critical stress intensity (in MPa m^0.5)
- stress intensity from thermal expansion mismatch at grain boundaries (in MPa m^0.5)
- pore fluid pressure (in MPa)
- net pressure (stress) resulting from rock hydration (in MPa)
- old crack size prior to hydration/dehydration (in m)
- old crack size prior to mineral dissolution/precipitation (in m)
- current crack size (in m).Outputs are zero outside of the core.
- xCrack_depth_WR.txt: The bulk water:rock mass ratio in the fractured zone. This file has three columns:
- time (in Gyr)
- depth below seafloor of the fractured zone (in km)
- water:rock ratio by mass in cracked zone.Outputs are zero if the core is not cracked or if there is no liquid.
- xHeats.txt: Cumulative heats (in erg) produced or consumed by endogenic and exogenic processes. The six columns describe:
- time (in Gyr)
- radiogenic heat
- gravitational heat
- heat of rock hydration
- heat consumed in rock dehydration
- heat from tidal dissipation.
- xOrbit.txt (only for simulations with a nonzero host planet mass and in which the moon's orbit is allowed to change): Orbital parameters. Columns list:
- time (in Gyr)
- semi-major axis (in km)
- osculating semi-major axis in km (0 if no resonance)
- eccentricity
- product of eccentricity and cosine of resonant angle
- product of eccentricity and sine of resonant angle
- resonant angle (in degrees)
- total tidal dissipation (in W)
- equivalent k2/Q for the moon (Segatz et al. 1988).
- xThermal.txt: There are n_zones rows for each grid zone, repeated total time/timestep times, i.e. for each time interval. Columns list, respectively, in each grid zone:
- grid zone radius (in km),
- grid zone temperature (in K)
- mass of rock (in g)
- mass of water ice (in g)
- mass of ammonia dihydrate (in g)
- mass of liquid water (in g)
- mass of liquid ammonia (in g)
- Nusselt number (if >1, convection)
- fraction of amorphous ice (always zero, a legacy of Desch et al. 2009)
- thermal conductivity (in W m^-1 K^-1)
- degree of hydration (0: fully dry; 1: fully hydrated)
- porosity
- integer indicating whether the grid zone is fractured, and by which process (Neveu et al. 2015, JGR): 0 = no cracks; 1 = cracks from thermal contraction; 2 = cracks from thermal expansion; 3 = cracks from hydration; 4 = cracks from dehydration; 5 = cracks from pore water dilation; 6 = mineral dissolution widening; 7 = mineral precipitation shrinking; -1 = mineral precipitation clogging; -2: clogging from hydration swelling
- tidal heating rate (in W).
In addition, each simulation with a nonzero host planet mass produces following files. Each of the last three files is read in N_moon x N_moon matrices, where N_moon is the number of moons. Matrices are symmetric since they describe interactions between pairs of moons. Element (x, y) represents interactions between the xth and yth worlds as specified in IcyDwarfInput. The first matrix is output at the first time step. Subsequent matrices are output following a time stamp that corresponds to the time at which pairs of moons get in and out of resonance.
- Primary.txt: Over time in Gyr (first column), the Q of the primary (second column) and the mass of any ring in kg (third column).
- Resonances.txt (for moon system): Values are integers j if the mean motions of the corresponding moons are commensurate in j+1:j ratios with j≤5, and if the migration of the moons is convergent (j dn_inner moon/dt ≤ (j+1) dn_outer moon/dt since dn/dt < 0 for expanding orbits). Values are 0 otherwise. If a moon is in resonance with only one other moon, the code computes moon-moon interactions (value in ResAcctFor below = j), otherwise interactions may be ignored (value in ResAcctFor = 0).
- ResAcctFor.txt: Stands for 'Resonances Accounted For'. A nonzero value in Resonance above is accounted for if a moon is in resonance with only one other moon. Otherwise, the code cannot compute the orbital evolution resulting from the interactions between more than two moons. In that case, the resonance accounted for is that between the pair of moons for which j is smallest (resonance for which the most moon-moon conjunctions occur per orbit). For equal values of j (e.g. for a 4:2:1 resonance, j would be 1 between the inner and middle moon, and also 1 between the middle and outer moon), the newer resonance is ignored. For moons with nonzero values, orbital evolution is computed by an averaged Hamiltonian subroutine (Meyer & Wisdom 2008). Otherwise, orbital evolution is computed solely due to effects from moon-primary and moon-ring interactions, ignoring moon-moon interactions.
- PCapture.txt: This output is not taken into account in computations, but provides an indicative probability of capture into resonance based on the equations of Borderies & Goldreich (1984). Whether or not capture occurs in a simulation depends on the outcome of orbital evolution computed with the averaged Hamiltonian routine. This matrix is not made symmetric, so usually the value of a coefficient in a position symmetric to that of a nonzero value is 0. In that case, only the nonzero value is meaningful.
Cryolava code
The cryolava routine outputs three files:
- Cryolava_molalities.txt (10 columns, n_ice_or_crust_grid_zones rows) shows the cryolava content in H2, CH4, CH3OH, CO, CO2, NH3, N2, H2S, SO2, Ar in mol per kg of liquid water
- Cryolava_partialP.txt, with the same layout as the molalities file, shows the partial pressure of each of these 10 species
- Cryolava_xvap.txt has the same amount of rows, but only 6 columns which show the depth under the surface (km), total gas pressure (bar), volumic vapor fraction x_vap (a dimensionless indicator of exsolution), fluid cryolava density (kg m-3), stress intensity K_I at the crack tip (Pa m^0.5), a boolean (0: no crack propagation; 1: crack propagation).
Compression code
The compression routine outputs one file, Compression.txt, which provides pressures and densities as a function of radius, both accounting for self-compression (output) and not accounting for it (output of the thermal code). The file structure, format, and units are explained in the file itself.
WaterRock_ParamExplor code
This routine outputs a file, ParamExploration.txt, that looks much like the PHREEQC selected output specified in the IcyDwarf/PHREEQC-3.1.2/io folder, with a few added columns at the beginning (starting T in celsius, P in bar, pH, pe, log fO2 at FMQ(T,P) buffer, pe-FMQ). The file is formatted for easy import into a spreadsheet, with each line describing a different simulation. Lines filled with zeros are PHREEQC simulations that did not converge.
The PHREEQC input file, IcyDwarf/PHREEQC-3.1.2/io/inputIcyDwarf, can be modified, but be aware that IcyDwarfPlot will plot results accurately only if the SELECTED_OUTPUT block is left unchanged.
If you wish to modify the code, set up your compiler and linker so that all the relevant flags are added.
Compiler setup
My compiling instructions look like this:
For IcyDwarf (gcc 8.3 on Mac OS 10.14 Mojave):
For IcyDwarfPlot (gcc 6.2 on Mac OS 10.12 Sierra):
You might need to specify the full path to gcc (e.g. /usr/local/bin/gcc) rather than simply the gcc alias.
Your include directories might be more simply found at -I/usr/include.
Email me if you have any issues.
If you communicate or publish scientific results using this code, please acknowledge one of the references listed below from newest to oldest. Each describes the development of one piece of the code. Thanks!
Neveu M., Rhoden A. (2019) Evolution of Saturn’s mid-sized moons. https://doi.org/10.1038/s41550-019-0726-y. (Fully coupled thermal-orbital evolution with moon-primary, moon-ring, and simplified moon-moon interactions, ability to simulate several objects simultaneously)
Neveu M., Desch S., Castillo-Rogez J. (2017) Aqueous geochemistry in icy worldinteriors: equilibrium fluid, rock, and gas compositions, and fate of antifreezes andradionuclides. Geochimica & Cosmochimica Acta 212, 324-371. https://doi.org/10.1016/j.gca.2017.06.023. (WaterRock_ParamExplor code)
Neveu M., Rhoden A. (2017) The origin and evolution of a differentiated Mimas. Icarus296, 183-196. https://doi.org/10.1016/j.icarus.2017.06.011. (Tidal dissipation with basic orbital evolution driven by moon-primary interactions)
Neveu M., Desch S. (2015) Geochemistry, thermal evolution, and cryovolcanism on Ceres with a muddy ice mantle. Geophysical Research Letters 42, 10197-10206. http://dx.doi.org/10.1002/2015GL066375. (Retention of part of rock in icy mantle)
Neveu M., Desch S., Castillo-Rogez J. (2015) Core cracking and hydrothermal circulation profoundly affect Ceres' geophysical evolution. Journal of Geophysical Research: Planets 120, 123-154. http://dx.doi.org/10.1002/2014JE004714. (Thermal evolution code in C, cracking subroutine, convective transfer in rocky core by hydrothermal situation)
Neveu M., Desch S., Shock E., Glein C. (2015) Prerequisites for explosive cryovolcanism on dwarf planet-class Kuiper belt objects. Icarus 246, 48-64. http://dx.doi.org/10.1016/j.icarus.2014.03.043. (Cryolava code)
Rubin M., Desch S., Neveu M. (2014) The effect of Rayleigh-Taylor instabilities on thethickness of undifferentiated crusts on Kuiper belt objects. Icarus 236, 122-135. http://dx.doi.org/10.1016/j.icarus.2014.03.047. (Refined treatment of ice-rock differentiation)
Desch, S., Cook, J., Doggett, T., Porter, S. (2009) Thermal evolution of Kuiper belt objects, with implications for cryovolcanism. Icarus 202, 694-714. https://doi.org/10.1016/j.icarus.2009.03.009. (Thermal evolution code in Fortran)
Borderies, N., Goldreich, P. (1984) A simple derivation of capture probabilities for the J+1:J and J+2:J orbit-orbit resonance problems. Celestial Mechanics 32, 127-136. https://doi.org/10.1007/BF01231120.
Meyer, J., Wisdom, J. (2008) Tidal evolution of Mimas, Enceladus, and Dione. Icarus 193, 213-223. https://doi.org/10.1016/j.icarus.2007.09.008.
Segatz, M., Spohn, T., Ross, M. and Schubert, G. (1988) Tidal dissipation, surface heat flow, and figure of viscoelastic models of Io. Icarus 75, 187-206. https://doi.org/10.1016/0019-1035(88)90001-2.
crack software download Datakit Cross Manager v2014.1 Cliosoft SOS v5.31 Linux CHAM Phoenics v2011 Vector Fields CONCERTO v6.0anwer8#nextmail.ru---- change '#' to '@' New Software everyday Update,Anything you need,You can also check here: ctrl + f3DFlow 3DF Zephyr Pro 1.012 Win643Muri v10Arqcom.CAD-Earth.v4.1.7 for.AUTOCAD.2007-2016.ZWCAD.BricsCADAspen Technology aspenONE v8.8AspenONE Engineering suite v8.8Batch Plot DWG 2.4CADSWES.RiverWare.v6.6.6.Win32_64CATIA V5-6R2015 SP2 Win32_64CD-Adapco Star CCM+ 10.04.008-R4 (single precision)CD-Adapco Star CCM+ 10.04.008-R8 (double precision)cncKad V10.066cnckad v10.498 x32x64Itasca.FLAC3D v5.01DesignBuilder.Software.Ltd.DesignBuilder.v4.2.0.054DNV GL AS Phast v7.11.33.0exocad v2015Archicad 19 build 3003 Win64MacOSXPackEdge 14.0.1 & Plato 14.0.1BySoft7 V7 .2.0.0Molecular Operating Environment 2014.0901 WinMaclnxjason v9.0 2015IAR Embedded Workbench for ARM 7.40Kongsberg.LedaFlow.Engineering.v1.7.248.921Altera Quartus II v15.0 x64Laker.OA.vJ-2014.09-SP1-4.Linux64laker_adp_v2015.03laker_v2015.03-1Lumerical 2015b build 501 win3264linux64macMapInfo.Professional.v12.5.4.Build.402.Win64Maplesoft.Maple.v2015.1.Win32_64Mentor.Graphics.ModelSIM.SE.v10.4.Win64microsoft_dynamics_gp_2015_r2_x86_x64Mimics.Innovation.Suite.v17.1.Medical.Win32_64nPower_PowerSurfacing_v2.20.0075_for_SW2012-2015_64bitPolar Speedstack 2013 v13.02Scanvec Amiable Enroute v5.0Scanvec Amiable Enroute v5.1Schrodinger.KNIME.Workflows.2015-2Seer3D v2.10Siemens SolidEdge ST8 v108.00.00.091 English Win64SIEMENS_EPACTOOL_V3.24solidThinking Suite (Evolve + Inspire) 2015.4940 Win64Structural.Vibration.Solutions.ARTeMIS.Modal.v4.0.0.2Synopsys.CustomExplorer.vK-2015.06.Winlinuxst32Thunderhead Engineering PyroSim 2015.2.0604 Win64WinAC.RTX.2010AcornPipe.v8-619Altium Designer v15.1.12 win7&8Altium Vault 2.5.9ANSYS 16.1 nCode DesignLife Win64&Linux64Aquaveo Groundwater Modeling System v10.0.9 Win64Bentley CivilStorm (SELECTSeries 5) V8i 08.11.05.58Bentley SewerCAD (SELECTSeries 5) V8i 08.11.05.58Bentley SewerGEMS (SELECTSeries 5) V8i 08.11.05.58Bentley StormCAD (SELECTSeries 5) V8i 08.11.05.58CAE Studio v3.24.25.0CFTurbo v10.0.2CorelCAD.2015.5.v15.2.1.2037.Win32_64&MacCPFD Barracuda VR v17.0.0 Win64Delcam.Crispin.Engineer.2015.R1.SP3.Win32_64Delcam_Crispin_ShoeCost 2015_R2_SP4Delcam PowerShape 2015 R2_SP3Digital Canal SolidBuilder v21.2Digital.Canal.JobTracker.v4.10.227.4EPLAN.Platform.v2.4.Hotfix.2.Win32_64GibbsCAM 2015 v10.9.19.0 Win64Global.Mapper.v16.2.1.b052915.Win32_64GNS Animator4 v2.1.2 WinLinux x64HDL.Companion.v2.8.R1 winlinuxLectra.Diamino Fashion.v6R1.SP4 Lectra.Optiplan.V3R3.SP3Helmel Engineering Geomet v7.01.182IBM.SPSS.Statistics.and.AMOS.v23.0 Win32_64IBM.SPSS.Statistics.v23.LinuxInstaCode v2014Spectrum.Micro-Cap.v11.0.1.2.Win32Intergraph PVElite 2015 SP1 v17.00.01Keil C51 v9.54Keil C166 v7.55LMS Samtech Tea Pipe rev15SL2 for CATIA V5 R18-R22 Win32_64Mastercam X9 v18.0.11898.10 for SW Win64Mastercam_X9_v18.0.11898.0_Eng_Win64Materialise 3-Matic v9.0 x32x64Materialise Mimics.Innovation.Suite.v17.1.Medical.x32x64MicroSurvey.inCAD.Premium.2015.v15.0.0.2180Noesis Optimus v10.14 Win3264NUMECA FINE Open with OpenLabs v4.3 Win32_64Oasys Suite v12.1 WinLINUXOptimal.Cutting.Optimization.Pro.v5.9.8.10Plexim Plecs.Standalone.v3.6.7 Win3264&linux3264PointWise v17.3 R2 build 0185201603PTC Creo ElementsPro 5.0 M270 x32x64Ricardo Suite 2015.1 Windows + LinuxSafe.Software.FME.Desktop.v2015.1.0.1.15482.Win32_64See Electrical 7R2Schrodinger Suites v2015-2 WinMacLnx & KNIME WorkflowsSiemens.NX Nastran.v10.2.Win64.&.Linux64SolidCAM 2015 SP2 HF4 for SW2012-2015 Win32_64Synopsys.Fpga.Synthesis.vJ-2015.03.WinlinuxTHE_FOUNDRY_MODO_V901 win64linux64Thunderhead.Engineering.Pathfinder.v2015.1.0520Virtutech Simics 3.0.31 Linux32_64Xilinx.Vivado.Design.Suite.v2015.1Zuken E3.series 2015Synopsys Saber vJ-2015.03 Windows Synopsys Saber vJ-2015.03 Linux Forsk.Atoll.v3.3.0.7383.Win32_64Golden Software MapViewer 8.2.277IBM SPSS Statistics v23.0 x86x64Magic.Bullet.v12.0.3.for.FCPX.WinmacStruSoft.FEM-Design.Suite.v14.00.004VERO ALPHACAM 2015 R2Vero Edgecam 2015 R2XYZ.Scientific.TrueGrid.V3.1.2 Win32_64AeroHydro.SurfaceWorks.v8.7.392.0.Win32_64Ansys SpaceClaim 2015 SP1 Win32_64Coreldaw.Graphics.Suite.X7.5.Win32_64CSI Safe 2014 v14.1.0Exelis.ENVI.v5.2.SP1.Win32_64FoamWorks v4.0GemCad v1.09Global.Mapper.v16.2.1.Build.052915.x86x64IBM SPSS Modeler v14.1 Win32_64Inpho.UASMaster.v5.6.3KBC Petro-SIM v6.0 Win32_64KESZ.ConSteel.csJoint.v9.0.004LimitState FIX v2.0.0.380 x86x64 Maplesoft.Maple v2015.1 (MapleSim)Mentor Graphics Capital 2014.1 Win64Nuhertz Filter Solutions 2015 version 14.0OMRON CX-ONE v4.33PV SOL premium v7.5 R4QPS.Fledermaus.v7.4.4a.Win32_64Siemens LMS TecWare v3.10 Win32_64Structural.Vibration.Solutions.ARTeMIS.Modal.v4.0.0.2Valentin.PVSOL.Premium.v7.0.R5VeriSTAR Info VeriSTAR Hull v5.8.1 Win32_64Wilcom E2.0T Portable