Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Current nondimensionalized time step of 1.0000e+00 exceeds the nondimensionalized stable time step of 1.5418e-01. #783

Open
yuxin0909 opened this issue Oct 29, 2024 · 4 comments
Labels
debugging-help Need help debugging

Comments

@yuxin0909
Copy link

Hi,I'm using the pylith 2.2.2 version and I'm having problems running the viscoelastic model“RuntimeError: Current nondimensionalized time step of 1.0000e+00 exceeds the nondimensionalized stable time step of 1.5418e-01.I have three viscoelastic blocks, two of which are 1e+20pas and one that is 1e+18pas, and then I get an error on run.

@yuxin0909 yuxin0909 added the debugging-help Need help debugging label Oct 29, 2024
@baagaard-usgs
Copy link
Contributor

In a quasi-static simulation with a Maxwell viscoelastic model, PyLith checks whether the time step is less than 0.2 * Maxwell time. Given your shear modulus and viscosity, PyLith computes that the nondimensional stable time step is 0.154, but your nondimesional time step is 1.0. This means you need to reduce your time step. Assuming you are using the default scales for nondimensionalization (default time scale is 1 year) means you need to reduce your time step from 1.0year to 0.15year or less. In general, we recommend a time step that is about 0.1 times the shortest Maxwell time to ensure that you sufficiently resolve the relaxation.

@yuxin0909
Copy link
Author

yuxin0909 commented Oct 30, 2024 via email

@yuxin0909
Copy link
Author

yuxin0909 commented Oct 30, 2024 via email

@yuxin0909
Copy link
Author

yuxin0909 commented Oct 30, 2024 via email

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
debugging-help Need help debugging
Projects
None yet
Development

No branches or pull requests

2 participants