Introduction to geodynamic modelling Introduction to DOUAR
|
|
- Lucas Evans
- 5 years ago
- Views:
Transcription
1 Introduction to geodynamic modelling Introduction to DOUAR David Whipp and Lars Kaislaniemi Department of Geosciences and Geography, Univ. Helsinki 1
2 Goals of this lecture Introduce DOUAR, a 3D thermomechanical numerical modelling program for creeping flows Present some of the important features in DOUAR that will be relevant for our use 2
3 DOUAR in a nutshell DOUAR (Braun et al., 2008) is a 3D finite-element code for modelling geodynamic processes It is designed to be run efficiently in parallel on computer clusters to be able to solve 3D problems at high spatial resolution We don t have time for a complete description of the numerical and computational aspects of DOUAR, but we ll see the highlights DOUAR is the word for Earth in the Breton language May 20,
4 <latexit sha1_base64="gxv2cpy2kibgecpa9euyg71uae0=">aaacs3icdvbnsxxbeo3zagjwe1dz9fk4cipkmzfacjag5jkjblefnxwp6e3dbeypobtgwib5x/kt/obccsgluecmhtkz7ki2hwvnv3qviqp6wa6kpzj+ejwerkw+fbb2vlm+8ellzmtr+8lbwnhr5vzzd5whf0oa0svjslzltqdolljmbj7u+uwtcf5ac07txpq1jo0csy4uqehru2owuwgph1qcvbdc/olkdhlbwses5ddlqpemtpfnvxuarxdqdu8hdrml40grhxeo4jrgb5b0zn/czvm4hbs+puplcy0mcyxe95i4p36jjirxomqmhrc58hsci16abrxw/xj2ewv7grncylrwdmgm/b2jro39vgehst7b/6nv5l+0xkgjd/1smrwgyfjjofghgczursjqosfjtqna7mtyffgehxikdi9nyxtvdkyslof/g4ujthj3krm37zpjut1rbiftsn2wslfshh1kp6zlopvmvrlv7ed0f/2m7qohx9jgno95xzaisfolp9oyrq==</latexit> <latexit sha1_base64="gxv2cpy2kibgecpa9euyg71uae0=">aaacs3icdvbnsxxbeo3zagjwe1dz9fk4cipkmzfacjag5jkjblefnxwp6e3dbeypobtgwib5x/kt/obccsgluecmhtkz7ki2hwvnv3qviqp6wa6kpzj+ejwerkw+fbb2vlm+8ellzmtr+8lbwnhr5vzzd5whf0oa0svjslzltqdolljmbj7u+uwtcf5ac07txpq1jo0csy4uqehru2owuwgph1qcvbdc/olkdhlbwses5ddlqpemtpfnvxuarxdqdu8hdrml40grhxeo4jrgb5b0zn/czvm4hbs+puplcy0mcyxe95i4p36jjirxomqmhrc58hsci16abrxw/xj2ewv7grncylrwdmgm/b2jro39vgehst7b/6nv5l+0xkgjd/1smrwgyfjjofghgczursjqosfjtqna7mtyffgehxikdi9nyxtvdkyslof/g4ujthj3krm37zpjut1rbiftsn2wslfshh1kp6zlopvmvrlv7ed0f/2m7qohx9jgno95xzaisfolp9oyrq==</latexit> <latexit sha1_base64="gxv2cpy2kibgecpa9euyg71uae0=">aaacs3icdvbnsxxbeo3zagjwe1dz9fk4cipkmzfacjag5jkjblefnxwp6e3dbeypobtgwib5x/kt/obccsgluecmhtkz7ki2hwvnv3qviqp6wa6kpzj+ejwerkw+fbb2vlm+8ellzmtr+8lbwnhr5vzzd5whf0oa0svjslzltqdolljmbj7u+uwtcf5ac07txpq1jo0csy4uqehru2owuwgph1qcvbdc/olkdhlbwses5ddlqpemtpfnvxuarxdqdu8hdrml40grhxeo4jrgb5b0zn/czvm4hbs+puplcy0mcyxe95i4p36jjirxomqmhrc58hsci16abrxw/xj2ewv7grncylrwdmgm/b2jro39vgehst7b/6nv5l+0xkgjd/1smrwgyfjjofghgczursjqosfjtqna7mtyffgehxikdi9nyxtvdkyslof/g4ujthj3krm37zpjut1rbiftsn2wslfshh1kp6zlopvmvrlv7ed0f/2m7qohx9jgno95xzaisfolp9oyrq==</latexit> <latexit sha1_base64="gxv2cpy2kibgecpa9euyg71uae0=">aaacs3icdvbnsxxbeo3zagjwe1dz9fk4cipkmzfacjag5jkjblefnxwp6e3dbeypobtgwib5x/kt/obccsgluecmhtkz7ki2hwvnv3qviqp6wa6kpzj+ejwerkw+fbb2vlm+8ellzmtr+8lbwnhr5vzzd5whf0oa0svjslzltqdolljmbj7u+uwtcf5ac07txpq1jo0csy4uqehru2owuwgph1qcvbdc/olkdhlbwses5ddlqpemtpfnvxuarxdqdu8hdrml40grhxeo4jrgb5b0zn/czvm4hbs+puplcy0mcyxe95i4p36jjirxomqmhrc58hsci16abrxw/xj2ewv7grncylrwdmgm/b2jro39vgehst7b/6nv5l+0xkgjd/1smrwgyfjjofghgczursjqosfjtqna7mtyffgehxikdi9nyxtvdkyslof/g4ujthj3krm37zpjut1rbiftsn2wslfshh1kp6zlopvmvrlv7ed0f/2m7qohx9jgno95xzaisfolp9oyrq==</latexit> Physics of DOUAR DOUAR calculates the flow of a highly viscous flow using the Stokes equation, which we have previously seen r (rv + rv T ) rp = g where η is the fluid viscosity, v is the fluid velocity, p is pressure, ρ is density, and g is the acceleration due to gravity. 4
5 <latexit sha1_base64="swjta1hfegwjvdasskhupi/5dti=">aaacgxicdvdlsgmxfm3uv62vuzdugkvw4zbtbf0ofny4rgaf0blkjpnpq5pmkgqkzejwn/ax3orenbh15dyvmdnwsd4ohhs4515u7gltrpv23xertls8srpwxq9sbg5t79i7ey2vzbktjk5yijshuorrqzqaaky6qssih4y0w+fv4bdhrcqaifs9tknauv/qmgkkjdsz4ynptheeoc9qyezbuakhh/j8nigxmo3zvdepuqxgb+i50+pwwrynnv3hrwnoobeam6ru13ntheriaoozmvt8tjeu4shqk66hangignx6yqqegswccslnexpo1e8toejkjxloojnsa/xtk8s/vg6m4/mgpylnnbf4tijognqjlgkbezueazy2bgfjzv8hhicjsdbhlwwj+ariqvm6hp5pwjxhcx3v5rrav5jhuwyh4bacaw+cgtq4bg3qbbjcgqfwcj6se+vzerfez60laz6zdxzgvx0crzqfwq==</latexit> <latexit sha1_base64="swjta1hfegwjvdasskhupi/5dti=">aaacgxicdvdlsgmxfm3uv62vuzdugkvw4zbtbf0ofny4rgaf0blkjpnpq5pmkgqkzejwn/ax3orenbh15dyvmdnwsd4ohhs4515u7gltrpv23xertls8srpwxq9sbg5t79i7ey2vzbktjk5yijshuorrqzqaaky6qssih4y0w+fv4bdhrcqaifs9tknauv/qmgkkjdsz4ynptheeoc9qyezbuakhh/j8nigxmo3zvdepuqxgb+i50+pwwrynnv3hrwnoobeam6ru13ntheriaoozmvt8tjeu4shqk66hangignx6yqqegswccslnexpo1e8toejkjxloojnsa/xtk8s/vg6m4/mgpylnnbf4tijognqjlgkbezueazy2bgfjzv8hhicjsdbhlwwj+ariqvm6hp5pwjxhcx3v5rrav5jhuwyh4bacaw+cgtq4bg3qbbjcgqfwcj6se+vzerfez60laz6zdxzgvx0crzqfwq==</latexit> <latexit sha1_base64="swjta1hfegwjvdasskhupi/5dti=">aaacgxicdvdlsgmxfm3uv62vuzdugkvw4zbtbf0ofny4rgaf0blkjpnpq5pmkgqkzejwn/ax3orenbh15dyvmdnwsd4ohhs4515u7gltrpv23xertls8srpwxq9sbg5t79i7ey2vzbktjk5yijshuorrqzqaaky6qssih4y0w+fv4bdhrcqaifs9tknauv/qmgkkjdsz4ynptheeoc9qyezbuakhh/j8nigxmo3zvdepuqxgb+i50+pwwrynnv3hrwnoobeam6ru13ntheriaoozmvt8tjeu4shqk66hangignx6yqqegswccslnexpo1e8toejkjxloojnsa/xtk8s/vg6m4/mgpylnnbf4tijognqjlgkbezueazy2bgfjzv8hhicjsdbhlwwj+ariqvm6hp5pwjxhcx3v5rrav5jhuwyh4bacaw+cgtq4bg3qbbjcgqfwcj6se+vzerfez60laz6zdxzgvx0crzqfwq==</latexit> <latexit sha1_base64="swjta1hfegwjvdasskhupi/5dti=">aaacgxicdvdlsgmxfm3uv62vuzdugkvw4zbtbf0ofny4rgaf0blkjpnpq5pmkgqkzejwn/ax3orenbh15dyvmdnwsd4ohhs4515u7gltrpv23xertls8srpwxq9sbg5t79i7ey2vzbktjk5yijshuorrqzqaaky6qssih4y0w+fv4bdhrcqaifs9tknauv/qmgkkjdsz4ynptheeoc9qyezbuakhh/j8nigxmo3zvdepuqxgb+i50+pwwrynnv3hrwnoobeam6ru13ntheriaoozmvt8tjeu4shqk66hangignx6yqqegswccslnexpo1e8toejkjxloojnsa/xtk8s/vg6m4/mgpylnnbf4tijognqjlgkbezueazy2bgfjzv8hhicjsdbhlwwj+ariqvm6hp5pwjxhcx3v5rrav5jhuwyh4bacaw+cgtq4bg3qbbjcgqfwcj6se+vzerfez60laz6zdxzgvx0crzqfwq==</latexit> <latexit sha1_base64="d4lq1wl/okwmhjnaosmql5e6cn4=">aaacehicdvdlssnafj3uv62vqasxbgal4cokrdcfqsgnywr2au0ok8mkhtozctotqgn5cx/bre7divv/wk1f4qst4ppamidz7uxee8kuuavd982qlc2vrk5v12sbm1vbo/buxkclmcskjrowyf6ifgfuklammpfekgniispdchxv+t0jkyom4lzpuxjwnbq0phhpiw3sa1+gkcho4yjr0a95pingjxqhdt11gm4j+jt4zux362cb1sb+96mez5wijrlsqu+5qq5yjdxfjbq1p1mkrximhqrvqeccqccfhvday6neme6keuldmfq1i0dcqskptsvheqr+eqx4l9fpdhwe5fskmsyczwffgym6gwuamkksym2mhiasqdkv4hgscgut2bcpis9qjptpy+h/pnnwpnfxbk7rzytfpfvwci7acfdagwica9acbybbae7ba3i07qwn69l6mzdwrexppvgg6/ud3eichq==</latexit> <latexit sha1_base64="d4lq1wl/okwmhjnaosmql5e6cn4=">aaacehicdvdlssnafj3uv62vqasxbgal4cokrdcfqsgnywr2au0ok8mkhtozctotqgn5cx/bre7divv/wk1f4qst4ppamidz7uxee8kuuavd982qlc2vrk5v12sbm1vbo/buxkclmcskjrowyf6ifgfuklammpfekgniispdchxv+t0jkyom4lzpuxjwnbq0phhpiw3sa1+gkcho4yjr0a95pingjxqhdt11gm4j+jt4zux362cb1sb+96mez5wijrlsqu+5qq5yjdxfjbq1p1mkrximhqrvqeccqccfhvday6neme6keuldmfq1i0dcqskptsvheqr+eqx4l9fpdhwe5fskmsyczwffgym6gwuamkksym2mhiasqdkv4hgscgut2bcpis9qjptpy+h/pnnwpnfxbk7rzytfpfvwci7acfdagwica9acbybbae7ba3i07qwn69l6mzdwrexppvgg6/ud3eichq==</latexit> <latexit sha1_base64="d4lq1wl/okwmhjnaosmql5e6cn4=">aaacehicdvdlssnafj3uv62vqasxbgal4cokrdcfqsgnywr2au0ok8mkhtozctotqgn5cx/bre7divv/wk1f4qst4ppamidz7uxee8kuuavd982qlc2vrk5v12sbm1vbo/buxkclmcskjrowyf6ifgfuklammpfekgniispdchxv+t0jkyom4lzpuxjwnbq0phhpiw3sa1+gkcho4yjr0a95pingjxqhdt11gm4j+jt4zux362cb1sb+96mez5wijrlsqu+5qq5yjdxfjbq1p1mkrximhqrvqeccqccfhvday6neme6keuldmfq1i0dcqskptsvheqr+eqx4l9fpdhwe5fskmsyczwffgym6gwuamkksym2mhiasqdkv4hgscgut2bcpis9qjptpy+h/pnnwpnfxbk7rzytfpfvwci7acfdagwica9acbybbae7ba3i07qwn69l6mzdwrexppvgg6/ud3eichq==</latexit> <latexit sha1_base64="d4lq1wl/okwmhjnaosmql5e6cn4=">aaacehicdvdlssnafj3uv62vqasxbgal4cokrdcfqsgnywr2au0ok8mkhtozctotqgn5cx/bre7divv/wk1f4qst4ppamidz7uxee8kuuavd982qlc2vrk5v12sbm1vbo/buxkclmcskjrowyf6ifgfuklammpfekgniispdchxv+t0jkyom4lzpuxjwnbq0phhpiw3sa1+gkcho4yjr0a95pingjxqhdt11gm4j+jt4zux362cb1sb+96mez5wijrlsqu+5qq5yjdxfjbq1p1mkrximhqrvqeccqccfhvday6neme6keuldmfq1i0dcqskptsvheqr+eqx4l9fpdhwe5fskmsyczwffgym6gwuamkksym2mhiasqdkv4hgscgut2bcpis9qjptpy+h/pnnwpnfxbk7rzytfpfvwci7acfdagwica9acbybbae7ba3i07qwn69l6mzdwrexppvgg6/ud3eichq==</latexit> Incompressible flow, almost It is further assumed that the fluid is nearly incompressible In an incompressible fluid the divergence of the velocity field is zero r v =0 In DOUAR, the slight compressibility of the fluid is used to determine fluid pressure (eliminating it from the Stokes equation) using the penalty method r v = p where λ is the penalty factor (typically 8 orders of magnitude larger than the shear viscosity) 5
6 <latexit sha1_base64="ijsrrn09blu8bhwxqidf9ieqzqe=">aaacmnicdvdlsgmxfm34tr6qlt0ei1axtjnfuefbconsxwqhrswt3mpojjmkd4plmg/xj/wft7rwnbr1i8zucj4vhhs45z5yjx9jydb1h52r0bhxicmp6dzm7nz8qn5x6cyesezq5aemdc1nbqrquewbemqrbhb4es797kgmn/dagxgqu+xh0azypridwrlaqpxfaqayukez1ercldbaisanhtmqgsfdlv4kxrmonrx1q8f1lbspy+rkdd1t5qtuqejmqx8drztiboem46ivf7gzeryaqi6zmxxpjbczmi2cs0hzjdhaxhixxuldqsucmm1kcgjk1yztpp1q26eqdtivhqkljokhvq0mgf6zn1pg/qxvy+xsnxohohhb8y9fnvhsdgnmf20ldrxl3wlgtbb/pfykacbruvptix+kowvk5+x0f3bwkxluytvelozvdu2ziitklrsjr7bipjkkr6rkolkhd+sepdi3zppz7lx+li44w55l8i2ct3fq36mo</latexit> <latexit sha1_base64="ijsrrn09blu8bhwxqidf9ieqzqe=">aaacmnicdvdlsgmxfm34tr6qlt0ei1axtjnfuefbconsxwqhrswt3mpojjmkd4plmg/xj/wft7rwnbr1i8zucj4vhhs45z5yjx9jydb1h52r0bhxicmp6dzm7nz8qn5x6cyesezq5aemdc1nbqrquewbemqrbhb4es797kgmn/dagxgqu+xh0azypridwrlaqpxfaqayukez1ercldbaisanhtmqgsfdlv4kxrmonrx1q8f1lbspy+rkdd1t5qtuqejmqx8drztiboem46ivf7gzeryaqi6zmxxpjbczmi2cs0hzjdhaxhixxuldqsucmm1kcgjk1yztpp1q26eqdtivhqkljokhvq0mgf6zn1pg/qxvy+xsnxohohhb8y9fnvhsdgnmf20ldrxl3wlgtbb/pfykacbruvptix+kowvk5+x0f3bwkxluytvelozvdu2ziitklrsjr7bipjkkr6rkolkhd+sepdi3zppz7lx+li44w55l8i2ct3fq36mo</latexit> <latexit sha1_base64="ijsrrn09blu8bhwxqidf9ieqzqe=">aaacmnicdvdlsgmxfm34tr6qlt0ei1axtjnfuefbconsxwqhrswt3mpojjmkd4plmg/xj/wft7rwnbr1i8zucj4vhhs45z5yjx9jydb1h52r0bhxicmp6dzm7nz8qn5x6cyesezq5aemdc1nbqrquewbemqrbhb4es797kgmn/dagxgqu+xh0azypridwrlaqpxfaqayukez1ercldbaisanhtmqgsfdlv4kxrmonrx1q8f1lbspy+rkdd1t5qtuqejmqx8drztiboem46ivf7gzeryaqi6zmxxpjbczmi2cs0hzjdhaxhixxuldqsucmm1kcgjk1yztpp1q26eqdtivhqkljokhvq0mgf6zn1pg/qxvy+xsnxohohhb8y9fnvhsdgnmf20ldrxl3wlgtbb/pfykacbruvptix+kowvk5+x0f3bwkxluytvelozvdu2ziitklrsjr7bipjkkr6rkolkhd+sepdi3zppz7lx+li44w55l8i2ct3fq36mo</latexit> <latexit sha1_base64="ijsrrn09blu8bhwxqidf9ieqzqe=">aaacmnicdvdlsgmxfm34tr6qlt0ei1axtjnfuefbconsxwqhrswt3mpojjmkd4plmg/xj/wft7rwnbr1i8zucj4vhhs45z5yjx9jydb1h52r0bhxicmp6dzm7nz8qn5x6cyesezq5aemdc1nbqrquewbemqrbhb4es797kgmn/dagxgqu+xh0azypridwrlaqpxfaqayukez1ercldbaisanhtmqgsfdlv4kxrmonrx1q8f1lbspy+rkdd1t5qtuqejmqx8drztiboem46ivf7gzeryaqi6zmxxpjbczmi2cs0hzjdhaxhixxuldqsucmm1kcgjk1yztpp1q26eqdtivhqkljokhvq0mgf6zn1pg/qxvy+xsnxohohhb8y9fnvhsdgnmf20ldrxl3wlgtbb/pfykacbruvptix+kowvk5+x0f3bwkxluytvelozvdu2ziitklrsjr7bipjkkr6rkolkhd+sepdi3zppz7lx+li44w55l8i2ct3fq36mo</latexit> Rheologies Materials in DOUAR are either viscous or plastic (no elasticity) Nonlinear viscosity is modelled using the equation for temperature-dependent nonlinear viscosity = 0 " 1/(n 1) exp (Q/nRT ). where η 0 is the viscosity pre-factor, ε is the strain rate, n is the power law exponent, Q is the activation energy, R is the universal gas constant, and T is temperature in Kelvins 6
7 <latexit sha1_base64="uwaube5dhfpds3uemonct4aykzi=">aaacn3icdvbnaxsxenwmtzqpjnhbyy6ijhakmbsh0b5scptsywp1evaamyupbwfjk6rrwsz7h/jrcj0lf6on3kkuvfrcrenc3tqphb7vztazr7bkekrt78nkk6era8/wnza3nm/v7lzevdz3zxacu6jupbsswkosbrsksegldqi6uhhrtd80/suxdf6w5jpnlpy1ji0csqeupuhrtr4sqedvuechppdyrgfq5rpo4nrl6prnbibndiihrxbaouob8ick68z/tm0wobu0fuxdugsnhoqc73tzaqlfgsmpfnabefboquxhjl1idwj0/wp+u833ozlko9lfz4jp1b87ktdez3qrk5tl/b9ei/7p6wuavetx0thaamt9offqnerebmsh0qegnysehjnxvy4m4ebqjhfpsqgxbqichtvear2d+pmgf5ych3wytjn9om6fniwiw2d77du7ybl7y07zr3bgukywk/avxbob5fvyi7ln7u5lv5jfzyu2hotnb0cnrta=</latexit> <latexit sha1_base64="uwaube5dhfpds3uemonct4aykzi=">aaacn3icdvbnaxsxenwmtzqpjnhbyy6ijhakmbsh0b5scptsywp1evaamyupbwfjk6rrwsz7h/jrcj0lf6on3kkuvfrcrenc3tqphb7vztazr7bkekrt78nkk6era8/wnza3nm/v7lzevdz3zxacu6jupbsswkosbrsksegldqi6uhhrtd80/suxdf6w5jpnlpy1ji0csqeupuhrtr4sqedvuechppdyrgfq5rpo4nrl6prnbibndiihrxbaouob8ick68z/tm0wobu0fuxdugsnhoqc73tzaqlfgsmpfnabefboquxhjl1idwj0/wp+u833ozlko9lfz4jp1b87ktdez3qrk5tl/b9ei/7p6wuavetx0thaamt9offqnerebmsh0qegnysehjnxvy4m4ebqjhfpsqgxbqichtvear2d+pmgf5ych3wytjn9om6fniwiw2d77du7ybl7y07zr3bgukywk/avxbob5fvyi7ln7u5lv5jfzyu2hotnb0cnrta=</latexit> <latexit sha1_base64="uwaube5dhfpds3uemonct4aykzi=">aaacn3icdvbnaxsxenwmtzqpjnhbyy6ijhakmbsh0b5scptsywp1evaamyupbwfjk6rrwsz7h/jrcj0lf6on3kkuvfrcrenc3tqphb7vztazr7bkekrt78nkk6era8/wnza3nm/v7lzevdz3zxacu6jupbsswkosbrsksegldqi6uhhrtd80/suxdf6w5jpnlpy1ji0csqeupuhrtr4sqedvuechppdyrgfq5rpo4nrl6prnbibndiihrxbaouob8ick68z/tm0wobu0fuxdugsnhoqc73tzaqlfgsmpfnabefboquxhjl1idwj0/wp+u833ozlko9lfz4jp1b87ktdez3qrk5tl/b9ei/7p6wuavetx0thaamt9offqnerebmsh0qegnysehjnxvy4m4ebqjhfpsqgxbqichtvear2d+pmgf5ych3wytjn9om6fniwiw2d77du7ybl7y07zr3bgukywk/avxbob5fvyi7ln7u5lv5jfzyu2hotnb0cnrta=</latexit> <latexit sha1_base64="uwaube5dhfpds3uemonct4aykzi=">aaacn3icdvbnaxsxenwmtzqpjnhbyy6ijhakmbsh0b5scptsywp1evaamyupbwfjk6rrwsz7h/jrcj0lf6on3kkuvfrcrenc3tqphb7vztazr7bkekrt78nkk6era8/wnza3nm/v7lzevdz3zxacu6jupbsswkosbrsksegldqi6uhhrtd80/suxdf6w5jpnlpy1ji0csqeupuhrtr4sqedvuechppdyrgfq5rpo4nrl6prnbibndiihrxbaouob8ick68z/tm0wobu0fuxdugsnhoqc73tzaqlfgsmpfnabefboquxhjl1idwj0/wp+u833ozlko9lfz4jp1b87ktdez3qrk5tl/b9ei/7p6wuavetx0thaamt9offqnerebmsh0qegnysehjnxvy4m4ebqjhfpsqgxbqichtvear2d+pmgf5ych3wytjn9om6fniwiw2d77du7ybl7y07zr3bgukywk/avxbob5fvyi7ln7u5lv5jfzyu2hotnb0cnrta=</latexit> Rheologies There are several options for different plasticity criteria The Mohr-Coulomb criterion is what we use most often = c n tan where τ is the shear stress, c is the cohesion, σ n is the normal stress, and φ is the internal angle of friction 7
8 <latexit sha1_base64="uuf6x3qlzjncpjyncxx0dwz/hqi=">aaace3icdvfdaxnbfj3dvq31k+qjlxedud8iu0xrb4wcl32sklsftaizk7vjkplyzu4wwrl/0hf/ia+c4gysuuphhweo59zdvxomqlqklgxfknrv/9btowd3d+/df/dwue/xk/pgai9xjj12/riqabwyocjfgi8rj8iugi+k5edov7hch5szq1pvodfiblwppkbitxuw+4udcvxjsufasy9kwyvhsqknw/yguwuvgremuwqby5kj4fyuwsaquffzbb2et9furl/edernn+d02utng+osk/gb5ip1nfxzts6mve985mrt0jluiorxnlu0abqdpmb2kncbkygxyo7jck0wgcbnopcwxkrmbqxz8vicnfu7oxemhjupyqcrtah/ah35l21cu/lh0ihb1yrwbgavdcziqrcyzjrhsxovgzbexv1blksmlujx7ewpzm4bmrqae8rlg4o6tgp+d86pb3k2yl+87z983ez2wj6x5+yi5ew9o2gn7iynmgrf2y9kl9lpfqb99fx6ztoajlvpu7zt6btf/pzaxq==</latexit> <latexit sha1_base64="uuf6x3qlzjncpjyncxx0dwz/hqi=">aaace3icdvfdaxnbfj3dvq31k+qjlxedud8iu0xrb4wcl32sklsftaizk7vjkplyzu4wwrl/0hf/ia+c4gysuuphhweo59zdvxomqlqklgxfknrv/9btowd3d+/df/dwue/xk/pgai9xjj12/riqabwyocjfgi8rj8iugi+k5edov7hch5szq1pvodfiblwppkbitxuw+4udcvxjsufasy9kwyvhsqknw/yguwuvgremuwqby5kj4fyuwsaquffzbb2et9furl/edernn+d02utng+osk/gb5ip1nfxzts6mve985mrt0jluiorxnlu0abqdpmb2kncbkygxyo7jck0wgcbnopcwxkrmbqxz8vicnfu7oxemhjupyqcrtah/ah35l21cu/lh0ihb1yrwbgavdcziqrcyzjrhsxovgzbexv1blksmlujx7ewpzm4bmrqae8rlg4o6tgp+d86pb3k2yl+87z983ez2wj6x5+yi5ew9o2gn7iynmgrf2y9kl9lpfqb99fx6ztoajlvpu7zt6btf/pzaxq==</latexit> <latexit sha1_base64="uuf6x3qlzjncpjyncxx0dwz/hqi=">aaace3icdvfdaxnbfj3dvq31k+qjlxedud8iu0xrb4wcl32sklsftaizk7vjkplyzu4wwrl/0hf/ia+c4gysuuphhweo59zdvxomqlqklgxfknrv/9btowd3d+/df/dwue/xk/pgai9xjj12/riqabwyocjfgi8rj8iugi+k5edov7hch5szq1pvodfiblwppkbitxuw+4udcvxjsufasy9kwyvhsqknw/yguwuvgremuwqby5kj4fyuwsaquffzbb2et9furl/edernn+d02utng+osk/gb5ip1nfxzts6mve985mrt0jluiorxnlu0abqdpmb2kncbkygxyo7jck0wgcbnopcwxkrmbqxz8vicnfu7oxemhjupyqcrtah/ah35l21cu/lh0ihb1yrwbgavdcziqrcyzjrhsxovgzbexv1blksmlujx7ewpzm4bmrqae8rlg4o6tgp+d86pb3k2yl+87z983ez2wj6x5+yi5ew9o2gn7iynmgrf2y9kl9lpfqb99fx6ztoajlvpu7zt6btf/pzaxq==</latexit> <latexit sha1_base64="uuf6x3qlzjncpjyncxx0dwz/hqi=">aaace3icdvfdaxnbfj3dvq31k+qjlxedud8iu0xrb4wcl32sklsftaizk7vjkplyzu4wwrl/0hf/ia+c4gysuuphhweo59zdvxomqlqklgxfknrv/9btowd3d+/df/dwue/xk/pgai9xjj12/riqabwyocjfgi8rj8iugi+k5edov7hch5szq1pvodfiblwppkbitxuw+4udcvxjsufasy9kwyvhsqknw/yguwuvgremuwqby5kj4fyuwsaquffzbb2et9furl/edernn+d02utng+osk/gb5ip1nfxzts6mve985mrt0jluiorxnlu0abqdpmb2kncbkygxyo7jck0wgcbnopcwxkrmbqxz8vicnfu7oxemhjupyqcrtah/ah35l21cu/lh0ihb1yrwbgavdcziqrcyzjrhsxovgzbexv1blksmlujx7ewpzm4bmrqae8rlg4o6tgp+d86pb3k2yl+87z983ez2wj6x5+yi5ew9o2gn7iynmgrf2y9kl9lpfqb99fx6ztoajlvpu7zt6btf/pzaxq==</latexit> Thermal model The full 3D advection-diffusion equation with heat production is solved in + v rt = r krt + H where c is the heat capacity, T is temperature, t is time, k is the thermal conductivity, and H is the heat production per unit mass. 8
9 Numerical approach Braun et al., 2008 The finite-element mesh used in DOUAR is Eulerian and based on the octree division of space An Eulerian mesh is one that is fixed in space with respect to the fluid flowing within/through it The octree division of space is based on subdivisions of a unit cube 9
10 Numerical approach Braun et al., 2008 The unit cube is at the octree level 0 resolution, meaning that there are 2 0 elements along the width of the cube A mesh at octree level 6 would have 2 6, or 64 x 64 x 64 elements The initial octree resolution is an important setting in the input file 10
11 The octree division of space Level 2 Level 1 As an example, here is a face of a DOUAR model (2D, not full cube) at octree level 2 This is 2 subdivisions of the unit cube (red) Level 0 11
12 Other aspects of DOUAR (continued on whiteboard) 12
13 References Braun, J., Thieulot, C., Fullsack, P., DeKool, M., Beaumont, C. and Huismans, R., DOUAR: A new three-dimensional creeping flow numerical model for the solution of geological problems. Physics of the Earth and Planetary Interiors, 171(1-4), pp
A Novel Approach to High Speed Collision
A Novel Approach to High Speed Collision Avril Slone University of Greenwich Motivation High Speed Impact Currently a very active research area. Generic projectile- target collision 11 th September 2001.
More informationAnimation of Fluids. Animating Fluid is Hard
Animation of Fluids Animating Fluid is Hard Too complex to animate by hand Surface is changing very quickly Lots of small details In short, a nightmare! Need automatic simulations AdHoc Methods Some simple
More informationIntroduction to Finite Element Modelling in Geosciences
Spring Semester 2012 651-4144-00L Introduction to Finite Element Modelling in Geosciences August 5-9, 2013 Dave May (dave.may@erdw.ethz.ch) Marcel Frehner (marcel.frehner@erdw.ethz.ch) Thomas Philippe
More information2.7 Cloth Animation. Jacobs University Visualization and Computer Graphics Lab : Advanced Graphics - Chapter 2 123
2.7 Cloth Animation 320491: Advanced Graphics - Chapter 2 123 Example: Cloth draping Image Michael Kass 320491: Advanced Graphics - Chapter 2 124 Cloth using mass-spring model Network of masses and springs
More informationSNAC: a tutorial. Eunseo Choi. July 29, Lamont-Doherty Earth Observatory
SNAC: a tutorial Eunseo Choi Lamont-Doherty Earth Observatory July 29, 2010 What is SNAC? StGermaiN Analysis of Continua 2 of 20 What is SNAC? StGermaiN Analysis of Continua 3-D Version of FLAC 2 of 20
More informationALE and Fluid-Structure Interaction in LS-DYNA March 2004
ALE and Fluid-Structure Interaction in LS-DYNA March 2004 Workshop Models 1. Taylor bar impact 2. One-dimensional advection test 3. Channel 4. Underwater explosion 5. Bar impacting water surface 6. Sloshing
More informationTrue 3D CAE visualization of filling imbalance in geometry-balanced runners
True 3D CAE visualization of filling imbalance in geometry-balanced runners C.C. Chien, * C.C. Chiang, W. H. Yang, Vito Tsai and David C.Hsu CoreTech System Co.,Ltd., HsinChu, Taiwan, ROC Abstract The
More informationCHAPTER 3. Elementary Fluid Dynamics
CHAPTER 3. Elementary Fluid Dynamics - Understanding the physics of fluid in motion - Derivation of the Bernoulli equation from Newton s second law Basic Assumptions of fluid stream, unless a specific
More informationSSR Polygonal Search Area
SSR Polygonal Search Area 22-1 SSR Polygonal Search Area In this tutorial, Phase2 is used to determine the factor of safety of a slope using the shear strength reduction (SSR) method. The SSR Polygon Search
More informationOP-029. INTERPRETATION OF CIU TEST
INTERPRETATION OF CIU TEST Page 1 of 7 WORK INSTRUCTIONS FOR ENGINEERS KYW Compiled by : Checked by LCH : TYC Approved by : OP-09. INTERPRETATION OF CIU TEST INTERPRETATION OF CIU TEST Page of 7 9. INTERPRETATION
More informationFEMLAB Exercise 1 for ChE366
FEMLAB Exercise 1 for ChE366 Problem statement Consider a spherical particle of radius r s moving with constant velocity U in an infinitely long cylinder of radius R that contains a Newtonian fluid. Let
More informationCS-184: Computer Graphics Lecture #21: Fluid Simulation II
CS-184: Computer Graphics Lecture #21: Fluid Simulation II Rahul Narain University of California, Berkeley Nov. 18 19, 2013 Grid-based fluid simulation Recap: Eulerian viewpoint Grid is fixed, fluid moves
More informationCGT 581 G Fluids. Overview. Some terms. Some terms
CGT 581 G Fluids Bedřich Beneš, Ph.D. Purdue University Department of Computer Graphics Technology Overview Some terms Incompressible Navier-Stokes Boundary conditions Lagrange vs. Euler Eulerian approaches
More informationVerification and Validation of Turbulent Flow around a Clark-Y Airfoil
Verification and Validation of Turbulent Flow around a Clark-Y Airfoil 1. Purpose 58:160 Intermediate Mechanics of Fluids CFD LAB 2 By Tao Xing and Fred Stern IIHR-Hydroscience & Engineering The University
More informationIsotropic Porous Media Tutorial
STAR-CCM+ User Guide 3927 Isotropic Porous Media Tutorial This tutorial models flow through the catalyst geometry described in the introductory section. In the porous region, the theoretical pressure drop
More informationComputation of Velocity, Pressure and Temperature Distributions near a Stagnation Point in Planar Laminar Viscous Incompressible Flow
Excerpt from the Proceedings of the COMSOL Conference 8 Boston Computation of Velocity, Pressure and Temperature Distributions near a Stagnation Point in Planar Laminar Viscous Incompressible Flow E. Kaufman
More informationGeodynamic models 20 years past and future
www.sci.monash.edu Geodynamic models 20 years past and future Louis Moresi, School of Earth Sciences, University of Melbourne 1994 1999 2004 2009 2014 1 Plate tectonics a multi-scale problem 2 Integrated
More informationTHERMAL EXPANSION OF A NAVIGABLE LOCK
THERMAL EXPANSION OF A NAVIGABLE LOCK 15 THERMAL EXPANSION OF A NAVIGABLE LOCK A navigable lock is temporarily 'empty' due to maintenance. After some time there is significant increase of the air temperature,
More informationGeoStudio Product Details
Product Details I. Purchasing Options Universal Professional Standard Basic Student (Free) Feature Comparison Finite Element Integration No limit 500 elements 500 elements Analysis Ordinary Yes Yes Yes
More informationNUMERICAL INVESTIGATION OF THE FLOW BEHAVIOR INTO THE INLET GUIDE VANE SYSTEM (IGV)
University of West Bohemia» Department of Power System Engineering NUMERICAL INVESTIGATION OF THE FLOW BEHAVIOR INTO THE INLET GUIDE VANE SYSTEM (IGV) Publication was supported by project: Budování excelentního
More informationPhase-field simulation of two-phase micro-flows in a Hele-Shaw cell
Computational Methods in Multiphase Flow III 7 Phase-field simulation of two-phase micro-flows in a Hele-Shaw cell Y. Sun & C. Beckermann Department of Mechanical and Industrial Engineering, University
More informationPrinciples and Goals
Simulations of foam rheology Simon Cox Principles and Goals Use knowledge of foam structure to predict response when experiments are not easy to isolate causes of certain effects to save time in designing
More informationDynamic Computational Modeling of the Glass Container Forming Process
Dynamic Computational Modeling of the Glass Container Forming Process Matthew Hyre 1, Ryan Taylor, and Morgan Harris Virginia Military Institute, Lexington, Virginia, USA Abstract Recent advances in numerical
More informationDevelopment of Hybrid Fluid Jet / Float Polishing Process
COMSOL Conference - Tokyo 2013 Development of Hybrid Fluid Jet / Float Polishing Process A. Beaucamp, Y. Namba Dept. of Mechanical Engineering, Chubu University, Japan Zeeko LTD, United Kingdom Research
More informationA numerical model of mass movements over complex topography
massmov2d A numerical model of mass movements over complex topography Description massmov2d is a simple fluid mechanics model for simulating the kinematics (run out and deposition) of fluid-like mass movements
More informationAbstract. Die Geometry. Introduction. Mesh Partitioning Technique for Coextrusion Simulation
OPTIMIZATION OF A PROFILE COEXTRUSION DIE USING A THREE-DIMENSIONAL FLOW SIMULATION SOFTWARE Kim Ryckebosh 1 and Mahesh Gupta 2, 3 1. Deceuninck nv, BE-8830 Hooglede-Gits, Belgium 2. Michigan Technological
More informationUsing the Discrete Ordinates Radiation Model
Tutorial 6. Using the Discrete Ordinates Radiation Model Introduction This tutorial illustrates the set up and solution of flow and thermal modelling of a headlamp. The discrete ordinates (DO) radiation
More informationFluid Simulation. [Thürey 10] [Pfaff 10] [Chentanez 11]
Fluid Simulation [Thürey 10] [Pfaff 10] [Chentanez 11] 1 Computational Fluid Dynamics 3 Graphics Why don t we just take existing models from CFD for Computer Graphics applications? 4 Graphics Why don t
More informationSimulation of RF HEat Test
Simulation of RF HEat Test Date: Tuesday, December 22, 2015 Designer: Solidworks Study name: Stress One Third Emissivity Analysis type: Nonlinear - Dynamic Description No Data Table of Contents Description...
More informationA brief description of the particle finite element method (PFEM2). Extensions to free surface
A brief description of the particle finite element method (PFEM2). Extensions to free surface flows. Juan M. Gimenez, L.M. González, CIMEC Universidad Nacional del Litoral (UNL) Santa Fe, Argentina Universidad
More informationENGR Socolofsky
ENGR 102-213 - Socolofsky Date : due 9/3/2018 at 12:40 p.m. Engineering Lab I - Computation Lab Assignment #01b Writing Your Own Programs Return your solution (one per group) as outlined in the activities
More informationPresentation of PAM-CRASH v2004. Part 1: Solver News
Presentation of PAM-CRASH v2004 Part 1: Solver News. 1 Overview New Options Elements Materials Others Quality Numerical precision and robustness 2 CRASH/SAFE 2G Evolution V2002: Basic reengineering Common
More informationNavier-Stokes & Flow Simulation
Last Time? Navier-Stokes & Flow Simulation Implicit Surfaces Marching Cubes/Tetras Collision Detection & Response Conservative Bounding Regions backtracking fixing Today Flow Simulations in Graphics Flow
More informationThe 3D DSC in Fluid Simulation
The 3D DSC in Fluid Simulation Marek K. Misztal Informatics and Mathematical Modelling, Technical University of Denmark mkm@imm.dtu.dk DSC 2011 Workshop Kgs. Lyngby, 26th August 2011 Governing Equations
More informationParticle-based Fluid Simulation
Simulation in Computer Graphics Particle-based Fluid Simulation Matthias Teschner Computer Science Department University of Freiburg Application (with Pixar) 10 million fluid + 4 million rigid particles,
More informationChapter 5: Introduction to Differential Analysis of Fluid Motion
Chapter 5: Introduction to Differential 5-1 Conservation of Mass 5-2 Stream Function for Two-Dimensional 5-3 Incompressible Flow 5-4 Motion of a Fluid Particle (Kinematics) 5-5 Momentum Equation 5-6 Computational
More informationPHYSICALLY BASED ANIMATION
PHYSICALLY BASED ANIMATION CS148 Introduction to Computer Graphics and Imaging David Hyde August 2 nd, 2016 WHAT IS PHYSICS? the study of everything? WHAT IS COMPUTATION? the study of everything? OUTLINE
More informationMAE 3130: Fluid Mechanics Lecture 5: Fluid Kinematics Spring Dr. Jason Roney Mechanical and Aerospace Engineering
MAE 3130: Fluid Mechanics Lecture 5: Fluid Kinematics Spring 2003 Dr. Jason Roney Mechanical and Aerospace Engineering Outline Introduction Velocity Field Acceleration Field Control Volume and System Representation
More informationPredicting Tumour Location by Modelling the Deformation of the Breast using Nonlinear Elasticity
Predicting Tumour Location by Modelling the Deformation of the Breast using Nonlinear Elasticity November 8th, 2006 Outline Motivation Motivation Motivation for Modelling Breast Deformation Mesh Generation
More informationSmoke Simulation using Smoothed Particle Hydrodynamics (SPH) Shruti Jain MSc Computer Animation and Visual Eects Bournemouth University
Smoke Simulation using Smoothed Particle Hydrodynamics (SPH) Shruti Jain MSc Computer Animation and Visual Eects Bournemouth University 21st November 2014 1 Abstract This report is based on the implementation
More informationSPC 307 Aerodynamics. Lecture 1. February 10, 2018
SPC 307 Aerodynamics Lecture 1 February 10, 2018 Sep. 18, 2016 1 Course Materials drahmednagib.com 2 COURSE OUTLINE Introduction to Aerodynamics Review on the Fundamentals of Fluid Mechanics Euler and
More informationParallelization study of a VOF/Navier-Stokes model for 3D unstructured staggered meshes
Parallelization study of a VOF/Navier-Stokes model for 3D unstructured staggered meshes L. Jofre, O. Lehmkuhl, R. Borrell, J. Castro and A. Oliva Corresponding author: cttc@cttc.upc.edu Centre Tecnològic
More informationNon-Newtonian Transitional Flow in an Eccentric Annulus
Tutorial 8. Non-Newtonian Transitional Flow in an Eccentric Annulus Introduction The purpose of this tutorial is to illustrate the setup and solution of a 3D, turbulent flow of a non-newtonian fluid. Turbulent
More informationExample 13 - Shock Tube
Example 13 - Shock Tube Summary This famous experiment is interesting for observing the shock-wave propagation. Moreover, this case uses the representation of perfect gas and compares the different formulations:
More informationSimulation of engraving process of large-caliber artillery using coupled Eulerian-Lagrangian method
Simulation of engraving process of large-caliber artillery using coupled Eulerian-Lagrangian method Zhen Li 1, Jianli Ge 2, Guolai Yang 3, Jun Tang 4 School of Mechanical Engineering, Nanjing University
More informationNumerical Investigation of Non-Newtonian Laminar Flow in Curved Tube with Insert
Numerical Investigation of Non-Newtonian Laminar Flow in Curved Tube with Insert A. Kadyyrov 1 1 Research center for power engineering problems Federal government budgetary institution of science Kazan
More informationNUMERICAL METHOD FOR FREE SURFACE VISCOUS FLOWS
STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume LIII, Number 3, September 2008 NUMERICAL METHOD FOR FREE SURFACE VISCOUS FLOWS TIBERIU IOANA AND TITUS PETRILA Abstract. In this paper we present a new algorithm
More informationComparison Between Numerical & PIV Experimental Results for Gas-Solid Flow in Ducts
Fabio Kasper Comparison Between Numerical & PIV Experimental Results for Gas-Solid Flow in Ducts Rodrigo Decker, Oscar Sgrott Jr., Henry F. Meier Waldir Martignoni Agenda Introduction The Test Bench Case
More informationNumerical Analysis of Shock Tube Problem by using TVD and ACM Schemes
Numerical Analysis of Shock Tube Problem by using TVD and Schemes Dr. Mukkarum Husain, Dr. M. Nauman Qureshi, Syed Zaid Hasany IST Karachi, Email: mrmukkarum@yahoo.com Abstract Computational Fluid Dynamics
More informationKEY WORDS: signed distance function, no-penetration condition, mesh adaptation, anisotropic mesh, level-set methods.
Using a signed distance function for the simulation of metal forming processes: formulation of the contact condition and mesh adaptation. From a Lagrangian approach to an Eulerian approach. J. Bruchon,
More informationSand Simulation. Abhinav Golas. COMP Physically Based Simulation. Final Project Presentation. May 6,
Sand Simulation Abhinav Golas COMP 768 - Physically Based Simulation Final Project Presentation May 6, 2009 1 Motivation Movies, games Spiderman 3 Engineering design grain silos Avalanches, Landslides
More informationCFD in COMSOL Multiphysics
CFD in COMSOL Multiphysics Christian Wollblad Copyright 2017 COMSOL. Any of the images, text, and equations here may be copied and modified for your own internal use. All trademarks are the property of
More informationModeling of Granular Materials
Modeling of Granular Materials Abhinav Golas COMP 768 - Physically Based Simulation April 23, 2009 1 Motivation Movies, games Spiderman 3 Engineering design grain silos Avalanches, Landslides The Mummy
More informationCompressible Subsonic Flow on a Staggered Grid
Compressible Subsonic Flow on a Staggered Grid by Michael Patrick Bonner B.Sc., California Polytechnic State University, San Luis Obispo, 00 A THESIS SUBMITTED IN PARTIAL FULFILMENT OF THE REQUIREMENTS
More informationVOLCANIC DEFORMATION MODELLING: NUMERICAL BENCHMARKING WITH COMSOL
VOLCANIC DEFORMATION MODELLING: NUMERICAL BENCHMARKING WITH COMSOL The following is a description of the model setups and input/output parameters for benchmarking analytical volcanic deformation models
More informationCFD modelling of thickened tailings Final project report
26.11.2018 RESEM Remote sensing supporting surveillance and operation of mines CFD modelling of thickened tailings Final project report Lic.Sc.(Tech.) Reeta Tolonen and Docent Esa Muurinen University of
More informationENHANCED STRUCTURE CAE SOLUTION WITH MOLDING EFFECT FOR AUTOMOTIVE PARTS
ENHANCED STRUCTURE CAE SOLUTION WITH MOLDING EFFECT FOR AUTOMOTIVE PARTS Allen Y. Peng*, Wen-Hsien Yang, David C. Hsu CoreTech System Co., Ltd., HsinChu, Taiwan, ROC Abstract An increasing number of automotive
More informationTerrain settlement analysis
Engineering manual No. 21 Updated: 02/2018 Terrain settlement analysis Program: File: FEM Demo_manual_21.gmk This example contains the solution to terrain settlement under surcharge loading using the Finite
More informationExample 24 Spring-back
Example 24 Spring-back Summary The spring-back simulation of sheet metal bent into a hat-shape is studied. The problem is one of the famous tests from the Numisheet 93. As spring-back is generally a quasi-static
More informationModeling of Granular Materials
Modeling of Granular Materials Abhinav Golas COMP 768 - Physically Based Simulation April 23, 2009 1 Motivation Movies, games Spiderman 3 Engineering design grain silos Avalanches, Landslides The Mummy
More informationDevelopment of New Method for Flow Computations in Vehicle Ventilation
2005:110 CIV MASTER S THESIS Development of New Method for Flow Computations in Vehicle Ventilation FRIDA NORDIN MASTER OF SCIENCE PROGRAMME Luleå University of Technology Department of Applied Physics
More informationLagrangian and Eulerian Representations of Fluid Flow: Kinematics and the Equations of Motion
Lagrangian and Eulerian Representations of Fluid Flow: Kinematics and the Equations of Motion James F. Price Woods Hole Oceanographic Institution Woods Hole, MA, 02543 July 31, 2006 Summary: This essay
More information1 Stokes equations with FD on a staggered grid using the stream-function approach.
Figure 1: Discretiation for the streamfunction approach. The boundary conditions are set through fictious boundary points. 1 Stokes equations with FD on a staggered grid using the stream-function approach.
More informationRealistic Animation of Fluids
1 Realistic Animation of Fluids Nick Foster and Dimitris Metaxas Presented by Alex Liberman April 19, 2005 2 Previous Work Used non physics-based methods (mostly in 2D) Hard to simulate effects that rely
More informationIntroduction to CFX. Workshop 2. Transonic Flow Over a NACA 0012 Airfoil. WS2-1. ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.
Workshop 2 Transonic Flow Over a NACA 0012 Airfoil. Introduction to CFX WS2-1 Goals The purpose of this tutorial is to introduce the user to modelling flow in high speed external aerodynamic applications.
More informationEssay 1: Dimensional Analysis of Models and Data Sets: Similarity Solutions
Table of Contents Essay 1: Dimensional Analysis of Models and Data Sets: Similarity Solutions and Scaling Analysis 1 About dimensional analysis 4 1.1 Thegoalandtheplan... 4 1.2 Aboutthisessay... 5 2 Models
More informationParallel High-Order Geometric Multigrid Methods on Adaptive Meshes for Highly Heterogeneous Nonlinear Stokes Flow Simulations of Earth s Mantle
ICES Student Forum The University of Texas at Austin, USA November 4, 204 Parallel High-Order Geometric Multigrid Methods on Adaptive Meshes for Highly Heterogeneous Nonlinear Stokes Flow Simulations of
More informationTHREE-DIMENSIONAL CAE OF WIRE-SWEEP IN MICROCHIP ENCAPSULATION
THREE-DIMENSIONAL CAE OF WIRE-SWEEP IN MICROCHIP ENCAPSULATION Wen-Hsien Yang *, David C.Hsu and Venny Yang CoreTech System Co.,Ltd., HsinChu, Taiwan, ROC Rong-Yeu Chang National Tsing-Hua University,
More informationCoupled Analysis of FSI
Coupled Analysis of FSI Qin Yin Fan Oct. 11, 2008 Important Key Words Fluid Structure Interface = FSI Computational Fluid Dynamics = CFD Pressure Displacement Analysis = PDA Thermal Stress Analysis = TSA
More informationParticle-Based Fluid Simulation. CSE169: Computer Animation Steve Rotenberg UCSD, Spring 2016
Particle-Based Fluid Simulation CSE169: Computer Animation Steve Rotenberg UCSD, Spring 2016 Del Operations Del: = x Gradient: s = s x y s y z s z Divergence: v = v x + v y + v z x y z Curl: v = v z v
More informationModeling External Compressible Flow
Tutorial 3. Modeling External Compressible Flow Introduction The purpose of this tutorial is to compute the turbulent flow past a transonic airfoil at a nonzero angle of attack. You will use the Spalart-Allmaras
More informationExamine 2D 8.0 is here!
Examine 2D 8.0 is here! Rocscience is pleased to announce the release of Examine 2D version 8.0, an upgrade to our very popular FREE stress analysis program. Examine 2D uses the boundary element method
More informationUnstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications
Unstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications Per-Olof Persson (persson@mit.edu) Department of Mathematics Massachusetts Institute of Technology http://www.mit.edu/
More informationGeoStudio Product Details
Product Details I. Purchasing Options Universal Professional Standard Vadose Basic Student (Free) Feature Comparison Finite Element Integration No limit 500 elements 500 elements Analysis Methods Ordinary
More informationCOMPUTATIONAL METHODS FOR ENVIRONMENTAL FLUID MECHANICS
COMPUTATIONAL METHODS FOR ENVIRONMENTAL FLUID MECHANICS Tayfun Tezduyar tezduyar@rice.edu Team for Advanced Flow Simulation and Modeling (T*AFSM) Mechanical Engineering and Materials Science Rice University
More informationDriven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
More informationStability with FE Stresses
1 Introduction Stability with FE Stresses In the finite element stresses stability analysis method, the stresses in the ground can be computed using SIGMA/W, and then SLOPE/W uses the SIGMA/W stresses
More informationSimulation of Turbulent Flow in an Asymmetric Diffuser
Simulation of Turbulent Flow in an Asymmetric Diffuser 1. Purpose 58:160 Intermediate Mechanics of Fluids CFD LAB 3 By Tao Xing and Fred Stern IIHR-Hydroscience & Engineering The University of Iowa C.
More informationFIDAP Fidap 8.5 Page 1 of 6
FIDAP 8.5 Page 1 of 6 FIDAP 8.5 is the CFD solver of choice for a wide variety of laminar and turbulent flows that arise in the polymer processing, thin film coating, biomedical, semiconductor crystal
More informationIMPROVING THE NUMERICAL ACCURACY OF HYDROTHERMAL RESERVOIR SIMULATIONS USING THE CIP SCHEME WITH THIRD-ORDER ACCURACY
PROCEEDINGS, Thirty-Seventh Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 30 - February 1, 2012 SGP-TR-194 IMPROVING THE NUMERICAL ACCURACY OF HYDROTHERMAL
More informationFree surface computations in mantle convection models
Free surface computations in mantle convection models Ian Rose Bruce Buffett Timo Heister Abstract Geodynamic simulations increasingly rely on simulations with a true free surface to investigate questions
More informationMore Animation Techniques
CS 231 More Animation Techniques So much more Animation Procedural animation Particle systems Free-form deformation Natural Phenomena 1 Procedural Animation Rule based animation that changes/evolves over
More informationFinite Volume Discretization on Irregular Voronoi Grids
Finite Volume Discretization on Irregular Voronoi Grids C.Huettig 1, W. Moore 1 1 Hampton University / National Institute of Aerospace Folie 1 The earth and its terrestrial neighbors NASA Colin Rose, Dorling
More informationHomogenization and numerical Upscaling. Unsaturated flow and two-phase flow
Homogenization and numerical Upscaling Unsaturated flow and two-phase flow Insa Neuweiler Institute of Hydromechanics, University of Stuttgart Outline Block 1: Introduction and Repetition Homogenization
More informationNumerical Study of Turbulent Flow over Backward-Facing Step with Different Turbulence Models
Numerical Study of Turbulent Flow over Backward-Facing Step with Different Turbulence Models D. G. Jehad *,a, G. A. Hashim b, A. K. Zarzoor c and C. S. Nor Azwadi d Department of Thermo-Fluids, Faculty
More informationAdvances in LS-DYNA Metal Forming (II)
Advances in LS-DYNA Metal Forming (II) Xinhai Zhu, Li Zhang & Yuzhong Xiao Livermore Software Technology Corporation Abstract Some of the new features developed since the last conference will be discussed.
More informationOpenSees2DPS (Beta 0.1): 2D Plane Strain Analysis
University of California, San Diego Department of Structural Engineering OpenSees2DPS (Beta 0.1): 2D Plane Strain Analysis User Manual Draft Ning Wang, Jinchi Lu, Zhaohui Yang and Ahmed Elgamal (Please
More informationIMAGE ANALYSIS DEDICATED TO POLYMER INJECTION MOLDING
Image Anal Stereol 2001;20:143-148 Original Research Paper IMAGE ANALYSIS DEDICATED TO POLYMER INJECTION MOLDING DAVID GARCIA 1, GUY COURBEBAISSE 2 AND MICHEL JOURLIN 3 1 European Polymer Institute (PEP),
More informationEulerian Techniques for Fluid-Structure Interactions - Part II: Applications
Published in Lecture Notes in Computational Science and Engineering Vol. 103, Proceedings of ENUMATH 2013, pp. 755-762, Springer, 2014 Eulerian Techniques for Fluid-Structure Interactions - Part II: Applications
More informationSimulation of Flow Development in a Pipe
Tutorial 4. Simulation of Flow Development in a Pipe Introduction The purpose of this tutorial is to illustrate the setup and solution of a 3D turbulent fluid flow in a pipe. The pipe networks are common
More informationSolidification using Smoothed Particle Hydrodynamics
Solidification using Smoothed Particle Hydrodynamics Master Thesis CA-3817512 Game and Media Technology Utrecht University, The Netherlands Supervisors: dr. ir. J.Egges dr. N.G. Pronost July, 2014 - 2
More informationThree Dimensional Numerical Simulation of Turbulent Flow Over Spillways
Three Dimensional Numerical Simulation of Turbulent Flow Over Spillways Latif Bouhadji ASL-AQFlow Inc., Sidney, British Columbia, Canada Email: lbouhadji@aslenv.com ABSTRACT Turbulent flows over a spillway
More informationThe Development of a Volume-of-Fluid Interface Tracking Method for Modeling Problems in Mantle Convection
The Development of a Volume-of-Fluid Interface Tracking Method for Modeling Problems in Mantle Convection Jonathan Robey 1 August 27, 2016 1 Funded in part by a GAANN Fellowship 1 Summary In problems involving
More informationA 3D VOF model in cylindrical coordinates
A 3D VOF model in cylindrical coordinates Marmar Mehrabadi and Markus Bussmann Department of Mechanical and Industrial Engineering, University of Toronto Recently, volume of fluid (VOF) methods have improved
More informationANSYS AIM Tutorial Turbulent Flow Over a Backward Facing Step
ANSYS AIM Tutorial Turbulent Flow Over a Backward Facing Step Author(s): Sebastian Vecchi, ANSYS Created using ANSYS AIM 18.1 Problem Specification Pre-Analysis & Start Up Governing Equation Start-Up Geometry
More informationNumerical Modelling of Engineering Structures subjected to Dynamic Seismic Loading
Paper 115 Numerical Modelling of Engineering Structures subjected to Dynamic Seismic Loading S.N. Polukoshko 1, V.F. Gonca 2 and E.V. Uspenska 2 1 Engineering Research Institute Ventspils International
More informationSimulation of Connector Assembly C
Simulation of Connector Assembly C Date: Sunday, March 6, 2016 Designer: Solidworks Study name: Horizontal Stress Test on C inner bend Analysis type: Static Table of Contents Model Information... 2 Study
More information3DEC 3DEC ITASCA ITASCA VERSION 5.0 VERSION 5.0
ITASCA Consulting Group, Inc. An Itasca International Company 3DEC VERSION 5.0 Advanced, Three Dimensional Discrete Element Modeling for Geotechnical Analysis of Rock, Blocky Structures and Structural
More informationDynamic CFD Modelling of Deploying Fins During Transitional Ballistic
Dynamic CFD Modelling of Deploying Fins During Transitional Ballistic Anton Jybrink Engineering Physics and Electrical Engineering, master's level 2018 Luleå University of Technology Department of Computer
More informationcuibm A GPU Accelerated Immersed Boundary Method
cuibm A GPU Accelerated Immersed Boundary Method S. K. Layton, A. Krishnan and L. A. Barba Corresponding author: labarba@bu.edu Department of Mechanical Engineering, Boston University, Boston, MA, 225,
More information